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Résumé 


Un nouveau complexe de transfert de charge comprenant une liaison hydrogène entre la 
B-cyclodextrine avec les deux  accepteurs 7  tétracyanoéthylène (TCNE) et 
dichlorodicyanobenzoquinone (DDQ) a été étudié théoriquement par la méthode B3LYP/6-31G 
(d,p) et expérimentalement par les techniques RMN du 'H et IR-TF. L'étude par spectroscopies 
IR-TF et RMN du 'H a confirmé que la réaction a précédé la formation d'un complexe de 
transfert de charge (CTC). 

La stoechiométrie des complexes CT obtenus a été déterminée comme étant de 1:1 pour 
les deux complexes. Les distinctions remarquées dans les bandes vibrationnelles observées du 
complexe et de leurs constituants pour les groupes fonctionnels, c'est-à-dire C = O, O-H, C=N, 
C-H et C-CI participant à la formation des liaisons Hydrogène, affirment la formation du 
complexe actuel. 

Les complexes CT ont été formés par transfert de proton de B-CD vers DDQ et TCNE 
montrant des liaisons N*--H-O° et C=0*--H-0:. Comprendre l'interaction entre le donneur et 
l’accepteur capable d'interpréter l'interaction donneur-accepteur et le mécanisme d'action entre 
ces COMPOsÉS. 

L'énergie d'optimisation, l'énergie de complexation, les paramètres géométriques, les 
charges atomiques de Mullikan, les énergies LUMO et HOMO ainsi que l'optique non linéaire 
(NLO) ont été calculés et interprétés ; ils ont fortement référé la grande stabilité du complexe à 
l'existence d'un transfert de charge à côté de la liaison hydrogène dans le composant. 

En outre, les calculs de la théorie quantique des atomes dans les molécules (QTAIM) et 
de l'orbitale de liaison naturelle (NBO) ont été analysés et discutés ; ils démontrent clairement 


qu'un transfert de charge se produisait entre les molécules DDQ, TCNE et la B-CD. 


Mots clés: Complexe, f-Cyclodextrine, Transfert de charge, DDQ, TCNE, Transfert de proton, 
QTAIM, Liaisons hydrogène. 


Abstract 


A novel charge transfer complex including hydrogen bonding between f-cyclodextrin 
with the two r-acceptors tetracyanoethylene (TCNE) and dichlorodicyanobenzoquinone (DDQ) 
have been investigated theoretically by B3LYP/6-31G (d,p) method and experimentally by ‘H 
NMR and FT-IR techniques. FT-IR and 'HNMR spectroscopy investigation confirmed that the 
reaction precceded through charge transfer complex (CTC) formation. 

The stoichiometrie of the obtained CT complexes was determinate to be 1:1 for all two 
complexes. The distinctions noticed in the observed vibrational bands of the complex and their 
constituents for functional groups 1.e. C=0O, O-H, C=N, C-H and C-Cl participating in H-bond 
formation affirms the formation of present complex. 

CT complexes was formed by proton transfer from fB-CD to DDQ and TCNE showing 
N°--H-O° and C=O*--H-O° bondings. Understanding the interaction between donor and 
acceptor able to interpret the donor-acceptor interaction and acting mechanism between these 
compounds. 

The optimization energy, complexation energy, geometrical parameters, Mullikan 
atomic charges, LUMO and HOMO energies as well as non-linear optical (NLO) were 
calculated and interpreted; they strongly referred the high stability of the complex to the 
existence of charge transfer beside hydrogen bonding in the formed. 

In addition, Quantum Theory of Atoms in Molecules (QTAIM) and Natural Bonding 
Orbital (NBO) calculations were analysed and discussed; they clearly demonstrate that a charge 


transfer was occurring between DDQ, TCNE and B-CD molecules. 


Key words: Complex, B-Cyclodextrin, Charge transfer, DDQ, TCNE, Proton transfer, QTAIM, 
Hydrogen bonding. 
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La chimie est une science de la nature expérimentale qui étudie la composition de la 
matière et ses transformations. Elle joue un rôle primordial dans notre compréhension des 
phénomènes matériels, dans notre capacité d’agir sur eux, de les modifier, de les contrôler et 
d’inventer de nouvelles expressions [1]. 

A l'interface entre la chimie, la biologie et la physique, nous trouvons la chimie 
supramoléculaire, comme elle est définie par l’un de ses principaux fondateur, Jean-Marie 
Lehn, lauréat du prix Nobel pour son travail dans ce domaine en 1987 ; « la chimie au-delà de 
la molécule « ou » la chimie des assemblages moléculaires et de la liaison intermoléculaire ». 
Elle porte sur les entités organisées, d’une complexité supérieure, et qui résultent de 
l’association de deux ou plusieurs espèces chimiques maintenues ensemble par des interactions 
non covalentes [2]. 

Les interactions non covalentes jouent un rôle constitutif dans la science des relations 
intermoléculaires. Les processus d'assemblage chimique et biologique, orchestrés par une 
liaison non covalente, sont dirigés par des expressions élégantes d'un comportement collectif 
à l'échelle moléculaire [3]. Dans la nature, ces interactions constituent le fondement du 
processus de vie lui-même, l’articulation de la fonction ultime, à la fois mécanique et 
cognitive. En chimie de synthèse, les interactions entre sous-unités moléculaires conçues de 
manière rationnelle conditionnent l'assemblage d'agrégats nanoscopiques dotés de fonctions 
ciblées. La recherche dans ce domaine s’inspire de tout, des mécanismes fondamentaux de la 
fonction aux caractéristiques structurelles des systèmes biologiques [4-20]. 

Une compréhension claire et une description précise de la totalité des interactions 
entre les molécules sont essentielles pour le développement de la chimie supramoléculaire, en 
particulier l’élucidation des mécanismes des fonctions biologiques et le développement de 
nouvelles applications synthétiques dans la catalyse, la science des matériaux et la médecine 
[21]. 

De nombreuses classes de molécules peuvent former des complexes par des 
interactions non covalentes. Ces dernières années, l’accent a été mis sur les cyclodextrines en 
tant qu’agents complexants (à savoir des ligands) et en tant que solubilisant dans les produits 
pharmaceutiques. Chaque année, les cyclodextrines font l’objet de plus de 3000 articles de 


recherche et de plus de 1500 brevets [22]. 
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Une recherche bibliographique sur les complexes d’inclusion avec la B cyclodextrine 
(BCD) a fait remarquer que ces associations formées suite à l’insertion partielle ou totale de 
l’invité dans la molécule hôte sont stabilisées par des différentes interactions non covalentes 
où la (BCD) joue le rôle dans la majorité des cas de donneur par l’intermédiaire des atomes 
d'oxygène. Dans ce cadre, il serait alors intéressant, d’estimer le rôle de cette interaction via la 
formation des complexes de transfert de charge avec les accepteurs x. 

La complexation par transfert de charge revêt une grande importance 
dans les réactions chimiques, y compris l'addition, la substitution, la condensation, les 
systèmes biochimiques et biologiques [23-25]. Elle trouve également une large application dans 
nombreuses applications et de nombreux domaines, les divers domaines de matériaux optiques 
non linéaires, les matériaux électriquement conducteurs, la chimie de surface, les 
photocatalyseurs, le stockage d'énergie solaire, les semi-conducteurs organiques et l'étude de 
processus redox [26-31]. 

Ce travail se propose de mener une étude physico-chimique de l’interaction entre la f- 
CD en tant que donneur et deux composés, à savoir, la 2,3-dichloro-5,6-dicyano-1,4- 
benzoquinone (DDQ) et le tétracyanoéthylène (TCNE) en tant accepteurs x. 

Dans la première partie menée expérimentalement, les complexes à l'état solide ont été 
synthétisés et caractérisés par la spectroscopie infrarouge par transformée de Fourier IR-TF, et 
par la technique de la résonance magnétique du proton RMN'H. 

La seconde partie a été consacrée à une étude par modélisation moléculaire au niveau 
DFT et à travers laquelle des informations complémentaires ont été obtenues : structure des 
complexes, propriétés électroniques et orbitales moléculaires. 

La nature ayant joué un rôle dans la stabilisation de ces complexes a été étudiée en 
menant les analyses des orbitales moléculaire naturelle (NBO), théorie des atomes dans les 


molécules (QTAIM) et l’approche du gradient de densité (RDG). 


Ce manuscrit s’articule en deux grandes parties : 
La première partie «Synthèse Bibliographique» de ce manuscrit se compose de trois chapitres. 


Un bref rappel bibliographique sur les complexes moléculaire, fait l’objet du premier chapitre. 
Des généralités sur le DDQ, le TCNE, les interactions non covalentes et le transfert de proton 
constituent le deuxième chapitre. Le troisième chapitre est consacré, à de brefs rappels des 
différents concepts et outils théoriques que nous avons utilisé. 

La deuxième partie intitulée «Calculs, Résultats et Discussions» présente les études théoriques 


et expérimentales proprement dites menées dans ce travail et la discussion des résultats obtenus. 


', 
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Elle comporte les chapitres quatre et cinq : Dans le quatrième chapitre sont décrites en détail, 
la synthèse et la caractérisation des complexes de transfert de charge entre la B-cyclodextrine 
et les accepteurs (DDQ et TCNE). Le cinquième chapitre rapporte une étude théorique par 
modélisation moléculaire de deux complexes : [(B-CD) (DDQ)] et [(B-CD) (TCNE)]. 


Enfin, nous clôturons notre manuscrit par une Conclusion Générale et une annexe. 
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COMPLEXES MOLECULAIRES : 


Inclusion et transfert de charge 


Chapitre I Complexes Moléculaires : Inclusion et Transfert de charge 
I. Les complexes moléculaires 


La complexation moléculaire impliquant deux molécules neutres. Il peut s'agir de 
l'inclusion si l'hôte est une macromolécule pouvant inclure partiellement ou totalement la 
molécule de l'invité, ou de transfert de charge si l'une des molécules est donneuse de charge et 


l'autre accepteuse. 


L.1 Les cyclodextrines 
L.1.1 Introduction 


La discipline de la chimie supramoléculaire fit son apparition dans les années 70, basée 
sur l'association entre deux ou plusieurs espèces chimiques, elle est édifiée non pas à l'aide de 
haisons covalentes (liaisons intramoléculaires connectant les atomes entre eux pour former des 
molécules) mais, par des interactions intermoléculaires. La majorité de ces interactions sont du 
type hôte-invité (host-guest) ou complexe d'inclusion. L'un des premiers chimistes qui s'orienta 
vers ces interactions intermoléculaires fut Jean-Marie Lehn [1], qui réussit en 1969 à synthétiser 
des molécules tridimensionnelles creuses capables d'enfermer des cations à l'intérieur d'une 
cavité « crypte », qui sont constituées d'atomes d'azote de carbone et d'oxygène. 

En 1987, D.J.Cram, J.M Lehn et C.J.Pedersen reçurent le prix Nobel de chimie pour 
leurs travaux sur le développement et les applications des molécules qui peuvent se reconnaitre 
entre-elles et choisir avec quelles autres molécules elles formeront un complexe. Les trois 
chercheurs travaillaient respectivement sur les Cryptands et les éthers couronnes. 

Les complexes sont formés à partir d’une espèce réceptrice (hôte) qui inclut plus ou moins 
profondément un substrat moléculaire ou ionique (invitée), de ce fait une ou plusieurs molécules 
peuvent être « encapsulées » dans une ou deux molécules de cyclodextrine (CD). Cette 
inclusion n’est, en général, pas accompagnée de la formation d’une liaison covalente ou de 
coordination, car leurs énergies sont en général bien plus faibles que celles des liaisons 
covalentes, les édifices intermoléculaires peuvent se faire et se défaire à une vitesse dépendant 
de leur nature, de leur nombre, et de leur disposition. Aïnsi, les processus d’interactions et de 


reconnaissances moléculaires sont gouvernés uniquement par des forces dites faibles. 


L.1.2 Historique 


Les cyclodextrines ont été isolées pour la première fois par Villiers [2]en 1891, grâce à 
l’expérience de la dégradation de l’amidon par une souche de micro-organismes (l’amylase de 


Bacillus macerans : cyclodextrinase). Villiers met en évidence deux produits (probablement 
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l’a- et la B-cyclodextrine) ayant des propriétés physico-chimiques proches de celles de la 
cellulose. Les cyclodextrines ont été caractérisées en 1903 par Schardinger [3] en tant que 
oligosaccharides cycliques ; c’est pour cette raison qu’elles sont nommées dextrines de 
Schardinger dans les premières publications traitant des cyclodextrines. En 1938 Freudenberg 
et al. [4,5] ont mis en évidence que les cyclodextrines sont construites à partir d’unités de D- 
glucoses liées entre elles par des liaisons à (1—4) glucosidiques. Freudenberg et al, ont 
découvert que les cyclodextrines étaient capables de former des complexes d’inclusion, et 
déterminent entièrement la structure de la y-cyclodextrine [6]. Dans les années 1950, le groupe 
de French [7] et celui de Cramer [8] ont intensément travaillé sur la synthèse et la purification 
de complexes de cyclodextrines. Le premier brevet concernant l’application des cyclodextrines 
pour la mise en forme d’un composé à activité biologique sera déposé par Freudenberg [9] en 
1953. A partir de ce moment, l’étude des cyclodextrines prend un essor considérable : 
fabrications industrielles, synthèse de cyclodextrines modifiées, synthèse de complexes 
d’inclusion, etc. Dans les années 1970-80, Szejtli [10,11], appelé aussi “godfather” des 
cyclodextrines apporte une importante contribution dans le domaine. Depuis 1970, on 
dénombre un peu plus de 130000 documents traitant des cyclodextrines (publications, brevets, 


abstracts). 


L.1.3 Propriétés des cyclodextrines 
L.1.3.1 Propriétés structurales 


Les CDs sont des oligosaccharides cycliques constitués d’unités D-glucosidiques 
obtenus par la dégradation enzymatique de l’amidon, plus particulièrement de sa composante 
linéaire, l’amylase, via l’enzyme CGTase. Cela permet l’obtention des a-, B- et y-CDs, aussi 
appelées CDs natives, dont les lettres grecques désignent leur nombre d’unités glucosiques : & 


pour6, B pour 7 et y pour 8 (Figure L.1) 
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CH0H ae, CH,0H 2 CH0H 


Amylose 


foyer 


Amidon 


CDs natives 


Figure L.1 CDs obtenues de la dégradation enzymatique de l’amidon par la 


CGTase. 


Les CDs sont composés d’unités glucopyranoses sous une conformation chaise [12] C1, 
liées entre elles par des liaisons à (1—4). C’est cette conformation qui conduit à une structure 
de type cône tronqué de l’anneau oligosaccharidique. Tous les groupements hydroxyles 
secondaires (liés aux C2 et C3) sont placés sur la face la plus large du cône, appelée face 
secondaire. Les hydroxyles primaires (liés aux C6) se trouvent quant à eux sur la face la plus 
étroite, la face primaire (Figure 1.2). Etant donné que les fonctions hydroxyles sont placées à 
l’extérieur de l’anneau, elles rendent hydrophile la surface externe des CDs qui, par conséquent, 
sont solubles dans l’eau. Pour sa part, l’intérieur du cône comprend les liaisons C-H 
(notamment des C3, C5 et C6) et les —-O glycosidiques, rendant la cavité de la CD plus 
hydrophobe (Figure L.2) [13]. 
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Figure L.2 Structure caractéristique des CDs. 


De plus, le groupement hydroxyle du C2 d’une unité glucose peut établir une liaison 
hydrogène avec le groupement hydroxyle du C3 du glucose adjacent (Figure [.3). Cela permet 
donc de stabiliser la molécule et l’influence de cette liaison hydrogène sera discutée dans le 


paragraphe suivant. 


Figure L.3 Représentation de la liaison hydrogène intramoléculaire entre le 
groupement hydroxyle du C2 d’un glucose et le groupement hydroxyle du C3 


du glucose adjacent d’une CD. 


Concernant les trois CDs natives, l’a-CD présente une ceinture de liaisons hydrogène 
incomplète puisqu’une de ses unités glucoses est dans une forme distordue. Ainsi, à la place de 
former 6 liaisons hydrogène, seulement 4 sont simultanément possibles. Dans le cas de la f- 
CD, cette ceinture de 7 liaisons hydrogène lui donne une structure rigide qui empêche son 
hydratation et diminue ainsi sa solubilité dans l’eau. En effet, la B-CD est la moins soluble des 
CDs natives (18,5 g/L contre 145 g/L pour l’a-CD). La y-CD, pour sa part, dispose d’une 
structure non coplanaire et plus flexible, se présentant ainsi comme la plus soluble des CDs 


natives (232 g/L contre 145 g/L pour l’a-CD) mais aussi la moins stable chimiquement [10]. 


Chapitre I 


L.1.3.2 Propriétés physico-chimiques 


Complexes Moléculaires : Inclusion et Transfert de charge 


Les principales propriétés physico-chimiques des cyclodextrines naturelles a, B et y 


sont regroupées dans le tableau L.1. 


Tableau [L.1 Propriétés physico-chimiques des cyclodextrines naturelles 


a-CD B-CD y-CD 
Nombre d’unités glucose 6 7 8 
Formule brute C36H60030 C42H70035 C48H80040 
Masse molaire (g/mol) 972 1135 1297 
Solubilité dans l’eau à 
25°C(g/L) 145 18,5 232 
Diamètre de la cavité (À) 4,7-5,3 6,0-6,5 7,5-8,3 
Diamètre de la périphérie 14,6 +0,4 15,4 +0,4 17,5 +0,4 
(A) 
Hauteur du cône (À) 7,9 +0,1 7,9 +0,1 7,9 +0,1 
Volume de la cavité (À?) 174 262 427 
Nombre de molécules 6-8 11-12 13-17 
d’eau dans la cavité 
Pouvoir rotatoire [œ]p à + 150,5 + 162,5 + 177,4 
25°C 


L.1.4 Domaines d’application des cyclodextrines 


La propriété d’inclusion des cyclodextrines, mise en évidence dans les années 1930 et 
largement acceptée à partir des années 1950, est à l’origine de la plupart des applications 
industrielles. Grâce à cette propriété, les cyclodextrines sont impliquées dans différents 
domaines industriels, tels que la pharmacie, la chimie analytique, l’environnement, la 
cosmétique, le textile, l’alimentation, etc. ... 

Dans le domaine pharmaceutique les cyclodextrines sont généralement utilisées pour 
améliorer la solubilité de médicaments [14]. En raison de leur capacité à changer les propriétés 
des molécules invitées par la formation d’un complexe d’inclusion, elles agissent généralement 


en tant que système de transport de molécules bioactives à travers les membranes biologiques. 
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On les trouve aussi pour le masquage des effets secondaires, ainsi que le stockage et 
l’absorption du médicament, les cyclodextrines étant dégradables par une enzyme, l’a-amylase, 
provenant de micro-organismes de la flore intestinale [15]. La plupart des médicaments à base 
de cyclodextrines sont administrés par voie orale (tablettes, dragées, sirops. . . ). Elles peuvent 
encore éliminer une odeur ou un goût désagréable. 

Elles jouent un grand rôle dans les sciences environnementales, notamment pour la 
solubilisation de contaminants organiques, pour l’enrichissement et le déplacement des 
polluants organique et des métaux lourds du sol de l’eau et de l’atmosphère. Les cyclodextrines 
sont utilisés dans le traitement de l’eau afin d’éliminer les agents contaminants comme le 
phénol, le p-chlorophénol ou le benzène [14]. Elles sont également présentes dans la 
formulation de certains pesticides [14]. 

Les cyclodextrines ont trouvé de nombreuses applications dans le domaine alimentaire, 
elles forment des complexes d’inclusion avec une variété de molécules comprenant des arômes, 
des graisses et des colorants. Elles sont utilisées pour masquer des 
composants indésirables, ou relarguer des constituants désirés tels que les arômes [16]. Dans 
l’industrie cosmétique, ses principaux avantages sont la stabilisation et le contrôle des odeurs, 
la diminution de la volatilité des parfums et les procédés permettant d’augmenter la conversion 
d’une substance liquide à sa forme solide par précipitation du complexe d’inclusion. On peut 
les trouver dans les dentifrices, les crèmes corporelles, les adoucissants. . . [14] 

Le « finissage » du textile est aussi un secteur où les cyclodextrines attirent l’attention. 
Elles peuvent être directement greffées sur les textiles dans le but toujours d’éliminer ou 
masquer des odeurs gênantes, comme la fumée de cigarette. Elles permettent aussi de renforcer 
la prise de colorants par le tissu afin de réduire la quantité de colorant perdue au lavage [17]. 


Les cyclodextrines sont également très utilisées dans les techniques séparatives 
(HPLC, CPG) ou elles servent d’additifs dans les phases stationnaires ou mobiles pour 
améliorer les séparations [18]. 


L.2 Complexe d’inclusion 
L2.1 Introduction 


La propriété la plus caractéristique des cyclodextrines est leur aptitude à former des 
complexes d’inclusion en solution aqueuse avec des molécules dont la taille et l’hydrophobicité 
correspond à celle de leur cavité [19]. Cette propriété a été mise en évidence dans les années 


1930 [20,21]. Le terme « complexe d’inclusion » a été présenté en 1950 [8]. 


Un complexe d’inclusion est une association moléculaire entre une ou plusieurs 
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molécules dont l'une est l'hôte (le récepteur) et l'autre l'invitée (le substrat), sous l’effet 
d'interactions faibles [22]. Donc, aucune liaison covalente n’est créée, ce qui permet une 
dissociation aisée du complexe formé. 

La molécule « invitée » est alors encapsulée de façon totale ou partielle, le récepteur 
jouant le rôle de molécule « hôte ». Les tailles respectives de la cyclodextrine et du substrat 


conditionnent généralement la stæœchiométrie du complexe [23,24]. 


L2.2 Principe de formation du complexe d’inclusion 


Les cyclodextrines peuvent donc se comporter comme des molécules « refuges » 
capables d'inviter dans leur espace interne, de façon réversible, des molécules et donner ainsi 


des complexes d’inclusion. 


Le phénomène de complexation est la résultante d'une multitude d'interactions 
(substrat/solvant, solvant/solvant et CD/solvant), qui conduisent à l'état thermodynamique le 
plus stable. 

La cavité apolaire de la CD est occupée par des molécules d’eau, énergétiquement 
défavorables (interactions polaire / apolaire). Ces molécules d’eau pourront donc être 
facilement substituées par une “molécule invitée” appropriée, moins polaire que l’eau (Figure 


L4). 


CD (Hôte) Complexe (hôte-invité) 


Figure L.4 Représentation schématique de la formation d’un complexe d’inclusion. 


La formation d’un complexe est un processus réversible gouverné par un équilibre 
thermodynamique entre l’invité libre, l’hôte libre et le complexe sous forme : 


CD+I ++ CD (1) 


Chapitre Complexes Moléculaires : Inclusion et Transfert de charge 


La constante d’association Ka est définie par l’équation (1.2), qui est le rapport entre la 
concentration du complexe CD/I et le produit des concentrations de la CD et de la molécule 
invitée I : 


Ka = [CD-I]/[CD] [I] (2) 


L.2.3 Stæchiométrie du complexe d’inclusion 
Les complexes formés entre la (les) cyclodextrine(s) “hôte(s)” et la (les) molécule(s) 
“invitée(s)” peuvent être de plusieurs types [25]: on distingue les complexes d’inclusion (1-6) 


(hôte : invité) et les complexes d’association (7). 


La Figure IS résume les nombreux exemples de complexes CD-invité, avec divers 


arrangements structuraux. 


4 (1:22) 


5 (2:1) 6 (2:2) 7 


Figure L.S Représentation schématique des différents types du complexe d’inclusion de 


cyclodextrine. 


L2.4 Forces impliquées dans les complexes d’inclusion des cyclodextrines 


Le phénomène de complexation est la résultante de nombreuses interactions mettant en 
jeu la molécule hôte, la molécule invitée et le solvant et conduisant à l’état thermodynamique 
le plus stable. La nature et l’importance relative des différents phénomènes responsables de la 
formation d’un complexe CD-soluté sont encore mal connues et font toujours l’objet d’une 


controverse. Les forces motrices les plus discutées sont les suivantes [26]: 


1. Interactions électrostatiques (ion-ion, ion-dipôle, dipôle-dipôle). 


2. Interactions de van der Waals (dipôle-dipôle induit, dipôle induit-dipôle induit). 
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3. Interactions hydrophobes. 

4. Liaisons hydrogènes (qui sont en fait de nature électrostatique). 

5. Soulagement des contraintes de conformation de la CD (tension). 

6. Exclusion des molécules d’eau initialement présentes dans la cavité. 


7. Interactions de transfert de charges (type d’interaction de van der Waals). 


L.2.5 Conséquences de la complexation et applications 


L’inclusion des molécules invitées dans la cavité des CDs constitue une encapsulation 
moléculaire et ces molécules voient leurs propriétés modifiées. Parmi ces modifications, on 


peut citer [27, 28]: 


(1) modification de l’activité chimique de la molécule invitée 
Ÿ” certaines substances réactives sont protégées par l’inclusion et peuvent ainsi être 
mélangées à d’autres substances à moindre risque, 
Ÿ” des réactions chimiques peuvent être réalisées sélectivement par inclusion de certains 
groupements fonctionnels (les CDs jouent alors le rôle de catalyseurs), 
Ÿ” certaines réactions peuvent être favorisées ou, à l’opposé, supprimées. 
(2) modification des propriétés physico-chimiques de la molécule invitée 
Y_ des substances peu solubles dans l’eau voient leur solubilité s’accroître par addition de 
CDs, 
Ÿ_ la couleur de certaines substances peut être altérée puisque, généralement, l’inclusion 
induit un changement des propriétés spectrales de l’invité, 
Ÿ” certains goûts déplaisants peuvent être éliminés. 
(3) stabilisation de substances sensibles à la lumière ou à l’oxygène 
(4) fixation de substances très volatiles 
Ÿ”_ le stockage et le maniement de certaines substances toxiques peuvent être améliorés, 
Ÿ”_ la quantité de substance volatile peut être économisée puisque l’évaporation est réduite, 
Y les quantités de substances aromatiques et actives physiologiquement peuvent être 


facilement quantifiées. 


L3 Complexes de transfert de charge 


Un complexe de transfert de charge (CTC) ou complexe électron-donneur-accepteur est 
une association de deux ou plusieurs molécules, ou de différentes parties d'une grosse molécule, 
dans lequel une fraction de la charge électronique est transférée entre les entités moléculaires. 


Le résultat électrostatique d’attraction fournit une force de stabilisation pour le complexe 
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moléculaire. La molécule de la source à partir de laquelle la charge est transférée est appelée le 
donneur d'électrons et l’espèce de réception est appelé accepteur d'électrons. La nature de 
l'attraction dans un complexe de transfert de charge n'est pas une liaison chimique stable, et est 
donc beaucoup plus faible que les forces covalentes. Beaucoup de ces complexes peuvent subir 
une transition électronique dans un état électronique excité. L'énergie d'excitation de cette 
transition se produit très fréquemment dans la région visible du spectre électromagnétique, qui 
produit la couleur intense caractéristique de ces complexes. 

Le terme complexe de transfert de charge (CTC) a été introduit par Mulliken [29,30]. 
Pour définir un nouveau type de produit d’addition pour expliquer le comportement de certaines 
classes de molécules qui ne sont pas conformes aux modèles classiques de type ioniques, 
covalentes, et la coordination des composantes de liaison hydrogène, bien que ces produits 
d’addition conservent largement une partie des propriétés des composantes. Certains 
changement sont évidents, par exemple: solubilité, diamagnétique et le paramagnétique. 
D’autres différences ont également été trouvées par des techniques électrochimiques. 

Mulliken [31,32] a également montré que les interactions de transfert de charge au sein 
d’un complexe moléculaire constitué d’un donneur d’électrons; D et un accepteur d’électrons; 


À, cause d’une résonance avec un transfert de charge de D à A: 
D + À — D*A°-— D'+ A (13) 


Les complexes de transfert de charge (CTC) qui proviennent d’une faible interaction 
entre donneur d’électrons et molécules accepteurs en équilibre biomoléculaire ou dans des 
composés modèles interaction intramoléculaire. 

En général, le transfert de charge complet conduit à une liaison ionique dans une interaction 
simple paire d’ion radical [33]. 

Un complexe de transfert de charge résulte d’un rapport partiel d'un électron à partir 
d'une molécule donneur ou, plus rarement, d’une entité atomique à une acceptation. Un tel 
transfert partiel ne doit pas être interprété comme signifiant qu'il est le seul, ni même le mode 
dominant par lequel le complexe est stabilisé. Le transfert de charge partielle, 1ci, signifie que 
le transfert d'électrons passe relativement plus de temps dans le voisinage de la molécule 


d'accepteur (A) que dans celui de la molécule donneuse (D). 


L3.1 Développement historique de la théorie de transfert de charge 


Le travail expérimental classique qui a mené à la théorie moderne du complexe de 


transfert de charge est celui de Benesi et de Hildebrand [34,35] qui ont étudié l'effet de divers 
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solvants sur les spectres d'absorption de l'iode moléculaire. Ils ont constaté qu'un mélange 
d'hydrocarbures aromatiques (benzène) et d'iode possède des maxima d'absorption ne figurant 
pas dans les spectres de benzène ou l'iode. Ils ont attribué cette nouvelle bande à la formation 
d'un produit d'addition entre les deux composés, et ont commencé à examiner la nature de ce 
nouveau produit appelé complexe, par modification des substituants sur le benzène. Ces 
modifications spectrales sont le résultat de l'addition d'électrons avec la libération des groupes 
au benzène. Ils ont conclu que ces complexes sont le résultat d'une interaction acide-base au 
sens de Lewis. Cette étude expérimentale systématique, fournie par Mulliken [36] a permis de 
rer une théorie acceptable de ce phénomène. 

Il est noté depuis longtemps que les explications théoriques des complexes sont fondées 
sur les descriptions de valence. On comprit alors que de nouvelles propriétés d'absorption 
qu'une nouvelle entité chimique a été présente et la seule façon preuve de la formation de la 
nouvelle substance chimique serait à travers de la formation d’une liaison. Il n'est donc pas 
surprenant que les premières théories de Sudborough [37] fassent la suggestion selon laquelle 
les complexes sont de nature covalente. L’idée que la formation du complexe a été liée à la 
formation de liaisons covalentes a persisté pendant trente ans. La structure a été d’abord 
expliquée en termes de modèle de l'atome de Bohr par Sidgwick [38] qui cherchait à tirer un 
ensemble de règles de valence générales. Dans la même année le sujet de valence a été abordé 
dans une réunion de la Société de Faraday, dans un article présenté par Fowler [39], la nature 
covalente de composés de coordination a été soulignée et 1l a été suggéré que la complexation 
est le résultat de la mise en commun ou l'association d'orbitales (à ce moment-là la mécanique 
quantique était appliquée aux problèmes de structure atomique). 

La supposition de la formation de liaison apparaît dans les discussions ultérieures où 
Fowler [40] soutient que la formation du complexe est impossible dans le cas où un seul électron 
peut être partagé; et ce, malgré une remarque contradictoire de Bragg [41] qui, sur la base 
d'études de diffraction des rayons X, a indiqué qu'il existe des forces entre les molécules 
organiques beaucoup plus faibles que les forces de valence. Il a noté que les distances entre les 
molécules complexées étaient plus importantes que celles entre les atomes dans les molécules 
et il a conclu que les forces entre les molécules complexantes étaient plus faibles. 

Pfeiffer [42] de sa part a indiqué que la complexation est le résultat de la force résiduelle de 
valence. La théorie de Pfeiffer a été considérée comme inacceptable par Bennett et Willis [43] 
qui ont estimé que la valence résiduelle ne pourrait pas expliquer la stœchiométrie de nombreux 
complexes ni leur stabilité en solution, et 1ls ont donné comme exemple le complexe entre le 


benzène et des composés nitrés : 
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+ __-0 
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No 


Figure L.6 Complexe entre le benzène et des composés nitrés. 


La vue la plus extrême jugeant les forces électrostatiques responsables de la formation 
du complexe était celle de Weiss [44] qui a considéré que les formes ioniques des molécules 


comme parties composées du complexe. Il a proposé : 
A +B — [AB]: — A*B: (1.4) 


Weiss a proposé que la fonction d'onde totale du complexe puisse être écrite comme une 
combinaison linéaire des fonctions d'onde qui décrivent le système de van der Waals (forces de 
dispersion de LONDON) et des interactions ioniques. XJN 

La détermination du coefficient de pondération est définie comme la différence I — £, 
où I est le potentiel d'ionisation et &Æ l'affinité électronique du donneur et de l’accepteur, 
respectivement. Dans le cas ou I- est petit, le complexe formé est de caractère ionique, ces 
idées ont été justifié plus tard par Benesi et Hildebrand. 

Le développement de la théorie des complexes de transfert de charge par la suite s’est 


grâce aux études de Benesi et Hildebrand. 


L.3.2 Théorie de Mulliken (théorie donneur-accepteur) 


Mulliken [45,46] s’est consacré à l’étude des complexes par transfert de charge, sa 
théorie s’applique au transfert électronique complet ou partiel d’un partenaire à l’autre du 
complexe de TC qui se traduit généralement par l’apparition d’une nouvelle bande d’absorption 
caractéristique en spectroscopie ultraviolette. La théorie de Mulliken s’applique seulement aux 
cas ou le donneur et l’accepteur sont des molécules neutres aux couches saturées. 

Mulliken a exprimé la théorie de transfert de charge en termes de mécanique quantique 


[36, 47,48]. Selon cette théorie, l’interaction d’un état fondamental non lié W, (DA) et d’un état 
excité polarisé Ÿ, (D* A’) conduirait à : 


> Un état fondamental stabilisé dont la fonction d’onde Y est donnée par : 
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Y, =Y,(DA)+XY, (D'A°) (LS) 
> Un état excité de TC dont la fonction d’onde est : 
Y, =, (D'A )+uY,(DA) (L.6) 
Les coefficients À et LL sont en général négligeables devant l’unité. 
Il considère que la bande de TC du complexe est associée à la transition électronique Y, — Y, 
. Ce calcul peut être généralisé et on peut inclure un troisième terme correspondant à la structure 
(D'A°) dont la fonction d’onde est, . 

Une autre terminologie, utilisant les coefficients a, b et c à la place de X et u, a également 
été utilisée pour exprimer la théorie de transfert de charge en termes de la mécanique quantique. 
Dans ce cas, lorsque le donneur D (base faible) et l’accepteur A (acide faible) sont dans leurs 
états fondamentaux symétriques, la fonction d’onde du complexe moléculaire Y,,(DA) 
s’exprime par l’équation : 

YF, (DA) = aF,(DA)+bY,(D'A )+c#,(D A°) (1.7) 

"la fonction d’onde non liante Ÿ, correspond aux structures dans lesquelles la liaison 
est dûe aux forces intermoléculaires classiques (dipôle-dipôle, dipôle-dipôle induit, etc) ou aux 
dispersions de London. 

"La fonction Y, se rapporte à la liaison semi-polaire qui existe quand un électron 
est transféré du donneur D à l’accepteur À, de sorte, qu’en plus des forces classiques reliées à 
,; 

"La fonction d’onde Y, correspond à la structure DA”. 


" Les coefficients a, b, c, sont tels que c << b et que le rapport b°/a° est généralement 
très petit. 
L'énergie E de l’état fondamental associée à la fonction d’onde totale du complexe 
est : 
(W, -EXW, -E) = (H,, -ES) (L.8) 
Wo: énergie associée à la structure (DA) et elle est donnée par l’intégrale : 
W, = ['PHçdr (L9) 
Wii: énergie associée à la structure (D*A° ) [49]. 
W,= [F,HY,dr (L.10) 


H : Hamiltonien d’un complexe formé d’un nombre entier les noyaux et les électrons 


Chapitre I Complexes Moléculaires : Inclusion et Transfert de charge 
Soi : recouvrement de fonctions correspondantes aux structures (DA) et (DA). 
Su = [dr (1.11) 
H,, : Energie d’interaction de (DA) et (DA ). 


Ho:= S Y1H Wodt (12) 


L.3.3 Liaison et énergie de formation de liaison dans les complexes de TC 


Les effets de transfert de charge sont proches des effets d’association. Ils sont dûs à la 
déformation des nuages électroniques par recouvrement lorsque deux molécules arrivent en 
proche contact, si une molécule du couple possède un faible potentiel d’ionisation tandis que 
l’autre a une forte affinité électronique ; il en résulte la formation d’un CTC entre le donneur 


d’électrons du couple et l’accepteur : 


D+A— D", A*] | 


Le transfert partiel d’électrons du donneur vers l’accepteur modifie les paramètres 


structuraux (longueurs et angles de liaison) et les charges des couples. L’intensité de l’effet de 
stabilisation supplémentaire dépend du recouvrement entre les orbitales du donneur et de 
l’accepteur. L’énergie d’interaction qui en résulte est intermédiaire entre celle des liaisons de 
van der Waals et celle des liaisons covalentes faibles. Elle augmente lorsque les orientations 
respectives des molécules permettent le plus grand recouvrement, ce qui confère à ces liaisons 
un rôle important dans la détermination de la structure cristalline pour certains cristaux 
moléculaires. 

Les distances séparant les donneurs et les accepteurs dans les complexes d’addition à 
l’état cristallin, sont très grandes devant celles correspondant aux liaisons de covalence [501]. 
Elles sont aussi beaucoup plus faibles que celles attendues pour des liaisons de van der Waals. 

Dans la littérature, il est admis que la formation d’un complexe moléculaire est liée à 
l’existence d’une liaison covalente entre le donneur et l’accepteur [51]; les progrès 
expérimentaux ont permis de rejeter cette hypothèse et d’élaborer différentes théories nouvelles. 

Weiss a raffiné la théorie cinétique du déplacement électronique accompagnant les 
collisions intermoléculaires en faisant l’hypothèse du transfert d’un électron du donneur à 
l’accepteur. Il admet ainsi la formation d’une liaison semi-polaire intermoléculaire schématisée 


de la façon suivante : 
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D:+A—[D.| [A 
A la lumière de ces considérations, les couleurs intenses des divers complexes 
d’addition sont attribuées aux molécules possédant des électrons célibataires, celles-ci 
nécessitent, en général, de faibles énergies d’excitation [52,53]. Cependant, les enthalpies de 
formation des complexes nitrés ou quinoniques sont beaucoup trop faibles pour qu’il y ait 
formation de sels [54]. Ceci a amené Brackmann [55] à considérer la structure des complexes 
comme un mélange de deux formes hybrides en résonance, dans lesquelles 11 ya contribution 


des deux liaisons : semi-polaire et ionique, par exemple : 


 D:A&D*':A". 


L.3.4 Propriétés des complexes de TC 


On peut considérer que la transition de transfert de charge W,, — Y', oblige un électron 


à « sauter » de la molécule du donneur D à la molécule de l’accepteur A. Cette transition est 
alors caractéristique du complexe, car elle concerne à la fois D et A et se traduit par un spectre 
électronique intense. On observe en général, une modification importante de la forme des 


bandes du donneur et de l’accepteur et des variations considérables des fréquences de vibration. 


Ceci se produit même dans le cas de faibles complexes d’addition. Le moment (l,,) associé 


à cette transition de transfert de charge dépend des distances moyennes entre les électrons des 
orbitales des molécules donneur et accepteur, d’une origine arbitraire et la position moyenne 
d’un électron dont la distribution de charge est due au recouvrement des fonctions d’onde des 


deux espèces constituantes. 
Le moment dipolaire (u,,) d’un complexe moléculaire est donné par la relation : 
bu =—ef EDP rdr (1-13) 


Avec e la charge électronique et ri la distance vectorielle de l’électron i mesurée à partir d’une 
origine convenable. 
Le spectre électronique du complexe de transfert de charge est considérer sous quatre aspects : 
polarisation, intensité, énergie et fréquence. 
a) Polarisation : l’absorption de lumière polarisée dans une direction perpendiculaire aux 
plans benzéniques est plus intense que celle polarisée dans le plan ; la position du 
maximum d’absorption des complexes formés est déplacée vers les basses fréquences 


d’une quantité plus grande que lorsque la lumière est polarisée dans le plan [56]. En 
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effet, dans le processus de transfert de charge le déplacement électronique du donneur 
vers l’accepteur ne peut se produire que lorsque la composante de la lumière vibre 
perpendiculairement aux noyaux benzéniques [48]. 

b) Intensité : la force théorique de l’oscillateur peut être évaluée au moyen de l’équation : 


fi=4,704x10 7%, x LEN (1-14) 


Où un est le moment dipolaire associé à la transition N —Eet v.. est le nombre d’onde au 


maximum d’absorption [57]. D’un autre côté la force de l’oscillateur observé expéripérimen- 
-talement est donnée par l’expression: 

fo =L35x10 Ex (Vrux — Vip) (1-15) 
Où v,,, est la largeur de la bande d’absorption correspondant à la moitié de l’absorbance et &. 
est le coefficient d’extinction molaire au maximum d’absorption. 

c) Energie : l’énergie d’une transition de transfert de charge est donnée par la relation de 
mécanique quantique : 

E=hv=W -W, (1-16) 
Où les termes W se rapportent aux énergies associées aux structures correspondantes. 

d) Fréquence : les spectres visibles et ultraviolets des complexes de TC concernent des 
produits organiques ont été étudiés dans l’hypothèse de Mulliken [47], à savoir que le 
transfert de charge photochimique est responsable de l’absorption de la lumière [58,59]. 
La fréquence v.. d’un photon nécessaire pour effectuer un transfert de charge est 
donnée approximativement par la relation : 

vcr= Ib- #-A (1-17) 
Æ: Affinité électronique de l’accepteur. 
I, : Potentiel d’ionisation du donneur. 
À : constante, c’est la différence entre les énergies de liaisons des composés dans l’état 
fondamental et excité. 

On peut aussi calculer la fréquence v.. de transfert de charge selon la théorie de 
Mulliken et Person [45,46], qui donnent lieu à une corrélation avec les potentiels d’ionisation 
Pl, des donneurs pour les complexes d’un même accepteur, ou avec les affinités électroniques 
À des accepteurs lorsqu’on travaille sur un seul donneur : 

vcr= 1Db-%+ (Go- G1) + (Bo° + B1° / Ib — %# + (Go- Gi) (1-18) 


G, : Energie des interactions non-liantes entre A et D. 
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G, : Energie d’attraction coulombienne entre A° et D*. 


Les intégrales B; ayant la dimension d’une énergie, peuvent être explicités par un calcul de 


perturbation du second ordre [601]. 
L3.5 Transition de transfert de charge 


Un composé transparent dans un domaine spectral peut devenir absorbant s’il est mis en 
présence d’une espèce avec laquelle il interagit par un mécanisme du type donneur-accepteur 
(D-A). Ce phénomène est lié au passage d’un électron appartenant à un orbitale liant du 
donneur (le partenaire nucléophile) vers une orbitale vacante de l’accepteur (l’électrophile), 
d’un niveau d’énergie proche. La transition est appelée transition de transfert de charge. 

La position de la bande d’absorption correspondante est en fonction du potentiel 
d’ionisation du donneur et de l’affinité électronique de l’accepteur ; la valeur de £ pour une telle 
transition est en général très grande. 

L’expression de transfert de charge possède réellement un sens absolu quand on l’utilise 
pour représenter une transition électronique qui est accompagnée par une migration de la charge 
d’une molécule à l’autre. Le terme n’a pas de sens, si l’état fondamental du complexe est par 
exemple décrit au moyen de la théorie des orbitales moléculaires ; il est cependant utilisé pour 
le passage d’une charge du donneur à l’accepteur. 

Dans les transitions intermoléculaires (Figure L8), 1l y a transfert électronique complet 
où partiel de la plus haute orbitale moléculaire occupée HOMO (HOMO = highest occupied 
molecular orbital en Anglais) à la plus basse orbitale moléculaire vide LUMO (LUMO = lowest 
unoccupied molecular orbital). 

La théorie de Mulliken s’applique seulement aux cas où le donneur et l’accepteur sont 
des molécules neutres aux couches saturées et sont tous les deux dans des états électroniques 
singlets symétriques. Les forces de transfert de charge ont des propriétés directionnelles, de 
telle sorte que les complexes sont formés dans des rapports moléculaires simples. De plus, ces 
propriétés directionnelles ont été confirmées par l’étude de la structure cristalline de nombreux 
complexes moléculaires, qui montre que les molécules des composants sont placées les unes 
au-dessus des autres, dans des plans parallèles, avec de très faibles variations de la distance 


entre les donneurs et les accepteurs. 
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Donneur Accépteur 
Orbitales Orbitales 
moléculaire a ——— } moléculaire 
anti-liante anti-liante 
LUMO 


: | HOMO | | 


Orbitales Orbitales 


moléculaire sr moléculaire 
liante liante 


Figure L.7 Mécanisme de transfert donneur-accepteur. 


L3.6 Classification des différents donneurs et accepteurs d’électron 


La notion de complexe à transfert de charge regroupe un grand nombre de types de 


complexes. Mulliken a tout d'abord prédéfini une classification des donneurs et des accepteurs. 
> Donneurs d’électrons 


Les donneurs sont ainsi divisés en fonction de leur structure électronique, c'est-à-dire 
ceux possédant un doublet libre (n), des électrons appariés (x et 6) ou des électrons célibataires 


(radicaux libres R) : 


Les donneurs (n) tels que les amines, les cétones, les alcools forment généralement des 
complexes très liés contrairement aux donneurs (6), tels que les hydrocarbures aliphatiques qui 
libèrent un électron de leur orbitale de liaison et forment alors de faibles complexes. 

Les donneurs (x), tels que les hydrocarbures et surtout les systèmes polycycliques ont tendance 
à perdre facilement un électron ce qui leur permet de former des complexes stables. Les 
donneurs (R) forment seulement des complexes avec les accepteurs (x) car avec les accepteurs 
à nombre pair des électrons (n ou x), les composés sont de très faible durée de vie ou très 


instables. Ils donnent alors souvent lieu à des réactions intermédiaires [61]. 
> Accepteurs d’électrons 


Les accepteur sont été classés, de manière analogue, par les composés possédant des 


orbitales vacantes (v), des électrons appariés (x et 6) ou des électrons (radicauxlibres Q). 
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Les accepteurs possédant des orbitales libres (v) regroupent les cations monoatomiques du type 
(Pd?*, Ag*, Cu?*), mais sont exclus les composés du type B(X)3, AI(X)3, qui forment des liaisons 
tellement fortes qu'il en résulte une liaison covalente. Les accepteurs (6) regroupent les 
halogènes (X=F, CI, Br ou I) ou les molécules halogénées du type D, Br, ICI, ICN. Les 
accepteurs (n) incluent tous les composés aromatiques possédant des substituants très 
attracteurs d'électrons tels que les cyano-, halo- ou nitro-aromatiques, ainsi que les quinones. 
De même que pour les donneurs (R), les accepteurs (Q) donnent lieu à des réactions 
intermédiaires avec les donneurs (n, o et x). Les accepteurs classiques les plus connus sont : 
2,5-Dichloro-3-dicyano p-benzoquinone (DDQ), Tétracyanoéthylène (TCNE), l'iode (1), 


l’acide picrique et le chloroaniline [611]. 


L3.7 Applications des complexes de transfert de charge 


L'interaction par transfert de charge joue un rôle important dans les systèmes 
biologiques (fixation de l'ADN) comme antibactérien, antifongique, insecticide, et dans le 
transfert des ions à travers les membranes lipophiles [62-69]. Les complexes de transfert de 
charge sont de nos jours d'une grande importance dans le domaine des semi-conducteurs 
organiques [70] et des photocatalyseurs [71]. 

En chimie analytique les complexes de transfert de charge sont intensément utilisés 
dans le dosage des substances actives dans les formulations pharmaceutiques en faisant appel 


aux diverses techniques analytiques [72]: spectrophotométrie, potentiomètrie, HPLC etc … 
L.4 Etude physico-chimiques des complexes moléculaires 


L.4.1 Spectres UV/Visible 


La spectrophotométrie UV/Visible nous permet de déterminer par mesure d'absorbance, 
l'énergie de la bande de transfert c'est-à-dire le saut d'absorption de l'état fondamental vers l'état 
excité. Cette énergie correspond à la longueur d'onde maximale de la bande de transfert obtenue. 
Ces bandes d'absorption apparaissent dans le domaine des longueurs d'onde du visible et on 
observe expérimentalement la formation du complexe par une variation de coloration des 


solutions ou des composés solides [73-74]. 


1.4.2 Constante de stabilité 


La plupart des constantes d'équilibre ont été déterminées par spectrophotométrie 


UV/Visible sur des complexes en phase liquide. En effet, nous observons non seulement 
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l'absorption des deux composés individuels mais aussi l'apparition d'une nouvelle bande 
caractéristique du transfert de charge du complexe formé. La mesure de la constante d'équilibre 
s'effectue alors dans le domaine de longueur d'onde où seul le complexe absorbe. 
Pour un complexe de stæchiométrie de 1:1, la relation d'équilibre prend la forme: 

D +A —= DA 


La constante d'équilibre thermodynamique K est alors définie par: 


K = 


[AD] 
Avec [AD] concentration à l'équilibre du complexe, [A ]concentration à l'équilibre de 


l'accepteur et [D] concentration à l'équilibre du donneur. 


Ces constantes sont fréquemment mesurées par différentes méthodes : Méthode de Benesi- 


Hildebrand, Méthode de Scott et Méthode de Foster. 


L.4.3 Les grandeurs thermodynamiques 


Les grandeurs thermodynamiques ; l’enthalpie (AH), l’entropie (AS,) et l’enthalpie 


libre (AG, ) liée à la formation des complexes et ont été évaluées à partir de relation : 
nKr=-——<+— (1-19) 

Où : 
R : Constante des gaz parfait (R=8.413 J/mol.K) ; 
T : Température en Kelvin ; 
Kr: Constante de stabilité des complexes ; 
Cette expression présente la variation des constantes de stabilité en fonction de l’inverse des 
températures « In Kr =f e) »; les valeurs de AH, et AS, sont calculées à partir de graphe 
tracé. 
Selon la relation de Van tHoff [75] on peut calculer la valeur d’enthalpie libre : 

AGRr= AH,-T AS,=-RTInK (1-20) 
R : Constante des gaz parfaits. 


T : Température en Kelvin. 


KT: Constante de stabilité des complexes de TC. 
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LS Techniques de détection des complexes moléculaires 


Les caractérisations des complexes moléculaires font appel à des techniques communes 


et parfois appropriées dépendant du mode de complexation. 


La formation d’un complexe peut être vérifiée à l’aide de différentes méthodes d’analyse. On 
trouve dans la littérature des caractérisations par fluorimétrie, diffraction de rayons X, 
spectroscopie UV-visible, spectroscopie Infrarouge (IR), Résonance Magnétique Nucléaire 
(RMN), Analyse Calorimétrique Différentielle (ACD ou DSC), Thermo Gravimétrie (TG), 
etc... 

Ces techniques ont des approches qualitatives (il s’agit juste de mettre en évidence qu’il 
y a eu complexation) et/ou quantitatives si elles permettent de déterminer la constante 


d’association K. 


L.5.1 Spectroscopie UV-VIS 


La spectroscopie d’absorption UV-visible est une technique fréquemment utilisée pour 
la détermination de la constante d’équilibre d’association dans les applications biochimiques. 
Elle est particulièrement intéressante pour la mise en évidence de la formation des complexes 
moléculaires pour les composés qui absorbent dans le domaine UV-visible. Souvent, la 
formation d’un complexe produit des modifications du spectre d’absorption de la molécule 
libre. Ces modifications se manifestent sous la forme d’un déplacement et/ou d’un 
élargissement de bande. Ce déplacement peut être attribué au changement d’environnement 


(effet s’apparentant à un effet de solvant) [76]. 
L.5.2 Résonance Magnétique Nucléaire 


La résonance magnétique nucléaire en milieu liquide est une technique utilisée pour 
l’analyse des structures de nombreuses molécules biochimiques. La RMN est aujourd’hui 
considérée comme une méthode de choix dans l’étude des complexes formés par les 
cyclodextrines [77]. 

En effet, la RMN, du proton notamment, permet d’observer des changements de 
déplacements chimiques dus à la complexation. Après inclusion, les protons situés à l’intérieur 
de la cavité de la B-cyclodextrine (H-3, H-S) voient leur environnement électronique modifié 
du fait de la présence d’une molécule invitée dans la cavité, tandis que les protons localisés à 
l’extérieur (H-2, H-4 et H-6) de la cavité restent pratiquement inchangés à cause de leur très 


faible interaction avec la molécule invitée. 
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La molécule invitée présente quant à elle, des modifications de déplacements chimiques au 
niveau de la plupart de ses protons. On peut considérer, en première approximation, que la 
molécule cage constitue le solvant de la molécule invitée. 

En ce qui concerne le complexe de transfert de charge, la RMN du proton permet de 


déterminer la géométrie exacte du complexe de TC. 


L5.3 Spectroscopie Infrarouge (IR) 


L’infrarouge est une technique qui peut être également utilisée pour caractériser les 
complexes moléculaires. D’après la littérature, pour le complexe d’inclusion, la bande OH 
caractéristique des cyclodextrines est influencée par la formation du complexe et les bandes de 
la partie incluse de la molécule invitée sont généralement déplacées ou leur intensité diminuée 
suite à la modification de l’environnement de l’invité. Pour le complexe de transfert de charge, 
les bandes d’absorption infrarouge de l’un ou des deux composés peuvent être modifiées, si 
dans le complexe de TC, les liaisons sont fortes avec des orientations relatives privilégiées, 
stabilisées par une forme ionique (D*-A°). Dans des complexes pour lesquels les liaisons sont 
faibles, les bandes d’absorption de vibration du spectre infrarouge sont très légèrement 


perturbées. 


L5.4 Spectroscopie Raman 


En spectrométrie Raman, l’analyse se fait par excitation du matériau. Porté à un niveau 
énergétique virtuel par une puissante source lumineuse monochromatique de type laser, il 
réémet ensuite une radiation qui est collectée puis analysée par un détecteur adéquat. 

Cette technique fournit des informations de natures différentes. Elle Permet d’identifier 
des composés chimiques ainsi que le type de liaison d’un composé et sa structure cristalline. 

Le but de cette approche est de chercher pour une signature particulière du complexe 
moléculaire et de le comparer à ses composants séparés. Les outils qui peuvent prévoir 
exactement les spectres et les interpréter en termes de structure moléculaire ou en termes de 
rapports structure propriété, favoriseraient considérablement l'utilité de cette spectroscopie pour 


la caractérisation. 


L.5.5 Analyse calorimétrique différentielle (DSC) 


L’analyse calorimétrique différentielle est une méthode thermo-analytique qui permet 
d’étudier les transformations physiques ou chimiques que subissent de nombreux matériaux 


pendant leur chauffage ou leur refroidissement. 
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La formation des complexes moléculaires peut être vérifiée par cette technique, par 
l’étude de l’effet thermodynamique résultant de la présence de la molécule libre (invitée ou 
accepteur). Néanmoins, cette technique est limitée aux composés ayant une température de 
fusion et d’ébullition en dessous de celle de dégradation de la cyclodextrine et du donneur. Il a 
été observé que la formation des complexes cyclodextrine-molécules invitées provoque la 
diminution du pic endothermique de fusion de la molécule invitée (cas d’un mélange physique), 


et sa disparition dans le cas d’une complexation complète [78]. 


L5.6 Rayons X 


L’étude des complexes d’inclusion par diffraction des rayons X est très souvent 
employée et permet, de plus, d’obtenir des informations sur la structure du complexe formé. 
Cette caractérisation se base sur les changements des structures cristallines de la cyclodextrine 
ou de l’invité lors du processus d’inclusion et est très utilisée lors de l’étude des complexes 
d’inclusion formés à partir de polymères [79,80]. 

Cette technique est utilisée aussi dans l’étude des structures cristallines pour les 


complexes de transfert de charge. 


1.5.7 Analyse thermogravimétrique (ATG) 


L’analyse thermique (analyse enthalpique différentielle, analyse thermogravimétrique) 
est aussi utilisée pour la détection et la caractérisation des complexes d’inclusion basés sur la 
cyclodextrine [81,82]. Les principaux effets observés sont souvent un décalage vers les hautes 
températures des différentes transitions (fusion, transition vitreuse) et une dégradation de la 
molécule invitée complexée à une température plus élevée que la température de dégradation 


de la molécule non complexée. 


L.5.8 Spectrométrie de masse 


Les spectromètres de masse sont très précis, donc très utiles pour l’analyse de mélanges 
de complexes. Ces dernières années la technique de l’Electron Spray Mass Spectroscopy (ES- 
MS) basée sur l’ionisation douce, semble être le moyen de choix pour la détection des 


complexes de transfert de charge [83]. 
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IT. Interactions non covalentes et Transfert de proton 


Le phénomène de transfert de proton et les interactions non covalentes faibles sont des 
sujets d'intérêt croissant ayant fait l'objet de nombreuses recherches scientifiques. Ils sont 
fondamentaux dans les divers domaines de la chimie et de la biologie. Il serait donc important 


de présenter quelques généralités sur ces deux concepts. 


IL.1 Les accepteurs 


Il y a eu différents types d’accepteurs décrits dans littérature pour les études de 
complexes de transfert de charge [1,2]. Récemment, beaucoup d'attention a été accordée à 
l'étude des complexes TC formés dans les réactions de composés organiques avec 6- ou x- 
accepteurs [3]. C'est en raison de leurs propriétés chimiques et physiques importantes. 

Dans des études antérieures [4,5] , nous avons rapporté la formation de complexes de 
transfert de charge stables formés à partir de la réaction de donneurs d’électrons, avec de l’iode 
et de r-accepteurs [(chloranil (CHL), acide picrique (PA), dihydroxy-p-benzoquinone, 2,4,4,6- 
tetrabromo-2,5-cyclohexadienone (TBCHD), acide 3,5dinitrosalicylique,  7,7,8,8,tetra- 
cyanoquinodimethane (TCNE), 2,3,5,6-tétrachloro-1,4-benzoquinone (LCH), 2,4,7-trinitro-9- 
fluorénone (TNF) et le tétracyanoéthyléne (TCNE), 2.3-dichloro-5.6-dicyano-1.4-benzoquin- 
one (DDQ)|]. 

Dans notre étude nous allons donner quelques caractéristiques de deux accepteurs 


d’électron, 2.3-dichloro-5.6-dicyano-1.4-benzoquinone (DDQ); tétracyanoéthyléne (TCNE). 


IL.1.1 Structure et propriétés des réactifs utilisés 
IL.1.1.1 Le 1.4 benzoquinone 


Comme les aldéhydes et les cétones, les quinones échangent l’atome d’oxygène avec 
celui de l’eau par un mécanisme d’addition-élimination avec une vitesse plus lente que les 
cétones saturées [5] (Schéma IL.1). 

D'un autre côté, la détermination de l’affinité électronique de quelques x-accepteur a 
fait l’objet d’études qui ont permis le calcul d’énergie des plus basses orbitales moléculaires 
inoccupées des accepteurs. 

Cette molécule est utilisée en synthèse organique comme accepteur de proton et comme 


agent oxydant. 
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Schéma IL.1 Mécanisme d’addition-élimination. 


Elle est, également, utilisée comme intermédiaire pour la synthèse des substances suivantes : 
Ÿ”_1,4-naphtoquinone, obtenue par oxydation du naphtalène en présence de trioxyde de 
chrome. 
Ÿ”_2.5-dichloro3.6-dicyano-1.4-benzoquinone (ou DDQ), un agent oxydant et 
déshydrogénant. 
Ÿ” Chloro-p-benzoquinone. 


Ÿ”_ Chloranile, 1,4-C6Cu4O, un puissant agent oxydant et déshydrogénant. 


IL.1.1.2 Le 2.5 dichloro 3.6 dicyano p-benzoquinone (DDQ) 


Les DDQ sont généralement connus comme accepteurs d'électrons (Fig.Il.1) et utilisés 
en synthèse chimique. Ils jouent un rôle essentiel dans les réactions biologiques [6,7]. Certaines 
de ces réactions sont simplement considérées comme modèles de processus biochimiques [8]. 
En raison de l'importance de la DDQ en tant que t-accepteur, son complexes de transfert de 


charge ont fait l'objet de nombreuses études avec des donneurs différents [9-14]. 


Figure IL.1 Structure de 2.3_dichloro_ 5,6_ dicyano p-benzoquinone (DDQ). 


En utilisant le DDQ comme agent de déshydrogénation des stéroïdes [15], ce réactif a été 


également utilisé comme agent d’oxydation des ethers [16], en passant par un complexe de 
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transfert de charge (CTC) comme étape intermédiaire. 


IL.1.1.2.1 Synthèse du 2.5 dichloro 3.6 dicyano p-benzoquinone (DDQ) 


DDR est disponible dans le commerce et sa première synthèse a été décrite par Thiele 
et Gunther en 1906 [17]. Il peut être préparé par un procédé en six étapes impliquant la 


cyanation et chloration de benzoquinone (Schéma IL.2). 


+ 2KCN + H}SO4 —> 


Schéma IL.2 Synthèse de 2,3-dichloro-5,6-dicyano-1,4-benzoquinone (DDQ). 
IL.1.1.2.2 Propriétés physiques 


a) Le point de fusion, solubilité et purification 


Le DD est un solide jaune brillant, Ts= 213-215, qui peut être gardé indéfiniment dans 
une atmosphère sèche. Il se décompose en présence d’eau malgré que la vitesse de 
décomposition soit considérablement retardée en présence d’acide fort. Le DDQ peut être 
purifié par simple recristallisation dans le dichlorométhane ou dans le mélange benzène- 
chloroforme (2/3). Il peut être purifié également par sublimation. 

Cependant, il n’y a pas des données valables dans la littérature concernant la solubilité 
de DDQ : 

Ÿ”_ Très soluble dans THF et l’acétate d’éthyle. 


Ÿ” Modérément solubles dans le dichlorométhane, toluène, dioxane, et acide acétique. 
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Ÿ”_ Insoluble dans l’eau. 
Le spectre UV de DDQ est difficile à déterminer puisque la quinone subit une auto 


décomposition dans la solution diluée, les solutions concentrées subissent une décomposition 
négligeable pendant plusieurs heures, le DDQ montre un maximum à 390 nm dans le dioxane 
et un autre à 410 nm dans le benzéne. 

Le spectre IR du DDQ montre une absorption faible à environs 2210 cm (CN) des 
absorptions fortes à environs 1680 (quinone) et 1555 cm '(C=C). 

Le DDQ ne montre aucune absorption (à plus de 3200 nm) quand il est pur, mais quand 
il est mélangé avec l’hydroquinone, un large signal radical attribuable à la formation semi 


quinone, apparaît immédiatement à 3300 nm. 


b) Oxydation, potentiel de réduction 


Comme il doit être attendu, les résultats dans les systèmes hydroxyliques présentent des 
anomalies, cela est particulièrement vrai pour les halogenoquinone. Dans ce cas, les variations 
de la solvatation et plus particulièrement les effets produits par liaison hydrogène, contribuent 
d’une manière imprévisible, cependant les potentiels de réduction sont liés à l’affinité 
électronique de la quinone. 

Les quinones sont de bons candidats pour les réactions éléctrocycliques car elles 
possèdent des doubles liaisons carbone-carbone et carbone-oxygène réactives, les p-quinones 


sont également des diénophiles pour les réactions de Diels-Alder [18]. 


IL.1.1.2.3 Propriétés chimiques 
a) Déshydrogénation 

L'introduction d’une double liaison par l’abstraction d’un ion hydrure et l’élimination 
du proton constitue la réaction primaire du DDQ. Quelques exemples sont rapportés dans la 
littérature montrant la formation de radicaux produits par la quinone. 

La réactivité du DDQ a fait l’objet de nombreuses études. Entre autres, la 
déshydrogénation des stéroïdes possédant une activité biologique désirable. 

Le DDAQ par sa grande affinité électronique, forme également des complexes de TC dans 


la majorité de ses réactions. 


IL.1.1.3 Le tétracyanoéthyléne (TCNE) 


Le tétracyanoéthyléne a été découvert dans les laboratoires de recherche de Du Pont de 


Nemours [19] en 1957 (Figure IL.2). 
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Figure IL2 Structure de tétracyanoéthyléne (TCNE). 


Ces industriels cherchaient, alors, à faire un polymère percyané analogue au téflon mais 
dont il était attendu que les propriétés soient encore meilleures, puisque ses liaisons sont encore 
plus polarisées [20]. Mais malheureusement, à ce jour, 1l n’a jamais été possible de le faire 
polymériser, ou même copolymériser. Toutefois, par les propriétés exceptionnelles mises à jour, 
ces chercheurs de Du Pont ont ouvert la voie à la chimie très riche des composés cyanocarbonés 
[21]. 

La synthèse du Tétracyanoéthylène procède par dibromation du malonitrile puis 


couplage du dibromomalonitrile en présence de cuivre (Schéma IL3) 


8Br>, KBr 
ice KBr[BrC(CN)|4 
H;,0. A 


Cu powder 


KBr [Br C(CN);| 4 


benzene, A 


Schéma IL3 Synthése de tétracyanoéthylène (TCNE). 


La réactivité exceptionnelle du tétracyanoéthylène vient de l’association de plusieurs 
propriétés. 

Il y a d’abord sa capacité à réagir comme un bon diénophile. Il réagit aussi facilement 
dans des réactions thermiques et non-concertées de type [2+2] avec les oléfines riches en 
électrons (schéma II.6). 

Autre réactivité apportée par l’importante électro déficience que lui confèrent ses quatre 
substituants cyano, il est très prompt à réagir avec les nucléophiles comme les amines, les 
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alcools, l’eau, les cétones énolisables et les anilines. Leurs attaques peuvent être suivies par 
l’élimination d’un groupement cyano, qui manifeste ainsi le rôle de bon groupe partant de ces 
derniers. C’est aussi un très bon accepteur pour les réactions de Michael. 

Enfin il a de remarquables propriétés d’oxydant du fait de la facilité avec laquelle il 
forme des radicaux et des diradicaux d’anions en acceptant réversiblement un ou deux électrons 


[221 (schéma IL). 


Schéma IL.4 Réactivité de tétracyanoéthylène : a) cycloaddition [2+2], b) avec les 


nucléophiles, c) comme oxydant. 


Le TCNE est un agent remarquable, employé souvent comme accepteur d’électron, et 
joue un rôle important dans l’étude des complexes de TC avec les composés aromatiques [23] 
riches en électrons, qui recristallise dans le chlorobenzène et s’évapore à 198- 200°C. 

Il est aussi capable de jouer un rôle de ligand et est capable de lier jusqu’à quatre métaux 
avec des applications dans les matériaux présentant des propriétés magnétiques. Il est employé 


pour préparer des nombreux supraconducteurs organiques [24]. 


IL.2 Interactions non covalentes 


En chimie et en biochimie, les interactions faibles sont un sujet d'intérêt majeur, 
notamment dans les domaines émergents de la reconnaissance moléculaire, de l'auto- 
assemblage, la chimie supramoléculaire et les interactions générales hôte-invité. Les 


interactions faibles sont fondamentales dans de nombreux aspects de la biologie, l'exemple par 
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excellence étant la double hélice élégante des interactions intermoléculaires faibles de l'ADN, 
résultant d'un jeu subtil des interactions faibles. 

Les interactions non covalentes de divers types, telles que la liaison hydrogène, dipôle 
—dipôle interactions et dispersion de London , auxquelles d’ajoutent les empilements r-x, les 
interactions CH/x , les liaisons halogène et les interactions répulsives, jouent un rôle clé dans 
l'influence des réactivités et structures moléculaires, dans la reconnaissance moléculaire et la 
haison inter et intramoléculaires, et dans la structure, la stabilité et les propriétés dynamiques 
des biomolécules [25-28]. 

Les forces non covalentes sont essentielles pour déterminer la structure secondaire, 
tertiaire ou quaternaire des grandes molécules, et donc leur fonction macromoléculaire. 
Comprendre l'ampleur relative de ces différents types d'interactions est important, en particulier 
telles qu'elles peuvent exister dans la coexistence ou la concurrence et ce qui a fait l'objet des 


efforts expérimental et théorique intenses [29-38]. 


Tableau IL.1 Liaisons faibles et énergies de liaisons correspondantes 


Interaction Energie de liaison (kJ/mole) 

Interaction de van der Waals 1-5 

Interaction 7-7 7-20 

Liaison hydrogène 10-20 
interaction hydrophobe 12-15 

Liaison ionique 12-20 

Liaison covalente de coordination 40-120 

Liaison covalente 150-1000 


La stabilité de ces systèmes moléculaires repose sur l'établissement des interactions non 
covalentes de faible énergie qui dépendent de la géométrie et de la nature des groupements 


chimiques assemblés [39,401]. 


Ces interactions intermoléculaires forment la base des processus hautement spécifiques et 
jouent un rôle fondamental dans des domaines variés, avec par exemple, la fixation d'un substrat 
à un récepteur protéique, les réactions enzymatiques, l'association immunologique antigène- 


anticorps et tous les autres aspects de la communication cellulaire. 
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IL.2.1 Liaisons hydrogène 


La Liaison hydrogène est une interaction attractive entre un atome d’hydrogène liéde 
façon covalente à un atome donneur fortement électronégatif tels que l’azote (électronégativité 
:4=3), l'oxygène (4 = 3,5) et le fluor (4 = 4) [41], et un atome accepteur similaire faisant partie 
ou non de la même molécule que celle du donneur. La longueur de la liaison hydrogène est 


définie comme la distance entre les centres de ces deux atomes (Figure IL.3). 


Liaison hydrogene 


Longueur de la liaison hydrogène 


Figure IL.3 Représentation schématique de la liaison hydrogène [42]. 


On définit la liaison hydrogène par trois variables (Figure IL.4) [43] : 
- La distance donneur-hydrogène : D-H. 
- L'interaction hydrogène-accepteur : H...A. 


- L'angle donneur-hydrogène-accepteur : D-H...A. 


Figure IL.4 Les trois paramètres géométriques. 


L'interaction entre le groupe D-H et l’atome A est de type électrostatique. Le nuage 
électronique de l'hydrogène est attiré par l’atome donneur qui est relativement plus 
électronégatif que l’atome d’hydrogène [44,451 créant ainsi une charge partielle positive sur 
l'hydrogène. Cette charge positive est attirée par la charge partielle négative portée par l’atome 


accepteur donnant ainsi naissance à une interaction désignée par pont hydrogène (Figure ILS). 
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Figure ILS Les charges partielles dans la liaison hydrogène [46]. 


L'interaction est toutefois suffisante pour que la distance entre l'atome d'hydrogène et 
l’atome accepteur soit inférieure à l’interaction de van der Waals. La plupart des molécules 
polaires, en particulier les molécules biologiques, peuvent donner des liaisons hydrogène et se 
lier entre elles ou avec d’autre molécules y compris les molécules d’eau; qui est l’exemple idéal 


pour comprendre cette liaison. 


Cas 1 : 2 molécules différentes Cas 2 : 2 molécules identiques 


liaison hydrogène liaison hydrogène 


Ÿ 


\ & ë F/ 


N— Hunt) ———p> D Hinmunnmmu() 
{ L' 


\ p2 
\ J 
R y 14 


\ / Sn 4 


Donneur de Accepteur de Donneur de Accepteur de 
liaison hydrogène liaison hydrogène liaison hydrogène liaison hydrogène 


Figure IL.6 Liaison hydrogène entre deux molécules : différentes et identiques. 


On peut mentionner les remarques suivantes : 

a) L’atome accepteur possède une paire d’électron libre ou des électrons x polarisables qui 
contribuent à stabiliser la liaison hydrogène. 

b) l'énergie de la liaison hydrogène est maximale lorsque les trois atomes impliqués (D)- 


(H)... (A) sont colinéaires. Les liaisons hydrogènes habituelles sont : 
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Figure IL7 Liaisons hydrogène habituelles. 


Dans certains cas particuliers, le groupe C-H peut se comporter comme donneur de 
haisons hydrogène. Cependant, les liaisons hydrogène impliquant ces atomes sont souvent très 
faibles. 

La liaison hydrogène à trois propriétés : 

> Comme les liaisons covalente, la liaison hydrogène est directionnelle : elle s’aligne 
suivant l’axe de la liaison covalente qui lui est associée. Cette propriété permet d'avoir 
des architectures moléculaires très bien définies comme dans la glace ou dans l'eau, ou 
encore dans les molécules biologiques telles que les protéines ou l'ADN où les deux 
brins sont liés entre eux par des liaisons hydrogène. 

> Les liaisons hydrogène peuvent se tordre, se rompre ou se restaurer en fonction de la 
température. Cette propriété donne aux architectures moléculaires assemblées par 
liaisons mêmes températures. Cette souplesse et ces possibilités d'évolution sont 
indispensables aux molécules biologiques. 

> Enfin, la liaison hydrogène est capable de transférer des ions H* entre les molécules 
qu'elle lie. Cette propriété est très importante car elle est à l’origine de la réactivité des 
milieux aqueux. Sans elle, ceux-ci seraient inertes et la vie ne serait pas possible, car les 
molécules biologiques doivent en permanence réagir pour rester actives. On soupçonne 
que ce sont les molécules d’eau qui permettent ces transferts d’ions H* en établissant 
certaines liaisons hydrogène spécifiques, en d'autres termes qu'elles donnent le feu vert 
pour que des molécules telles que les protéines puissent réagir. C'est peut-être là leur 


rôle fondamental, et ce serait grâce à cela que la vie est née dans l'eau et s'y poursuit ! 


IL.2.1.1 Les différents types de liaisons hydrogène 


La distance entre les deux atomes électronégatifs est en rapport inverse avec le type de 
liaison, cette dernière est faible et tend vers l’interaction de van der Waals quand la distance 


est grande. La plus courte liaison hydrogène tend vers la liaison covalente. 
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Les liaisons hydrogène peuvent être classées dans trois catégories différentes : liaisons fortes, 


liaisons modérées et liaisons faibles (Tableau IL.2). 


Tableau II.2 Propriétés des liaisons hydrogène [45]. 


Liaisons Fortes Modérées Faibles 
Type d'interaction Majoritairement | Majoritairement Electrostatique 
D-H...A covalente électrostatique 
Distances D-H = H...A D-H < H...A D-H << H...A 
Langueur H...A (À) 1.2-1.5 1.5-2.2 2.2-3,2 
Langueur D... A (À) 2.2-2.5 2.5-3,2 3.2-4.0 
Angle D-H...A (°) 175-180 130-180 90-150 
Energie de la liaison 14-40 4-15 <4 
(kcal. mol'!). 
-Dimères en phase | - Acides. -Dimères en phase 
gazeuse d'acides | - Alcools. gazeuse d'acides faibles 
Exemples forts ou de bases | - Phénols ou de bases faibles. 
fortes. - Hydrates. -Liaisons type C-H..O/N. 
-Pseudohydrates. -Molécules -Liaisons type O/N- 
-Complexes biologiques ELST 
d'acide 
fluorhydrique. 


IL.2.1.1.1 Liaisons hydrogène fortes 


Les liaisons hydrogène fortes font partie de la catégorie de liaisons hydrogène la plus 
intrigante. Elles sont formées quand le proton est partagé entre deux fortes bases, ou entre ions 
et molécules dans le cas où 11 y a une déficience dans la densité électronique pour le groupement 
donneur ou bien un excès de densité électronique pour le groupement accepteur [41,45]. Elles 


sont connues dans ce cas sous le nom de liaisons hydrogène fortes. 


IL.2.1.1.2 Liaisons hydrogène modérées 


Les liaisons hydrogène modérées, conventionnelles, ou normales [47] sont formées 
entre un donneur et un accepteur généralement neutres. L'atome donneur est relativement plus 
électronégatif que l'atome d'hydrogène est l'atome accepteur possède un doublet libre 
d'électrons. Ce sont les liaisons les plus communes en chimie et dans la nature [45]. On les 


considère comme des liaisons hydrogène normales. 


IL.2.1.1.3 Liaisons hydrogène faibles 


Une liaison hydrogène faible est généralement formée quand l’atome d’hydrogène 


Chapitre II Interactions non covalentes et Transfert de proton 


est lié par covalence à un atome légèrement plus électronégatif, tel que l’interaction C- H...A. 
Elle peut aussi se former lorsque l’accepteur n’a pas de doublet d’électron libre, mais des 
électrons x. L’énergie de ce type de liaison est comparable aux interactions de Van der Waals 
[41,45]. Les meilleurs exemples pour ce type de liaisons sont celles observés dans les phases 
gazeuses comme dans les réactions d’addition de HF, HCI, HBr et HCN avec N2, CO, OCS et 
CO qui ont une énergie d’environ 3 kcal.mole’! [48]. 

Dans certains cas l’atome accepteur participe avec des électrons x comme dans un cycle 
aromatique. 
Les exemples de groupes donneurs et de groupes accepteurs engageant une liaison hydrogène 


forte, modérée ou faible sont donnés dans le tableau (tableau IL.3). 


Tableau IL.3 Groupes donneurs et accepteurs des différentes catégories de liaisons hydrogène 


[45] 
Liaisons hydrogènes fortes 
Donneurs et accepteurs 
[F..H..F] Ion symétrique d'hydrobifiluorure 
[F..H..F]; (ion hydrobifluorure). 


Anions dans les composés contenant HF. 
Anions organiques contenant de l'hydrogène. 


[O — H..0”] Hydrogène phisphates et sulfates, ions carboxylates. 
[0O*H..0] 

[N*H.. N] Ions hydronium, pseudo hydrates, éponges de protons. 
[N..H..N] 


Liaisons hydrogènes modérées 
Donneurs et accepteurs 


O..H, P-O-H, H-Ow-H Eau, hydrates, alcools, acides carboxyliques, phénols, 
carbohydrates, oligo et polysaccharides nucléosides, 
nucléotides et acides nucléiques. 


(0 N Amines secondaires, amides, carbamates, hydrazides, 
NH NH purines, pyrimidines,  barbiturates,  nucléosides, 
C C nucléotides, peptides, protéines (chaîne principale et 


chaîne latérale). 
Donneurs seulement 


NH Sels d'ammonium. 


-NH; 
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NH Acides aminés sous forme de zwitterions. 
S-H Cystéine. 
a Protéines (chaîne latérale, acides nucléiques à bas pH) 
c_ 
C-N(DH Amines primaires, pyrimidines, purines, barbiturates. 


Accepteurs seulement 


C Ethers, carbohydrates, oligo- et polysaccharides. 
à. 
ce” 
Carboxylates, Acides aminés sous forme de zwitterions, 
Xe acides carboxyliques, cétones, esters, N-oxides, 
nu pyrimidines, purines, nucléosides, nucléotides, acides 
nucléiques, peptides, protéines (chaîne principale). 
X=0O Oxyanions, nitrates, chlorates, sulfates, phosphates. 
NN Amines tertiaires. 
T 
be Purines,  pyrimidines,  barbiturates,  nucléosides, 
DA nucléotides, acides nucléiques. 
N=O Composés nitriques aromatiques. 
e éthionine. 


Liaisons hydrogène faibles 
Donneurs 
C-H 


Accepteurs 


C=C (© 


IL.2.1.2 Caractéristiques des liaisons hydrogène 


Dans les structures cristallines, l’empilement des molécules est déterminé par leur forme 
ainsi qu’une variété de forces intermoléculaires dont les premiers sont les liaisons hydrogène. 


On peut classer les liaisons hydrogène selon leurs attachements en six catégories : 


IL.2.1.2.1 Liaison d’hydrogène à quatre centres 


Les liaisons à quatre centres (à trois groupements accepteurs), sont rarement observées 
dans les structures cristallines (<5%). Les distances H...A sont généralement plus longues que 
les liaisons d’hydrogène à trois centres et tous les angles D-H...A doiventêtre supérieurs à 90° 


[43,47] (Figure IL8). 
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Figure ILS8 Liaison hydrogène à quatre centres. 


IL.2.1.2.2 Liaison hydrogène à trois centres 


L’atome d’hydrogène est lié à trois atomes ; avec une liaison covalente et deux liaisons 
hydrogènes, Ces trois atomes (D, Al, A2) se trouvent dans un plan ou la somme des angles est 
proche de 360° [43,47]. Les liaisons à trois centres sont usuellement non symétriques où rl # 


r2 (Figure IL.9). 
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Figure IL.9 Liaison hydrogène à trois centres. 


IL.2.1.2.3 Liaison hydrogène à deux centres 


La liaison hydrogène forte est presque linéaire avec un seul atome accepteur, elle est 


dite à deuxcentres [49-52] (Figure IL. 10). 


D —1--—- A 


Figure IL.10 Liaison hydrogène à deux centres. 


IL.2.1.2.4 Liaison hydrogène chélate 


L’atome d’hydrogène est lié à deux atomes ; une liaison covalente et une liaison 
hydrogène ou à trois atomes ; une liaison covalente et deux liaisons hydrogène avec deux 


atomes accepteurs qui sont liés au même atome X [53] (Figure IL.11). 
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Figure IL.11 Liaison hydrogène chélate. 


IL.2.1.2.5 Liaison hydrogène bifurquée 


La combinaison de deux configurations à trois centres conduit à une configuration 
bifurquée, où ces liaisons à trois centres ne sont pas coplanaires et font un angle de l’ordre de 


50° [54] (Figure IL.12). 


Figure IL.12 Liaison hydrogène bifurquée. 


IL.2.1.2.6 Liaison hydrogène tandem 


La liaison hydrogène tandem est surtout observée dans la structure de la molécule d’eau 


oules atomes hydrogène sont désordonnés [54] (Figure IL.13). 


Figure IL.13 Liaison hydrogène tandem. 


IL.2.2 Forces van der Waals 


Le traitement théorique des forces interagissant entre les atomes et les molécules, 
maintenant connu sous le nom de van der Waals ou forces de dispersion, remonte aux années 
1920. Suite à la prise de conscience que des forces attractives doivent exister entre les atomes 
ou les molécules, des tentatives ont été faites pour fournir une compréhension de leur origine. 


Les premiers travailleurs ont dérivé les forces d'interaction entre les atomes individuels dans 
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une vapeur, où les principales interactions se font par paires [55-57]. Mais il s'est avéré difficile 
de généraliser ces résultats aux milieux condensés car il existe des interactions simultanées 
entre tous les atomes ou molécules. Plus tard, une approche différente a été développée [58,591]. 
Des expressions ont été trouvées pour l'interaction entre les corps macroscopiques sur des 
distances qui étaient longues par rapport aux espacements interatomiques, et 1l a pu être 
démontré que les résultats antérieurs pour les particules dans la vapeur étaient un cas particulier 
de la nouvelle théorie. 

Du côté expérimental, ces méthodes ont donné une prescription claire pour recueillir les 
informations nécessaires au calcul précis des forces de van der Waals. En 1948, la soi-disant 
théorie DLVO à été établie pour expliquer la stabilité des dispersions colloïdales [60], et c'est 
toujours la base des théories modernes, affirmant que l'attraction électrodynamique s'oppose à 
la répulsion électrostatique. 

Les forces électrodynamiques sont de plus longue portée et dominent donc à une séparation 


suffisamment grande (Figure IL. 14). 
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Figure IL.14 Potentiel d'interaction ;; de deux atomes en fonction de la distance r;; 


La portée des forces d’attraction de van der Waals varie entre des distances atomiques jusqu'à 
des distances colloïdales ce qui implique que l’influence attractive d’une particule peut 
s’étendre sur des distances considérables à partir de sa surface. Ces forces dites de longue 
portée, sont à l’origine du phénomène de floculation qui se traduit par l’agglomération des 
particules, cependant il existe divers mécanismes qui permettent d’écranter ces forces 
d’attraction comme par exemple en favorisant les répulsions électrostatiques ou par l’adsorption 
de polymères à la surface des particules. Par la suite, on va caractériser ces interactions de van 


der Waals sous forme d’énergie potentielle d’interaction. Tout d’abord une approche 
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microscopique sera adoptée pour ensuite considérer le cas des interactions entre corps 


macroscopiques. 


IL.2.2.1 Approche microscopique 


Pour comprendre l’origine des forces d’attraction entre particules colloïdales, il faut 
d’abord se situer au niveau moléculaire et considérer les interactions individuelles entre 
molécules. Ici, on va considérer trois types de forces intermoléculaires qui se traduisent par 
l'attraction entre particules colloïdales. Ensemble ces trois interactions sont connues comme les 
interactions de type van der Waals et ont lieu entre les dipôles (induits ou permanents) de 
chaque molécule: 

- Interaction dipôle permanent — dipôle induit 

- Interaction dipôle permanent — dipôle permanent 

- Interaction dipôle induit — dipôle induit 

Les interactions entre dipôles (permanent ou induit) sont le résultat du champ électrique 

produit par un dipôle, agissant sur un deuxième dipôle. Un traitement approfondi et détaillé de 
ces interactions implique que l’on définisse l’expression du champ électrique produit par un 
dipôle ainsi que l’effet de ce champ électrique sur une particule voisine (sans dipôle). On se 
limitera dans ce chapitre à donner les expressions de chacune des trois interactions sous forme 
d’énergie potentielle notée ®. Si nécessaire des informations supplémentaires peuvent être 


trouvées ailleurs [61]. 
i) Interaction dipôle permanent — dipôle induit 

L’expression de l’énergie potentielle totale pour deux molécules interagissant entre elles 
(une avec un dipôle permanent et l’autre avec un dipôle induit) est : 


= (ae, ui + C2) d”® 
(4TE0)? 


D) (I.1) 


où &,est la polarisabilité électronique, 4 le moment dipolaire, d la distance de séparation entre 
les deux molécules et £, la permittivité du vide. L’indice 1 correspond à la molécule avec le 
dipôle permanent et l’indice 2 à la molécule avec le dipôle induit. Cette interaction dite 
interaction d’induction ou de polarisation, s’effectue entre une molécule polaire et une non - 
polaire qui se polarise (modification du nuage électronique de la molécule) sous l’effet du 
champ électrique de la molécule polaire. L’équation (IL.1) est connue aussi comme l’équation 


de Debye [62,63] (indice D). 
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ii) Interaction dipôle permanent — dipôle permanent 


L’énergie potentielle pour ce type d’interaction s’écrit : 


Dy=— 2H GE (11.2) 
3 kBT(4TEo)? 

L’équation (IL.2) est connue comme l’équation de Keesom [62,63]. Ce type d’interaction 
a lieu entre molécules polaires et elle tient compte de l’orientation (rotation) des molécules. En 
l’absence de toute contrainte extérieure, les dipôles ont tendance à s’aligner mais cet alignement 

peut être perturbé sous les effets de la température et de l’agitation thermique. 
Cette interaction est également connue comme interaction d’orientation, elle peut être attractive 
ou répulsive puisque la rotation des molécules implique qu’il y a un changement de signe. En 


général dans le cas des liquides et des gaz à température ambiante, cette interaction est faible. 


iii) Interaction dipôle induit — dipôle induit 


Pour ce troisième type d’interaction appelée interaction de London [62,63], l'expression 


pour l’énergie potentielle est issue de la mécanique quantique et elle s’écrit: 


3 hVAe;lez; ,_6 


Dre 4 (4TE0)2 


(IL3) 


Où A représente la constante de Planck et v la fréquence de vibration associée au mouvement 
des électrons autour du noyau atomique. L’interaction de London ou de dispersion, s’effectue 
entre deux molécules non-polaires; chaque molécule non polaire possède un moment dipolaire 
moyen nul, pourtant à chaque instant il existe un moment dipolaire non nul associé aux positions 
instantanées des électrons par rapport aux protons situés dans le noyau des atomes qui forment 
la molécule [64]. Ce dipôle instantané va générer un champ électrique qui à son tour va polariser 
un atome d’une molécule voisine et ainsi induire un moment dipolaire. L’interaction résultante 
entre les dipôles donne naissance à une force instantanée d’attraction entre les deux molécules. 

Une fois définies les trois expressions pour chacune des interactions 1l est important de 
faire quelques remarques : tout d’abord il faut dire que les expressions correspondent au cas où 
deux particules interagissent dans le vide. Deuxièmement, pour les trois équations la 
dépendance vis-à-vis de la distance de séparation est en 1/ d. Les paramètres moléculaires tels 
que le moment dipolaire ou la polarisabilité jouent le rôle de facteurs de proportionnalité. De 
même, on constante que les interactions d’induction et de dispersion sont toujours attractives 


contrairement à l’interaction d’orientation qui peut être attractive ou répulsive. Il faut également 
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signaler que les forces de van der Waals sont anisotropes (à l’exception des interactions entre 
deux atomes appartenant à des gaz parfaits) du fait que la polarisabilité des molécules change 
de valeur en fonction de l’orientation des molécules ( &,,, Æ &e,), donc la polarisabilité est 
aussi anisotrope. L’énergie potentielle d’interactions de van der Waals pour une paire de 


molécules identiques s’écrit : 


Dypw — ®, + Dx + D; (IL. 4) 
_ _ __1 1 LAN. je 2) y-6_ __ Cvow 
Dypw = TERRE (2 au + 7 + . hvaë) d °= pa (IL.S) 


Où la constante Cvpw pour un système donné, s’exprime en J.m. 

Parmi les trois types d’interactions qui contribuent à l’énergie potentielle de van der 
Waals, en général l’interaction de London (appelée aussi de dispersion) prédomine. Ceci 
s’explique en partie du fait que les deux autres interactions (Debye et Keesom) ne sont pas 
toujours, contrairement à l’interaction de London qui s’applique pour toutes les molécules. 
Cependant lorsqu'on a affaire à des solvants très polaires comme l’eau, l’interaction de London 
n’est pas forcément dominante. 

Les interactions de van der Waals jouent un rôle important dans divers phénomènes tels 
que l’adsorption ou l’agrégation de particules en phase aqueuse, de même elles ont une 


influence sur les propriétés des gaz et des liquides ainsi que sur leur solubilité. 


IL.2.2.2 Cas particuliers 


11.2.2.2.1 Interactions de recouvrement 


Lorsque deux atomes (molécules) s’approchent l’un(e) de l’autre, il y a interpénétration 
de leurs nuages électroniques. D’après le principe d’exclusion de Pauli, des spins de même sens 
ne peuvent pas occuper la même position dans l’espace. Aïnsi, ils se repoussent mutuellement. 
Ce phénomène correspond à la répulsion de Pauli. Lorsque la distance entre les deux systèmes 


devient très faible, alors la répulsion de Pauli prédomine sur l’interaction électrostatique. 


IL.2.2.2.2 Interactions hydrophobes 


Dans une solution aqueuse, une molécule apolaire tend à interagir avec d’autres 
molécules apolaires plutôt que de s’entourer de molécules d’eau. On nomme cette propriété 
hydrophobicité et cela conduit à des interactions entre molécules qu’on appelle des interactions 


hydrophobes. L’origine thermodynamique de l’interaction hydrophobe peut s’expliquer de la 
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façon suivante. Considérons l’enthalpie libre d’une solution contenant de l’eau et des molécules 


non polaire (équation IL.6) : 
AG = AH -TAS (11.6) 


Comme pour tout processus, il faut examiner les contributions enthalpiques (AH) et 
entropiques (AS) pour comprendre les éléments favorables ou défavorables à une interaction. 
Du point de vue enthalpique, l’énergie d’interaction est forte entre les molécules d’eau (réseau 
de liaisons hydrogènes) alors qu’elle est plus faible entre l’eau et les molécules apolaires, ce 
qui conduit à la faible solubilité des molécules apolaires dans l’eau. La présence d’une molécule 
apolaire dans de l’eau est donc défavorable enthalpiquement. Du point de vue entropique la 
situation est également défavorable: paradoxalement, la présence d’une molécule apolaire dans 
l’eau accroît la structuration de l’eau, conduisant à une diminution de l’entropie. 

On peut comprendre ce phénomène en observant que l’insertion d’une molécule apolaire 
dans de l’eau nécessite de créer une cavité dans le solvant en rompant des liaisons hydrogènes, 
d’insérer la molécule apolaire dans cette cavité et de réorganiser le réseau de liaison hydrogène 
pour contourner cet « obstacle ». Pour maximiser les liaisons hydrogènes dans cette dernière 
étape, les molécules d’eau doivent adopter une géométrie très précise pour envelopper la 
molécule apolaire. C’est ceci qui conduit à une structuration locale plus importante de l’eau et 
donc à une diminution de l’entropie. Dans ces conditions, les deux termes, enthalpie et entropie, 
sont défavorables à la présence de molécules apolaires dans l’eau : par conséquent, l’eau pousse 
au regroupement de ces molécules apolaires entre elles, reliées par des liaisons de type van der 
Waals. 


Il est important de noter que la plus grande partie de l’énergie d’interaction est due non 
pas à ces interactions de van der Waals (qui sont faibles) mais aux contributions défavorables 
à la solvatation de ces molécules apolaires dans l’eau. Donc en absence d’eau, comme par 
exemple dans le vide du spectromètre de masse, seule subsiste la part de l’énergie d’interaction 
due aux interactions de van der Waals. 


IL3 Transfert de proton 


Le transfert de proton à travers la liaison hydrogène est une des plus simples et 
fondamentales réactions en chimie et en biochimie [65-74]. IL intervient dans plusieurs 
réactions d’oxydoréduction [75-77]. 

Le transfert de proton en phase gazeuse vers une molécule M est défini par la réaction 
IL. 7, dans laquelle tous les réactifs et produits sont en phase gazeuse. La réaction IL.7 représente 
un transfert de proton non dissociatif, c'est-à-dire que le proton saute simplement d'une 
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molécule à une autre sans autre changement. Cependant, nous verrons plus loin que ce 
processus de transfert de proton aboutit parfois à une fragmentation de la molécule réceptrice, 
appelée transfert de proton dissociatif. 

La discussion ci-dessous traite d’abord de la thermodynamique des réactions de transfert 


de protons en phase gazeuse. 
XH* +M — MH*+X (1.7) 


IL.3.1 Thermodynamique du transfert de proton 


La thermodynamique peut être utilisée pour prédire si une réaction peut ou non se 
produire dans des conditions données. La quantité clé requise est le changement d'énergie de 
Gibbs pour la réaction, AG, où G représente l'énergie de Gibbs et À fait référence à la différence 
de valeurs entre les produits et les réactifs. La réaction (117) sera thermodynamiquement 
spontanée si le changement d'énergie de Gibbs pour le processus est négatif. Une façon utile 
d'interpréter cette quantité pour une réaction de transfert de protons est de reconnaître qu'elle 
dérive d'une différence dans les changements d'énergie de Gibbs pour les deux réactions 
suivantes : 


M(g) + H°(g) — MH*(g) (ILS) 
X(g) + H'(8) — XH*(g) (IL.9) 


Puisque l'énergie de Gibbs est une fonction de l'état, la loi de Hess s'applique et donc le 


changement d'énergie de Gibbs pour la réaction (11.7) peut être exprimé comme : 
AG? (IL.7) = AGŸ (IL8) — AGŸ (IL.9) (IL.10) 


Où T est la température et le O en exposant fait référence à l'état standard. 
Le changement d'énergie de Gibbs est lié à la constante d'équilibre, Keq, pour la réaction 2.3 


par les expressions suivantes : 


AG9 =-RT In Keq (L.11) 


[MH*] [x] 
= MIRE] (IL.12) 


Dans les équations ci-dessus, T est la température absolue, R est la constante universelle 


des gaz et les crochets autour des composés chimiques se réfèrent à leurs concentrations en 
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phase gazeuse. 

En acceptant des protons, les réactifs M et X dans les réactions (IL8) et (11.9) agissent 
comme des bases. La volonté d’une substance d’accepte un proton est commodément 
représentée par sa basicité, qui est définie comme le négatif du changement d’énergie de Gibbs 
pour le processus d’acceptation du proton. Imaginez que vous mélangez des quantités à peu 
près égales de M et de XH* dans un récipient fermé et que vous obteniez des quantités à peu 
près égales de MH et XH* lorsque l’équilibre est atteint. Cela correspond à une constante 
d’équilibre avec une valeur proche de l’unité puisque le numérateur et le dénominateur de 
l'équation IL.12 seront presque égaux, ce qui implique à son tour que AGŸ # 0. Si le composé 
M a une basicité supérieure à celle du composé X, alors AGT est négatif pour la réaction IL.7 et 
la position d’équilibre est inclinée vers les produits. Cela peut être quantifié à l’aide de 
l’équation IL.11 Par exemple, si le composé M a une basicité d’environ 10 kJ mol ! au-dessus 
de celle du composé X, la substitution dan l’équation IL.11 donne Ke Æ 60 à une température 
de 298 K. Il est clair qu’une différence relativement faible dans les basicités entraînera un biais 
important de la réaction vers des produits ou des réactifs. 

Des tableaux complets de basicités en phase gazeuse pour une large gamme de composés est 
fournis par Hunter et Lias [78] et une liste de basicités pour certaines classes illustratives de 
composés est fournie dans le tableau IT.4. Ces basicités en phase gazeuse peuvent être utilisées 
pour déterminer la spontanéité d’une réaction de transfert de proton particulière, mais il est plus 
courant d’utiliser à la place d’affinité du proton. L’affinité pour les protons d’un composé est 
définie comme le négatif du changement d’enthalpie pour la réaction d’acceptation des protons 
(par exemple, réaction IL8 ou IL9). Ainsi, si l’acceptation des protons est un processus 
exothermique, l’affinité des protons est une quantité positive. La basicité en phase gazeuse (B) 
et l’affinité protonique (AP) sont liées par : 

AP=B-TAS? (11.13) 
Qui découle de l’expression bien connue 

AG? = AH$ -T AS£ (IL.14) 
Où H est l’enthalpie et S est l’entropie. L'hypothèse derrière l’utilisation d’un changement 
d’enthalpie (comme dans l’affinité du proton) au lieu du changement d’énergie de Gibbs 
(comme dans la basicité) pour juger de la spontanéité de la réaction est que le changement 
d’entropie, ASŸ, dans l’équation 11.13 est négligeable. Dans la plupart des cas, il s’agit 
probablement d’une bonne approximation car un proton est une entité simple dont le transfert 


d’une molécule à une autre est peu susceptible de modifier l’entropie nette de système de 
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manière substantielle. Ainsi, dans presque tous les cas pratiques, le processus de transfert de 
proton représenté dans la réaction IL7 sera thermodynamiquement spontané si l’affinité 


protonique de M est supérieure à celle de X. 


Tableau IL.4 Exemples de basicités en phase gazeuse et l’affinité protonique. 


Classification Molécule Basicité Affinité protonique 
(kj mor!) (kj mol°!) 
Gaz inorganique O2 396 421 
N2 465 494 
CO2 516 541 
O3 626 626 
H20 660 691 
NE 819 854 
Alcanes Méthane 521 544 
Ethane 570 596 
Propane 608 626 
Cyclopropane 722 750 
Alcènes Ethène 652 681 
Propène 723 732 
Alcynes Acétylène 617 641 
Propyne 723 748 
Hydrocarbure aromatique | Benzène 725 750 
Toluène 756 784 
Alcools Méthanol 725 754 
Ethanol 746 776 
Phénol 786 817 
Ethers Ether diméthylique 765 792 
Aldéhydes Formaldéhyde 683 713 
Acétaldéhyde 137 769 
Cétones Acétone 782 812 
Butanone 796 827 
Acides carboxylique Acide formique 710 742 
Acide acétique 753 784 
Acide propanoïque 766 797 
Esters Acétate de méthyle 791 822 
Nitriles Acétonitrile 748 119 
Amines Méthylamine 865 899 
Aniline 851 883 


IL.3.2 Transfert de proton vers et depuis le carbone 


L'une des réactions en solution les plus simples qui peuvent être étudiées cinétiquement 
est l'ajout ou le retrait d'un proton du carbone. Cela peut se produire de deux manières. Les 


hydrocarbures substitués qui sont suffisamment acides peuvent perdre un proton pour former 
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un carbanion (IL.15) ou un système insaturé peut être protoné pour donner un ion carbonium 


(IL.16). 


Ne , Bat (1.15) 


N 
7C—H + B = 


(L.16) 


N 7 N 7 
Æ=X. + BH JC —CTH + B 


Ces réactions couvrent une gamme inhabituellement large de réactivité et les transferts 
de protons à partir d’un acide aussi faible que le toluène (pKa environ 41) ou aussi fort que 
l’ion benzene-1-ylium (pKa environ 23) ont été étudiés. Cette revue concerne les taux de 
réaction (IL.15) et (IL.16) en solution. 

L'intérêt du transfert de protons vers et depuis le carbone découle en partie du fait que 
ce processus se produit comme une étape élémentaire dans les mécanismes d’un nombre de 
réactions importantes. Des réactions catalysées par un acide ou une base se produisent souvent 
à travers des ions carbonium intermédiaires ou des carbanions qui sont produits selon les 
réactions (IL.15) et (11.16). Connaître les propriétés acido-basiques des ions carbonium ou les 
carbanions peuvent également aider à comprendre les réactions dans lesquelles ces espèces sont 
présentes en tant qu’intermédiaires réactifs, même lorsqu'elles sont générées par des processus 
autres que le transfert de protons. 

Les études cinétiques des réactions simples telles que le transfert de protons sont 
également importantes dans le développement des théories de la cinétique. Etant donné que les 
taux et l’équilibre a peuvent souvent être mesurées pour (IL.15) et (IL.16), ces réactions ont été 
utiles dans l’étude des corrélations entre les coefficients de taux et constantes d’équilibre 
(relations d’énergie libre linéaires). 

Les réactions de transfert de protons vers et depuis le carbone ont également fourni des 
données pour tester les théories des effets cinétique des isotopes. La grande différence de masse 
entre l'hydrogène et le deutérium signifie que de grands effets cinétique sont observé avec ces 
isotopes et le transfert de protons vers et depuis le carbone est donc une réaction appropriée 
pour mesurer ces effets. 

Divers aspects du transfert de protons vers et depuis le carbone ont été passés en revue 
précédemment. Le travail classique de Bell [79] couvre l’ensemble du domaine des transferts 
de protons et discute en détail les mécanismes des réactions catalysées par acide et une base. 
La structure et la stabilisation carbanions ont été décrites [80] et une revue générale de 
l’ionisation des acides carbonés a été publiée [81]. Le rôle des transferts de protons dans les 
mécanismes des réactions chimiques et biologiques a également été décrit [82]. 
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IL.3.2.1 Méthodes expérimentales 


Des transferts de protons vers et depuis le carbone ont été observés dont la vitesse varie 
de ces réactions se produisant en moins d’une microseconde à celles qui prennent plusieurs 
années pour s’achever, couvrant ainsi presque toute la gamme des vitesses expérimentalement 
accessibles pour les réactions en solution. Un large éventail de techniques a été appliqué au 
problème du suivi des transferts de protons vers et depuis le carbone. Une sélection des 
méthodes les plus utilisées ces dernières années est décrite ici. Les mesures ne sont pas limitées 
aux conditions dans lesquelles l’acide carbonique est sensiblement dissocié. Les vitesses de 
transfert de protons pour des acides très faibles ou très forts dans des solutions aqueuses peuvent 
être obtenues avec des techniques tout à fait classiques. 

La discussion des méthodes expérimentales est divisée en deux sections : Sect. 
IL3.2.1.1, où des mesures directes de transfert de protons telles que (IL.15) et (IL.16) sont 
décrites et Sect. IL.3.2.1.2, qui concerne les mesures cinétiques d’une réaction globale avec un 


mécanisme dans lequel le transfert de protons est l’étape déterminant la vitesse. 


IL.3.2.1.1 Mesure directe des taux de transfert de proton 


Lorsque les réactions (IL. 15) ou (11.16) se produisent sur une échelle de temps de minutes 
ou d’heures, une méthode conventionnelle peut être utilisée pour suivre l’apparition du produit 
ou la disparition du réactif. L’ionisation des nitroparaffines en solution aqueuse montrée dans 


l’équation (11.17). 


RCHNO> + B —»> RCH:NO>: + BH* (L.17) 


Où B peut être OH, H20 ou une espèce tampon, a été suivi de mesurer le changement de pH 
[83], Conductance [84], ou absorbance UV [85] qui se produit lors de l’ionisation. 

Les transferts de protons comme (IL.15) et (IL.16) sont souvent trop rapides pour 
permettre une mesure utilisant une technique cinétique conventionnelle. Dans un nombre limité 
d’exemples, l’élargissement de raie RMN a été utilisé pour des réactions avec une demi-vie 
(t12) d'environ 107? sec. 

Deux protons distincts (disons un dans le solvant et l’autre dans l’acide carbonique) 
donneront deux raies dans le spectre RMN s1 la fréquence à laquelle ils s’échangent est petite 
par rapport à la différence entre leurs fréquences d’absorption. Si l’échange est rendu plus 
rapide (par exemple, en augmentant la concentration d’un catalyseur), les lignes fusionneront 
en une large ligne qui s’accentuera à mesure que la fréquence d’échange augmentera. Un 
comportement similaire sera observé si un proton échangeant se couple avec d’autres protons 
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de la molécule pour produire un multiplet. Si le taux d’échange du proton avec le solvant devient 
grand par rapport à la constante de couplage, le multiplet s’effondrera en un singulet. 

L'utilisation de ces deux approches a été faite dans une étude de l’échange de protons 
d’acétylène avec l’ion hydroxyde dans un mélange eau/alcool et une description complète de 
la technique est donnée [86]. Le transfert de protons catalysé par des ions hydroxydes à partir 
de 1,1-biséthylsufonyl éthane, (CH, CH, SO) CHMe, a été observé en solution aqueuse en 
utilisant la seconde approche [87]. 

Cependant, cette technique a été beaucoup plus largement utilisée dans les études sur 
l’échange protonique d’acides oxygénés et azotés [88]. 

Les méthodes de flux arrêté et de flux continu [89] ont été utilisées pour suivre les 
réactions de transfert de protons avec des demi-vies de l’ordre de la milliseconde. La méthode 
à flux arrêté qui est plus populaire est essentiellement un dispositif pour mélanger les réactifs 
rapidement (généralement en une milliseconde) avec certains moyens d’observer la réaction 
rapide qui suit. Le transfert de protons du cyanure de p-nitrobenzyle à l’ion éthoxyde dans des 
mélanges éthanol/éther à -77 °C a été étudié de cette manière [90]. La réaction a été suivie par 
spectrophotométrie. 

La réaction la plus rapide s’est produite avec t12 de 2 X10°?S bien que l’équipement soit 
adapté pour suivre les réactions avec t12 environ 2 X10*S. Une méthode similaire a été utilisée 
pour mesurer les taux de transfert de protons entre les acides carbonés faibles (par exemple, le 
triphénylméthane) et les bases (par exemple, les ions alcoxyde) dans le diméthylsulfoxyde [911]. 
Un appareil à flux continu avec détection spectrophotométrique a été utilisé [92] pour mesurer 
les taux d’ionisation des azulénes substitués en solution aqueuse. Réactions (IL.18) pour 


lesquelles des demi-vies comprises entre 2 et 70 msec ont été observées. 


HQ_H 
© 


+ HO  (ILI8) 


Les réactions se produisant plus rapidement que dans la gamme des millisecondes ne 
peuvent pas être suivies en utilisant des techniques d’écoulement en raison du temps fini 
nécessaire pour mélanger les solutions de réactifs. Ce problème est résolu de manière élégante 
en utilisant des méthodes de relaxation chimique [93]. Les réactifs et les produits sont autorisés 
à s’équilibrer et une perturbation est appliquée à la solution de réaction qui perturbe l’équilibre 


chimique. 
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Une perturbation en une seule étape peut être appliquée (par exemple, un changement 
soudain de température ou de pression) comme dans les méthodes transitoires ou une 
perturbation périodique peut être utilisée (par exemple, un champ électrique oscillant ou une 
onde ultrasonore) comme dans les méthodes stationnaires. Pour les réactions de transfert de 
protons en solution, une seule élévation progressive de la température se produisant en quelques 
microsecondes ou moins a été utilisée (méthode du saut de température). Une élévation de 
température de quelques degrés peut être obtenue en déchargeant un condensateur à travers la 
solution qui doit donc être électriquement conductrice (par exemple, une solution aqueuse ou 
non aqueuse contenant des sels dissous). 

Après le saut de température, la réaction chimique vers une nouvelle position d’équilibre 
peut être observée en surveillant le changement d’absorbance UV ou visible de la solution. Les 
conditions dans lesquelles la réaction est étudiée doivent généralement être telles que des 
concentrations finies à la fois de réactif et de produit soient présentes. Seuls quelques transferts 
de protons à partir du carbone ont été étudiés à l’aide de cette technique, en grande partie parce 
que ces réactions peuvent généralement être étudiées à l’aide de techniques conventionnelles. 
De plus, ce n’est que récemment que la fabrication commerciale d’instruments à saut de 
température a permis d’étudier des réactions de l’ordre de la microseconde [94]. Transfert de 
protons de l’acétylacétone (IL. 19) vers un grand nombre de bases (B) a été étudié [95] et d’autres 


exemples comprennent l’ionisation de la diacétylacétone [96] et du brormomalononitrile [97]. 


| 
H;C-C—CH—C—CH, + By* (19) 


Î Ï 
H;C—C—CH—C—CH + B 


ES 


Les études de transfert de protons à partir du carbone avec ti2 moins d’une 
microseconde sont assez rares. L’ionisation du cyanure d’hydrogène dans une solution alcaline 
aqueuse a été étudiée [98] en utilisant la méthode stationnaire ultrasonique qui est applicable 
aux réactions avec t12 dans la gamme de 10% à 10°? secs. 

Les études de transfert de protons vers et depuis le carbone impliquant la mesure directe 
de réactions comme (IL.15) et (11.16) sont moins courantes que les études d’une réaction globale 
dans laquelle le transfert de protons est une étape. C’est à ce dernier type d’étude que nous nous 


intéressons maintenant. 


IL.3.2.1.2 Etapes de transfert de proton déterminant la vitesse 


De nombreuses réactions en chimie organique montrent une catalyse acide ou basique 


générale et ce sujet a été bien étudié [79-82]. La catalyse acide générale signifie que la catalyse 
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à la fois par l’ion hydronium H30* et les acides tampons peut être observée et la catalyse 
basique générale fait référence à la catalyse par l’ion hydroxyde et les bases tampons. Une 
catalyse acide ou basique spécifique est utilisée pour décrire les cas où la seule catalyse 
détectable est due à l’ion hydronium ou à l’ion hydroxyde, respectivement. 

Dans de nombreux exemples où la catalyse générale est observée, il a été démontré que 
le mécanisme implique un transfert de protons déterminant la vitesse vers ou depuis le carbone 
et dans de nombreux cas, le coefficient de vitesse pour la réaction globale est identique au 
coefficient de vitesse pour l’étape de transfert de protons. Par conséquent, les coefficients de 
vitesse pour un simple transfert de protons vers ou depuis le carbone peuvent être obtenus en 
observant la réaction globale et quelques exemples de ceci sont décrits ci-dessous. Lorsqu’une 
catalyse spécifique est observée, le mécanisme généralement n’implique pas un seul transfert 
de protons déterminant la vitesse. 

Un mécanisme consistant en une protonation ou une dissociation rapide à l’équilibre 
d’une molécule de substrat, suivie d’une ou plusieurs étapes lentes qui n’implique pas de 
catalyseurs acides ou basiques, donne lieu à une catalyse spécifique et alors le coefficient de 
vitesse global ne se réfère pas à la vitesse d’un simple transfert de protons. Même lorsqu’une 
catalyse générale est observée, il faut veiller à ce que l’information mécaniste soit sans 
ambigüité. Par exemple, deux mécanismes peuvent être écrits pour une réaction montrant une 
catalyse acide générale. 

Ceux-ci sont représentés dans les réactions (11.20) et (1121) dans lesquels S se 


décompose aux produits en présence d’un acide général (BH) et d’une base (B). 


Rext 
S+BH*——> SH* + B Lente 
SH* — Produits Rapide (11.20) 


Taux = — d [S]/ dt = d [Produits] / dt = #gy+ [S] [BH*] 


S+BH*< SH* + B Rapide 

SH*+ B — Produits Lente (1.21) 
K =[SH*][B]/[S][BH*] 
Taux = — d [S]/ dt = d [Produits] / dt = # K [S] [BH*] 


Pour le mécanisme (11.20), le coefficient de vitesse observé (#gn+) est le coefficient de 
vitesse pour une étape de transfert de proton, alors que dans (11.21), le coefficient de vitesse 
observé ne se réfère pas à une seule étape de transfert de proton. Les exemples suivants illustrent 


les manières dont les coefficients de vitesse pour un simple transfert de protons vers ou depuis 
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le carbone ont été obtenus en mesurant la vitesse d’une réaction globale composée d’un certain 
nombre d’étapes. 


i. Echange isotopique de l’hydrogène 


Avec le développement des compteurs à scintillation liquide qui permettent de mesurer 
avec précision les niveaux de traceur de tritium, les réactions d’échange d’isotopes sont 
maintenant généralement effectuées avec du tritium plutôt qu’avec du deutérium. Cependant, il 
est souvent utile de comparer les taux d’échange pour les trois isotopes et l’échange de 
deutérium peut être suivi en utilisant, par exemple, la spectrométrie de masse [99]. Un exemple 
de l’application de l’échange d’isotopes au transfert de protons à partir du carbone est montré 
en équations (IL22) et (1123) pour l’échange catalysé par des ions hydroxydes de 


phénylacétylène en solution aqueuse [1001]. 


hi 
PhC=CT + OH PhC=C: + HTO (IL.22) 
E—————— 
hk-1Sec" 1 
ASec71 
PhC=C + H20 —> PhC=CH + OH (IL23) 


Etant donné que la réaction est effectuée dans une plage de pH où le phénylacétylène sera 
dissocié de manière négligeable (par exemple, à pH 8 en supposant une valeur pk de 21 pour le 
phénylacétylène [101] environ 1 partie sur 10! sera dissociée) le carbanion intermédiaire 
n’atteint que des concentrations à l’état d’équilibre. Etant donné que le marqueur tritium est 
présent en quantités de traceur, on suppose également que #2 > #1 où #1 sont écrits sous forme 
de coefficients de taux du premier ordre y compris les termes de concentration pour HTO et 


H20. Dans ces conditions, l’expression du taux se réduit à : 
—d [PhC=CT] /dt = #; [OH] [PhC=CT] (IL.24) 


La réaction d’échange a été observée [100] en suivant la perte d’activité de premier 
ordre du phénylacétylène marqué au tritium en présence d’ion hydroxyde en excès par 
extraction de la solution aqueuse à des intervalles de temps et l’activité a ensuite été mesurée 
avec un compteur à scintillation liquide. La comparaison du taux d’échange des différents 
isotopes peut être obtenue en mesurant les taux de détritiation et de perte de deutérium dans 
l’eau (#7 et #?) ou en mesurant le taux de détritiation dans D20 et le taux d’absorption du 
deutérium de D20 dans des composés non marqués. Substrat (#T et #4 en D20) [99]. Dans des 
conditions où l’expression de vitesse (IL.24) s’applique (#2 © #1), la vitesse d’échange 


isotopique mesure le coefficient de vitesse (#1) de l’étape de transfert de proton (11.22). Taux 
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de transfert de protons pour un grand. Une variété d’acides carbonés a été étudiée en solution 
aqueuse en utilisant cette technique et de nombreux exemples seront discutés tout au long de 
cette revue. Les réactions sont généralement catalysées par une base générale et par conséquent, 
on suppose généralement que le taux d’échange isotopique mesure le taux de proton. Transfert 
de l’acide carbonique à la base générale. 

Dans certains cas, notamment dans les solvants à faible constante diélectrique ou la 
formation de paires d’ion est importante, l’hypothèse que #2 © #1 n’est pas toujours valable. 
L’intermédiaire formé par réaction entra l’acide carbonique (RCT) et la base (B) peut retourner 


aux réactifs à une vitesse comparable à la vitesse à laquelle 1l passe aux produits (1.25). 
ha . 
RCT+B * RC°...TB* + Echange (IL.25) 
hi 


On dit alors que la réaction implique un retour interne [102]. Dans ce cas, l’expression du 
taux est donnée par : 


—d [RCT]/ dt = {#1 #2/ (k1+ #2) }[RCT] [B] (11.26) 


ii. Halogénation des acides carboniques 


L’halogénation catalysée par une base générale d’acides carbonés implique souvent un 
transfert de proton à partir du carbone déterminant la vitesse, et l’halogénation est une 
technique qui est fréquemment utilisée pour obtenir des vitesses de transfert de protons. La 
technique a été développée notamment par Bell [79] pour mesurer les taux de transfert de 


proton à partir des cétones. D’autres acides carbonés ont cependant été étudiés, par exemple 


les esters [103], les sulfonates [104], les nitro paraffines [105] et cyanocarbures [106].dans 
des conditions choisies qui peuvent varier selon le substrat, le mécanisme illustré pour 


l’acétone en équations (11.27) et (11.28). 


O © 

RCE + B de HC-C=CH, + BH* Lente (1127) 
(e O 

HC-C=CH, + X  — H,C-Ù—cH,x pee RADEON) 


d [CH:. CO . CH2X] / dt = d [X] / dt = #8 [B] [CH3COCH:] 


La réaction est d’ordre zéro en halogène (X)) et de premier ordre en catalyseur basique. 
Les concentrations d’halogène doivent être maintenues assez élevées pour que l’étape (11.28) 
soit beaucoup plus rapide que l’inverse de la réaction (11.27) [106]. Si l’inverse de la réaction 
(IL.27) ne peut être négligé, la réaction ne sera pas d’ordre zéro dans l’halogène. 


a. 
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Les cas où les mécanismes (11.27) et (11.28) tiennent avec halogène et catalyseur basique 
en excès par rapport à l’acide carbonique, la disparition de l’halogène ou l’apparition du produit 
bromé sera premier ordre. Avec un excès d’acide carbonique par rapport à l’halogène, la 


disparition de l’halogène ou l’apparition du produit sera d’ordre zéro. 


Des réactions ont été étudiées dans ces différentes conditions et la disparition de 
l’halogène a été suivie [103] par spectrophotométrie ou potentiomètre. Les produits peuvent 
être suivis par titrage du dérivé bromé [107] ou à partir de la conductance de l’ion halogénure 
libéré [108]. L’halogénation des cétones est également généralement catalysée par un acide. 
Le mécanisme consiste généralement en une protonation rapide de pré-équilibre du groupe 
cabonyle suivi d’un lent transfert de protons du carbone vers le catalyseur basique [1091]. L’énol 


ainsi produit réagit rapidement avec l’halogène. 


Le mécanisme global est similaire au mécanisme (1121) décrit précédemment et le 
coefficient de vitesse observé est un produit de la constante d’équilibre pour protonation du 
groupe carbonyle et le coefficient de vitesse pour le transfert de proton à partir du carbone et 


donc ne se réfère pas à un seul proton étape. 
ii. Transfert de protons vers les oléfines 


Il existe un certain nombre de réactions qui impliquent un transfert lent de protons vers 
le carbone oléfinique comme première étape suivie d’autres étapes rapides qui dépendent de 


l’oléfine étudiée [110]. Le mécanisme est représenté dans l’équation (IL.29), à savoir. 


k,lente Rapide . 
S + HA ——> SH* —— Produits (11.29) 


Taux = —d{[S] /dt = d [Produits] / dt = À [S] [AH] 
L’hydratation des oléfines avec Ri= Ar, R2 = H [43] et R1 = R2= Me [44] suit ce mécanisme 
comme indiqué en (IL.30), à savoir. 
hk,lente Rapide 
R1R2C=CH + H30* ——> RR2C'CH3 —— RR2C(OH)CH; (1.30) 
—d [R1R2C=CHb] /dt = d [R1R2C(OH)CHä] /dt = À [[R1R2C=CHb] [H30*] 


Une catalyse acide générale serait prédite pour une réaction avec le mécanisme (11.30) 
puisque d’autres acides peuvent remplacer H30* dans la première étape. 

Cependant, cette catalyse n’est pas toujours observable expérimentalement car elle peut 
être masquée par la catalyse par l’ion hydro. La catalyse acide générale a été détectée dans 


l’hydratation du p-méthoxy-a-méthylstyrène (R1= p-MeOC6H4, R2 =H) dans des solutions 


Cu 
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tampons de formiate [113] mais n’a pas pu être détectée pour l’hydratation de l’isobutène (R1= 
R2 = Me) [112]. Une légère catalyse générale était détectable dans l’hydratation du trans- 
cyclooctène(Il) et du 2,3-diméthyle-2-butène(IID) dans des tampons aqueux H3PO4H2PO'4 et 
HSO‘4S0"4 [1141]. 

L'observation de la catalyse acide générale établit le mécanisme comme (11.29) et est 
incompatible avec un mécanisme pour l’hydratation des oléfines montrées dans (11.31) qui 
implique la formation rapide et réversible d’un complexe, suivie d’un effondrement déterminant 


la vitesse en un ion carbonium et d’une hydratation rapide [113], à savoir. 


H Me 
(I) (HD) 
H* 
\ / Rapide _N t / Lente | 7H O | | 
C=C + H—% CC > H—C—Ct > H—C—C—OH(H31 
PRN  . EN ne 


La catalyse acide générale est également clairement observable [115] dans l’hydratation du 


diméthylacétal de cyanocétène et d’autres acétals de cétène (11.32), à savoir. 


NC /0Me 
ETC + H0 ——# NC.CH,.CO;Me + MeOH (IL32) 
H OMe 


Dans l’hydratation de l’isobutène, un mécanisme alternatif (1133) est écartée en 
montrant que lorsque la réaction est effectuée dans D20, une seule liaison carbone-hydrogène 


dans l’alcool t-butylique est produit était deutérée [112]. 
Me2C=CH2 + H30* + MeOH = Me:C* + H20 + CH3OH Rapide (1.33) 


Me:C* + H50 — Me3COH + H* Lente 


Cette découverte est compatible avec le mécanisme (11.34). 


: Lente > D,0 Rapide 
Me=CH2 + D30* —> MeC*CH2D ———— Me2CO(OD)CH:D (1134) 


Les preuves du fonctionnement du mécanisme (11.34) dans l’hydratation des oléfines et d’autres 


réactions apparentées ont été examinées [110]. 
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III. Méthodes théoriques de la modélisation moléculaire 


La modélisation moléculaire est une technique permettant, non seulement de représenter 
les propriétés, les forces ainsi que les réactions chimiques mais aussi de représenter les modèles 
des structures en deux ou trois dimensions. L'important développement au cours de ces 
dernières années des moyens informatiques a permis à la chimie de s'enrichir d'outils 


informatiques spécialement dédiés à la représentation des molecules. 


L'utilisation d'un modèle pour la prédiction des propriétés d'un système moléculaire, 
nécessitent des hypothèses ainsi que des approximations dont se sert la modélisation 
moléculaire et ceci en utilisant des méthodes de calculs théoriques telles que la mécanique 
quantique (QM) (ab initio ou semi empirique), la mécanique moléculaire (MM) et la dynamique 
moléculaire (DM). 

La mécanique quantique s'est imposée de nos jours, comme une théorie fondamentale 
permettant l'interprétation de maints phénomènes du monde matériel. Cela est dû au fait que 
la mécanique newtonienne, c'est à dire l'application de l'équation fondamentale de la 
dynamique, ne peut décrire correctement le mouvement des particules microscopiques. 

Les calculs quantiques représentent des outils de plus en plus utiles et fiables pour 
l'étude de la structure et de la réactivité des composés chimiques [1,2]. 

Comme toute discipline, la chimie quantique se base sur des postulats et sur des 
équations. Leurs résolution permet la compréhension de phénomènes qui, par voie 
expérimentale, seraient très difficiles, voir impossible. L'un des principes fondamentaux de la 


mécanique quantique repose sur la résolution de l'équation de Schrôdinger [3] 


IIL.1 Equation de Schrôdinger 


L'état d’un système à N noyaux et n électrons est d’écrit en mécanique quantique par 


une fonction d’onde Y satisfaisant l’équation de Schrôdinger [3]. 


HY =EY (IL. 1) 
Y : fonctions propres de H. 
E : valeurs propres de H. 
L'hamiltonien H total d’une molécule comportant N noyaux et n électrons, est défini par 
la somme de cinq termes, dont celui de la cinétique des électrons, de la cinétique des noyaux, 
des répulsions électrons-électrons, des répulsions noyaux-noyaux et celui des attractions 


electrons-noyaux. 
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L'approximation de Born-Oppenheimer [4] a été avancée en 1927 par Max Born et Robert 
Oppenheimer. Elle consiste à séparer les mouvements des noyaux et ceux des électrons. Elle se 
base sur le fait que la masse du noyau est beaucoup plus grande que celle d'un électron 
(mp/meT1836), alors que les forces agissant sur ces particules sont du même ordre. Les électrons 
vont donc s'adapter instantanément à tout changement de positions des noyaux. Les 


mouvements de chaque particule pourront ainsi être étudiés séparément. L'hamiltonien H peut 


se réduire à la forme suivante: 


2 n N nn Z e? n e° 
rs SA, D Te > (IIL.3) 
2m, = Kai Ry il 


La résolution exacte de l’équation (1) n’est possible que pour l’atome d’hydrogène et 
les systèmes hydrogénoïdes. Pour les systèmes poly-électroniques, il est nécessaire de faire 
appel aux méthodes d’approximation pour résoudre l’équation de Schrôdinger d’une manière 
approchée. 

Les propriétés moléculaires qui peuvent être calculées par la résolution de l'équation de 
Schrôüdinger sont multiples. On peut citer entre autres : 
Ÿ”_ Structures et énergies moléculaires 
Energies et structures des états de transition 
Fréquences de vibration, spectres IR et Raman 
Propriétés thermochimiques 
Energies de liaison 
Chemins réactionnels 
Orbitales moléculaires 
Charges atomiques 
Moments multipolaires 
Déplacements chimiques RMN et susceptibilités magnétiques 
Affinités électroniques et potentiels d’ionisation 


Polarisabilités et hyperpolarisabilités 


D 


Potentiels électrostatiques et densités électroniques 


La résolution exacte de l’équation de Schrüdinger n'est possible que pour l'atome de 
l'hydrogène et les systèmes mono-électroniques. Pour les systèmes poly électroniques, on fait 


appel aux méthodes d'approximation, pour la résolution de cette équation. 
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e Les modèles semi-empiriques. 
e Les méthodes ab-initio (Hartree-Fock). 


e Les approximations de champ moyen (DFT). 


IIL.2 Les méthodes semi-empiriques 


Les méthodes semi-empiriques [5] de chimie quantique sont basées sur le formalisme 
de Hartree-Fock , mais certains paramètres sont introduits pour simplifierles calculs. Ces 
paramètres proviennent, soit de données expérimentales, soit decalculs ab initio. Les méthodes 


| 


semi-empiriques permettent de modéliser les “ gros ” systèmes moléculaires grâce à deux 
approximations: La première consiste à ne prendre en compte que la couche de valence. La 
seconde consiste à négliger la plupartdes intégrales de répulsion électronique à plusieurs 
centres. On peut donc considérer les méthodes semi-empiriques comme intermédiaires entre les 
méthodes quantiques les plus sophistiquées et la mécanique moléculaire. L'avantage des 
méthodes semi-empiriques par rapport aux méthodes quantiques ab initio est la diminution des 
temps de calcul, car le nombre d’intégrales à calculer est réduit. Plusieurs méthodes semi- 
empiriques sont disponibles avec un degré de précision et une rapidité différents Pople, 1970. 
Pour nos calculs théoriques, nous avons choisi la méthode PM3 (parametric method 3), car 
de nombreuses études ont montré qu’elle donnait une meilleure descriptiondes géométries 
et propriétés électroniques moléculaires ainsi que les interactionsintermoléculaires dans les 
complexes d’inclusions et de transfert de charge. 

Et nous avons utilisé aussi la méthode nommée PM6 pour obtenir les minimas globaux ou 
locaux correspondent à des géométries « stables » de ces complexes. 

La méthode PM6 a été proposée par Stewart [6]. Et contrairement aux autres méthodes semi- 
empiriques, elle a adopté les paramètres d’interaction diatomiquenoyau-noyau de Voityuk’s [7] 
et utilise l’approximation de Thiel [8] en ajoutant des orbitales «d» des éléments hypervalents 


du groupe principal de la classification périodique. 


IIL3 Les méthodes ab initio 


Elles ne font appel à aucun paramètre ajusté à l’expérience. Les calculs sont 
généralement plus complexes, mais cette complexité ne signifie pas que les résultats soient 
exacts et non critiquables. En effet, ces méthodes utilisent toutes des approximations. Les 
calculs ab initio proviennent soit des méthodes d’Hartree Fock (et post Hartree Fock), soit de 
celles de la théorie de la fonctionnelle de la densité. Ces deux approches diffèrent dans la 


résolution de l’équation de Schrüdinger. 
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Contrairement à la méthode d’Hartree Fock, la description du système quantique dansle cadre 
de la théorie de la fonctionnelle de la densité est basée non pas sur la fonctiond’onde mais sur 


la densité électronique p(r). 


III.4 Théorie de la fonctionnelle de densité (DFT) 


La théorie de la fonctionnelle de la densité (DET) est actuellement l'une des méthodes 
les plus efficaces pour traiter les systèmes chimiques au niveau ab initio. Plusieurs raisons 


majeures ont contribué à sa popularité pour les applications chimiques: 


* Cette théorie inclut dans son formalisme une grande part de la corrélation électronique. 
*__ La méthode peut être appliquée à tout type de système: covalent, ionique ou métallique. 
+ Les ressources informatiques nécessaires sont moins importantes que pour les calculs 
"post-HF"', rendant ainsi accessible les études de systèmes moléculairesde plus grande 
taille. 
+ __L’aspect mono déterminant permet une interprétation « chimique » de la fonctiond'onde 
issue de ce formalisme. 
L'idée principale de la DFT est ancienne. Il a fallu attendre les années soixante pour 
qu’apparaisse la théorie de la fonctionnelle de la densité, basée sur les théorèmes de Hohenberg 
et Kohn ( Hohenberg, 1964 [9] ). Les premiers travaux permettant d’exprimer l’énergie d’un 


système en fonction de sa seule densité électronique p(r) [10]. 


IILS Différents types de fonctionnelles 
IIL.5.1 Fonctionnelles d’échange-corrélation 


Kohn et Sham ont en décrit 1965 une procédure de résolution d’un calcul DFT très 
similaire à la méthode HF. Les électrons sont ici considérés comme des particules 
indépendantes qui n’interagissent pas entre elles et les équations sont résolues de manière auto- 
cohérente [11]. L'énergie du système s’écrit comme la somme de quatre termes, dont trois sont 
calculables de manière exacte, à savoir, l’énergie cinétique des électrons, l’énergie potentielle 
d’attraction électrons-noyaux et une partie de l’énergie de répulsion électron-électron. Le 
dernier terme rassemble les corrections au modèle ainsi que l’ensemble des termes qu’on ne 
peut pas calculer sans approximation ; il est appelé énergie d’échange-corrélation et on la note 
Exc. Les équations de la DET étant similaires à celles des calculs HF (Hartree-Fock), les temps 
de résolution sont proches. Cependant, dans un calcul DFT la corrélation est prise en compte 


alors qu’elle ne l’est pas dans un calcul HF. Pour le même temps de calcul, on aura donc un 
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résultat meilleur. C’est pour cette raison que la DFT est autant utilisée à l’heure actuelle. 


Précisons tout de même que certains problèmes ne peuvent être résolus par la DET [12]. 


IIL.5.2 Différents niveaux d’approximation 


Le terme d’échange-corrélation est souvent décomposé en deux termes : un terme 
d’échange électronique et un terme de corrélation (on note alors Exc = Ex + Ec). Lors de 
l’utilisation de la DFT, le problème principal est de trouver la fonctionnelle d’échange 
corrélation la mieux adaptée au problème étudié. Il existe en effet de nombreuses 
fonctionnelles, calibrées sur un type de molécules données et avec des niveaux 
d’approximations différents. Selon le type de système considéré et de résultat recherché 
(thermochimie, spectroscopie, réactivité..etc.), on s’orientera vers l’une ou vers l’autre des 
fonctionnelles. On rencontre plusieurs niveaux d’approximation dans l’écriture des 
fonctionnelles d’échange-corrélation : 

Le premier niveau part du modèle du gaz uniforme d’électrons et consiste à ne considérer que 
la densité électronique locale (et pas ses dérivés). Il se nomme LDA (pour Approximation de 
la Densité Locale) ; la fin des années 1980 a vu le développement des fonctionnelles dites GGA 
(pour Approximation du Gradient Généralisé) dans lesquelles les fonctionnelles dépendent non 
seulement de la densité mais aussi de son gradient (dérivée première de la densité par rapport à 
l’espace) ; le niveau d’amélioration suivant consiste à utiliser en plus de la densité et de son 
gradient, le laplacien de la densité (dérivée seconde de la densité par rapport à l’espace) [13,14]. 
Pour des raisons techniques, il a été préféré à ces méthodes l’utilisation de la densité d’énergie 
cinétique des électrons qu’on nomme fonctionnelles méta-GGA. Considérer des dérivées 
troisièmes de la densité n’a pas d’intérêt puisqu’une telle fonction n’a pas de signification 
physique. Les fonctionnelles développées aujourd’hui cherchent principalement à prendre en 
compte de façon satisfaisante les interactions de van der Waals. Pour cela plusieurs voies sont 
suivies: il est possible d’ajouter un terme empirique de dispersion [15], de décomposer le terme 
d'interactions inter-électronique selon la distance entre électrons [16] ou encore d’introduire 
une part de résultat issue de calculs post-HF dans les fonctionnelles [17]. Bien évidement, si 


une complexité élevée apporte un meilleur résultat, elle nécessitera un calcul plus long. 


IIL.S5.3 Fonctionnelles hybrides 


En plus de ces trois niveaux d’approximation une autre famille de fonctionnelles est 
apparue, les fonctionnelles hybrides [18, 19]. Il s’agit ici de décomposer l’énergie d'échange 


en une partie calculée de manière similaire à la méthode Hartree-Fock, qu’on appelle échange 
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exact, et une partie issue des fonctionnelles précédemment citées. Parmi les premières 
formulations efficaces on trouve une écriture de la fonctionnelle d’échange corrélation sous la 
forme : EXC = aEXHF + bEXC DFT où a et b sont fixés [20]. Cependant si les énergies relatives 
obtenues sont correctes, les énergies absolues sont très mauvaises. Becke a donc proposé une 
formulation plus élaborée qui dépend de trois coefficients semi-empiriques pour donner lieu à 
la fonctionnelle B3PW91 dans laquelle la part d’échange exacte est de 20% [21] L’échange 
DFT vient de la fonctionnelle B88 proposée par Becke et ajustée sur l’énergie d’échange des 
gaz rares [22] alors que la corrélation est celle de la fonctionnelle PW91 [23, 22]. Les trois 
paramètres ont été choisis pour reproduire au mieux les énergies d’atomisation, les potentiels 
d’ionisation, les affinités protoniques et les énergies atomiques d’un ensemble de molécules. 
En remplaçant la fonctionnelle PW91I par LYP (développée par Lee, Yang et Parr et basée sur 
l’énergie de corrélation de l’hélium [24]), Stephens et coll. ont proposé la fonctionnelle hybride 
B3LYP qui est assurément la fonctionnelle la plus connue [25]. B3LYP ayant été abondamment 


utilisée en chimie organique. 


IIL.5.4 La fonctionnelle hybride B3LYP 


La fonctionnelle hybride B3LYP (Becke 3-paramètres Lee-Yang-Parr) est une 
fonctionnelle à trois paramètres combinant les fonctionnelles d’échange local, d’échange de 
Becke et d’échange HF (exact), avec les fonctionnelles de corrélation locale (VWN) et corrigée 
du gradient de Lee, Yang et Parr [26] La fonctionnelle d'échange-corrélation B3LYP [12, 13] 


est l’une des fonctionnelles les plus populaires: 


EXé =(1—a) ESA + a EST HD ESS + (1— CE + cEjyp (IL.4) 


ELyp C'est la partie de corrélation définie par Lee, Yang et Paar, EP58 


celle d'échange définie 
par Becke, et les paramètres a, b et c sont déterminés par ajustement par rapport à des données 
expérimentales. Généralement, les valeurs typiques sont a = 0,2 ; b = 0,7 et c = 0,8. 
L'introduction d'échange exact permet d'améliorer les résultats des calculs, même si la fraction 
d'échange à utiliser dépend du système étudié. Généralement, a = 0,2 est une valeur correcte 


pour de nombreux systèmes [27]. 


IIL.S.5 Fonctionnelles du Minnesota 


Les fonctionnelles GGA ont plusieurs défauts majeurs, en particulier elles sousestiment 
les barrières d’activation. Ceci peut être résolu en partie avec les fonctionnelles méta-GGA. De 


nouvelles fonctionnelles ont été développées récemment, du type hybride méta-GGA, qui 
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permettent de mieux résoudre ce type de problèmes. Parmi les plus populaires, on trouve les 
fonctionnelles de Zhao et Truhlar (MO05 [28], M06 [291, MO8 [30] et leurs dérivées). Elles ont 
été développées à l’université du Minnesota, d’où l’acronyme M, elles ont étés développées en 
2005, 2006 et 2008 respectivement. De l’aveu même de leurs auteurs, les fonctionnelles MOS 
sont à considérer comme une version préliminaire des fonctionnelles M06, qui sont ensuite été 
légèrement améliorées avec les fonctionnelles MO8 [31]. 

Outre les travaux de Zhao et Truhlar pour promouvoir leurs fonctionnelles, de nombreux 
articles de la littérature ont prouvés la qualité des résultats obtenus par les fonctionnelles M06 
[32,33]. Les fonctionnelles MO8 sont ainsi moins utilisées, bien qu’elles soient plus récentes. 
Leurs absences dans certains codes de chimie quantique est également une cause majeure de 
leurs assez faible utilisation. 

Au sein de la famille M06, on trouve plusieurs fonctionnelles (MO06-L [29], M06 
[33,34], M06-2X [33, 34], MOG6-HEF [35]), chacune étant plus spécifiquement adaptée à un 
certains types de problèmes (physique du solide, chimie organométallique, chimie organique, 
spectroscopie ...etc.). Si la fonctionnelle B3LYP dépend seulement de trois paramètres ; les 
fonctionnelles M06 dépendent d’environ de 40 paramètres (38, 39, 38, 35 respectivement). La 
part d’échange type Hartree-Fock est de 27% dans M06. Le 2X de la fonctionnelles M06-2X 
signifie qu’il y a deux fois plus d’échange de type Hartree-Fock que dans la fonctionnelle M06, 


ce paramètre a donc été fixé à 54% par Zhao et Truhlar, et tous les autres ont été recalibrés [31]. 


IIL.5.6 La série de fonctionnelles M06 


La série de fonctionnelles M06 est un ensemble de quatre fonctionnels métahybrides 
GGA. Ils sont construits avec un ajustement empirique de leurs paramètres, mais contraignant 
le gaz électronique uniforme. La famille comprend les fonctionnelles M06-L, M06, M06-2X et 
MO6-HF, avec une quantité différente d'échange exacte sur chacune d'elles. MO06-L est 
entièrement local sans échange HF, donc elle ne peut pas être considérée comme hybride. M06 
possède 27% d'échange HF, M06-2X contient 54% et MO6-HF comprend 100%. 

L’avantage de la fonctionnelle M06-L est dû au fait qu’est rapide, elle est bonne pour 
les calculs réalisés sur les métaux de transition, en chimie inorganiques et organométalliques. 
La fonctionnelle M06 est efficace pour la chimie organométalliques, cinétique chimique et 
l’étude des liaisons non covalentes. La fonctionnelle M06-2X est appropriée pour l’étude des 
éléments du groupe principal du tableau périodique des éléments ainsi que la cinétique 
chimique. En ce qui concerne la fonctionnelle MO6-HF elle est surtout commode pour les 


problèmes de transfert de charge. Les forces de dispersion ainsi estimées corrigent l'une des 
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plus grandes lacunes des méthodes DET. Le facteur d'ajustement relatif à la correction de 
dispersion à longue portée de Grimme est de 0,20, 0,25 et 0,06 pour M06-L, M06 et M06-2X 
respectivement. Alors qu’il prend les valeurs de 1.05 pour B3LYP et 0.70 pour PBEO [35]. 


IIL.6 Fonction de base 


Les orbitales atomiques et moléculaires résultent d’une combinaison linéaire de 
fonctions de base monoélectroniques, sachant qu’une base d’orbitales atomiques est définie 
pour chaque atome de la molécule. De surcroît, chaque orbitale (atomique ou moléculaire) peut 
être également développée sous la forme d'une combinaison linéaire de fonctions 


monoélectroniques 4, appelées fonctions de base: 


Les termes cui sont généralement désignés sous le nom de coefficients d'expansiondes 
orbitales moléculaires Pi. 


Il y a deux sortes de fonctions de bases qui sont d’un usage commun : 


IIL.6.1 Les fonctions de type Slater 


Ce sont des fonctions notées par STO (STO= Slater type orbitals en Anglais) [36] : sont 


caractérisées par une discontinuité à l'origine donnant une densité électronique correcte : 


S (6, r) = Nr" l exp (-Cr) Yem (0, d) (IIL.6) 
Les Y tm sont les harmoniques sphériques, N est une constante de normalisation, & l'exposant de 


l'orbitale, n et £ et m les nombres quantiques: principal, orbital et magnétique successivement, 


et r le vecteur position de l'électron en prenant comme origine la position du noyau. 


IIL.6.2 Les fonctions gaussiennes 


Ces fonctions ont été proposées par Boys (notée par GTO [37] (GTO= Gaussian type 
orbitals en Anglais) ; elles sont utiles pour la simplification du calcul des intégrales 


biélectroniques. La forme générale d'une gaussienne est : 

g (a, r) = NX" VU" Z'exp (-ar?) (IL. 7) 
Où a est une constante qui spécifie la taille ou extension spatiale de la fonction et Nest une 
constante de normalisation, Des combinaisons linéaires de primitives gaussiennes sont utilisées 


pour former les fonctions de base ®, et sont appelées gaussiennes contractées. 


L'utilisation des orbitales gaussiennes est cependant beaucoup plus répandue (logiciels 
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Gaussian 98, Games US, Molpro) car les calculs sont beaucoup plus faciles avec les fonctions 
gaussiennes et, de plus, une fonction Slater peut être correctement approximée en faisant une 
combinaison de fonctions gaussiennes. 

Pour notre étude, nous utiliserons des bases de gaussiennes, c’est-à-dire des bases qui 
utilisent des fonctions d'orbitales atomiques de type gaussiennes. En pratique, plus la base est 
grande et complète, plus flexible est la description des orbitales. En effet, moins de restrictions 
sont imposées aux positions des électrons dans l'espace. 

Lorsqu'une seule fonction de base est utilisée par orbitale atomique, la base est dite 
minimale. Si deux fonctions de base sont utilisées, la base est dite double -£ (DZ). On peut aussi 
ajouter des fonctions de polarisation, c’est-à-dire, des fonctions de base de nombre quantique 
de moment angulaire supérieur. 

De telles fonctions accroissent la flexibilité de la base en autorisant les charges 
électroniques a se délocaliser. De plus, on utilise des orbitales d'exposants petits qui diminuent 
progressivement lorsque le rayon augmente afin de décrire les régions de l'espace éloignées du 
noyau. Ces orbitales sont dites diffuses; elles sont importantes pour l'étude d'anions et d'états 


excités comme des états de Rydberg. 


Base minimale 


En pratique, les orbitales atomiques de type Slater (STO) sont approchées par une 
combinaison de plusieurs orbitales atomiques Gaussiennes (GTO). Les plus communément 
utilisées sont celles développées par Pople et collaborateurs [38]. Elles sont nommées STO-nG, 
encore appelées bases minimales, n représente le nombre de Gaussiennes primitives ajustées à 
une orbitale de type Slater (STO).La base minimale la plus simple est la base STO-3G, Ceci 
signifie que les orbitales de type Slater sont représentées par trois fonctions Gaussiennes 


primitives. 
Base double-zêta ou N-X1G 


Le niveau suivant développé par Pople comprend les bases de type valence séparée SV 
(SV= split-valence). La nomenclature N-X1G, N primitives sur les couches intérieures, X1 
nombre des primitives pour les contractions des couches de valence. Xest la notation pour les 
ensembles de bases de qualité double-zêta (DZ). 
Par exemple 3-21G, 6-31G, 6-311G donnent de bons résultats pour les calculs de géométrie et 
pour les propriétés des électrons de valence, tout en prenant en compte les électrons de cœur. 


Elles permettent donc de bien décrire la liaison chimique. 
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° La base 3-21G est une split valence-double zêta (SV-DZ), où chaque orbitale atomique 
des couches internes est décrite par une contraction de 3 gaussiennes primitives. Les 
orbitales de la couche de valence sont reparties en deux groupes: les orbitales proches du 
noyau sont décrites par une contraction de deux gaussiennes primitives, et les orbitales 


éloignées par une seule gaussienne primitive. 


+ La base 6-31G décrit une orbitale atomique par la combinaison de six fonctions 
gaussiennes par orbitale de cœur, de trois autres pour la description des électrons de 
valence et d’une dernière pour description des électrons de valence les plus éloignés du 
noyau (externe). 

* La base 6-311G est une split valence-triple zêta (SV-TZ) dans laquelle les orbitales de 
cœur (couches internes) sont exprimées par une contraction de 6 gaussiennes primitives. 
Les orbitales de la split couche de valence sont exprimées par des contractions de 3, 1 et 1 
primitives respectivement. Décrit une orbitale atomique par la combinaison de six 
fonctions gaussiennes par orbitale de cœur, de trois autres pour la description des électrons 
de valence et d’une dernière pour description des électronsde valence les plus éloignés du 


noyau (externe). 


Base polarisée 


Les bases polarisées contiennent des OA de polarisation qui permettent une description 
fine des systèmes plus ou moins polarisés, c’est à dire, des fonctions de base de nombre 
quantique de moment angulaire supérieur. De telles fonctions accroissent la flexibilité de la 
base en autorisant les charges électroniques à sedélocaliser, par exemple des systèmes 
ioniques. 

En pratique, des orbitales virtuelles sont ajoutées afin de prendre en compte la 
déformation des orbitales lorsqu'elles participent à une liaison. Par exemple, des orbitales p 
sont ajoutés pour l’atome d’hydrogène et des orbitales d sont ajoutés pour les atomes de la 
deuxième période de la classification périodique. 

Par exemple, dans la base 6-31G** des fonctions polarisées sont ajoutées sur les atomes 


lourds (première étoile) et sur les hydrogènes (deuxième étoile). La base 6- 31G** peut être 


notée également 6-31G (d, p). 


Base diffuse 


Les bases diffuses sont indispensables à inclure dans la base d’orbitale atomique chaque 


80 


Chapitre III Méthodes théoriques de la modélisation moléculaire 


fois que le phénomène physique décrivant la propriété étudiée nécessite une bonne description 
de l’espace situé au-delà des orbitales de valence (espace diffus).Ils sont utilisées pour décrire 
des régions éloignées des noyaux, qui augmentent la taille du nuage électronique. Pour les 
espèces ayant des doublets libres et les espèceschargées (anions), la présence d’orbitales 
diffuses sont nécessaires pour les systèmes ayant de faibles interactions intermoléculaires 
(exemple: liaisons hydrogène). Ces bases profitent aussi à la description des liaisons chimiques 


classiques, affinant la description de ces interactions. 


Ces bases sont désignées par un signe + (atomes lourds) ou ++ s1 la molécule possède des 


atomes d’hydrogène. 


IIL.7 Succès et limites de la DFT 


Le succès et la popularité de la DFT sont relativement récents [39]. W. Kohn a obtenu 
le prix Nobel pour ses travaux sur cette théorie en 1998. Depuis le début des années 90, le 
nombre de publications scientifiques dans différents domaines de la chimie et de la physique 
utilisant la DET a connu une ascension remarquable. 

Les développements récents utilisent un formalisme dépendant du temps (TD-DFT de l’anglais 


Time Dependant Density Functional Theory) [40] qui permet de décrire les états excités. 


IILS8 Etats excités - théorie de la fonctionnelle de la densité dépendante du temps (TD- 
DFT): 


Jusqu'ici la tâche principale des calculs de TD-DET à été de déterminer et d'interpréter 
des spectres électroniques d'absorption. Lorsque l’on veut étudier les états excités d’une 
molécule en DET, il n’est paspossible d’étendre simplement les équations de Kohn-Sham. Pour 
y remédier, 1l faut utiliser une autre approche qui est la théorie de la fonctionnelle de la densité 
dépendante du temps (TD-DET). 

Le théorème de Runge et Gross [41] établit de façon similaire au premier théorème de 
Hohenberg et Kohn, qu’il existe une relation directe entre la densité électronique dépendante 
du temps et le potentiel externe auquel il est soumis. Dès lors, le potentiel et la fonction d’onde 


sont dépendants de la densité. 


IIL.9 Modèles théoriques de solvatation 


Le rôle de l'eau est primordial pour la structuration et la fonction des édifices biomolé- 


culaires [42]. Pour tenir compte de ses effets en modélisation moléculaire, deux grandes 
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méthodes sont couramment utilisées. Dans la première, les molécules du solvant sont 
représentées de manière explicite à l'échelle microscopique, tandis que dans la seconde seuls 


les effets macroscopiques du solvant sont modélisés. 


IIL.9.1 Solvant explicite 


L'utilisation de molécules d'eau explicites est certainement la meilleure représentation 
du solvant disponible actuellement en modélisation moléculaire. Elle permet en effet de réaliser 
des simulations stables de dynamique moléculaire à l'échelle de plusieurs dizaines de 
nanosecondes. De plus, la représentation au niveau atomique de ces molécules rend possible 
l'observation de la formation de liaisons hydrogène entre le soluté et le solvant. Ainsi, 
l'utilisation d'un solvant explicite est indispensable lorsque l'eau est en interaction directe avec 


la biomolécule et que son rôle n'est pas limité à ses effets électrostatiques. 


IIL.9.2 Solvant implicite 


L'utilisation de modèles de solvant implicite est une alternative de plus en plus crédible 
à l'utilisation de solvants explicites. Dans ces modèles, le soluté est considéré comme une cavité 
de faible constante diélectrique &i au sein du solvant représenté comme un milieu continu de 
constante diélectrique &, élevée. Ils permettent notamment de dépasser certaines limitations des 
modèles de solvant explicite en étant plus avantageux en termes de temps de calcul et en 


permettant l'obtention de grandeurs thermodynamiques. 


IIL.10 Méthodes d'analyses de la densité électronique 
IIL.10.1 Orbitales Frontières HOMO-LUMO 


Lorsque deux atomes se combinent, deux orbitales atomiques donnent naissance à deux 
orbitales moléculaires : une orbitale d’énergie inférieure (orbitale liante) et une autre orbitale 
d’énergie supérieure (orbitale antiliante). Si na et ng sont les nombres des OM de A et 
B, le calcul des OM revient alors à traiter na ng problème à deux orbitales. 

Une première simplification peut être faite en remarquant qu’une interaction à 2 électrons est 
stabilisante et une interaction à 4 électrons est déstabilisante. Fukui (prix Nobel 1981) a alors 
introduit en 1952 une approximation [43] puisque plus les orbitales sont proches, plus leur 
interaction est forte, on ne considèrera que les interactions correspondantes aux deux paires 
d’OM occupées et vacantes les plus proches, et on négligera tout le reste. 

La plus haute OM occupée (HOMO) et la plus basse vacante (LUMO) sont appelées par Fukui 


orbitales frontières car elles séparent les orbitales occupées des orbitales vacantes (Figure IL. 1) 
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Figure IIL.1 Approximation des orbitales frontières HOMO-LUMO. 


Le problème devient encore plus simple dans le cas d’une réaction à fort caractère 
ionique. En effet, un réactif nucléophile (ou n’importe quel composé riche en électrons) donne 
facilement des électrons. 

Il réagit donc essentiellement par son HOMO), qui possède une énergie élevée. Inversement, un 
réactif électrophile (composé pauvre en électrons) tend à recevoir des électrons qu’il ne pourra 
loger que dans une orbitale vacante. 

Un électrophile réagit par conséquent essentiellement par sa LUMO qui a une énergie 
basse. La stabilisation dépend de l’écart énergitique HOMO-LUMDO et du recouvrement entre 


les orbitales en interaction . 


IIL.10.2 Les descripteurs de réactivité globaux (4, u, n, ©) 


La théorie de la fonctionnelle de densité (DET) constitue actuellement une vraie source 
de concepts chimiques comme le potentiel chimique électronique, l’électronégativité, la dureté, 
la mollesse, l’électrophilicité, ...etc. La DFT est fondée sur le principe variationnel. En effet, 
l’énergie d’un système est une fonctionnelle de la densité électronique. 

Les descripteurs globaux qui sont entre autre l’électronégativité (y), le potentiel 
chimique (L), la dureté chimique (n) et l'indice de L’électrophilie (@) ont été calculés en utilisant 
les énergies des orbitales frontières Exomo (comme l'énergie d'ionisation) et ELumo (comme 
l'affinité électronique). 

Ces descripteurs sont indépendants de l’espace, donc ils prennent la même valeur en 
tout point du système auquel ils se rapportent. Ils permettent donc de décrire la stabilité ou la 
réactivité globale d’une molécule. 

Le potentiel électrochimique (1) [44,45] qui est un descripteur global, est lié à 
l’électrophilicité et la nucléophilie d’une molécule. Il est l’inverse de l’électronégativité de 


Mulliken (}) [46]. 
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En tenant compte des énergies de la HOMO et de la LUMO correspondant au théorème 
de Koopmans [47], le potentiel électrochimique peut être en fonction de l’énergie du potentiel 


d’ionisation (PI) et de la première affinité électronique (AE) comme suit : 
1 1 
= [PI+AE] = ; [Exomo + ELumol = - % (IIL.8) 


La dureté globale (n), exprime la résistance de la molécule au transfert de charge qu’il 


s’agisse d’un gain ou d’une perte, elle mesure donc la stabilité de la molécule [48]. 
1 
n = : [ELumo- Exomol] (IL.9) 


La question posée est ainsi de savoir combien un électrophile pouvait acquérir 
d’électrons, s’il était immergé dans une mer d’électron libre. La valeur du nombre d’électrons 
acquis et la stabilisation énergétique qui en découlerait, permettrait de comparer le pouvoir 
électrophile de deux espèces [49,501. Soit donc un électrophile, la variation de son énergie 
électronique lors d’un transfert de charge à potentiel externe constant.La mollesse (s) est donnée 


par l’équation: 


— (IIL.10) 
n 


Des études expérimentales menées en 1998 [51,52], ont conduit à penser que l’indice 
d’électrophilicité (©) pouvait être un bon descripteur de la réactivité d’une molécule. 

L’électrophilicité (eq IL. 11) est définie par la capacité d’une molécule à se lier fortement 
à une entité nucléophile par un transfert d’électrons. Autrement-dit c’est la capacité d’acquérir 


des électrons pour se stabiliser. 


= — (IL. 11) 


IIL.10.3 L'analyse NBO (Natural Bond Orbital) 


L’idée originale de la méthode NBO est de proposer un découpage de la densité 
électronique d’un système moléculaire pour décrire la structure électronique de la molecule 
dans un schéma de type Lewis [53]. Ainsi la densité est décrite sur et entre les atomes avec des 
doublets de cœur (CR), des doublets non-liants (LP), des doublets de liaison (BD) résultant de 
la mise en commun par deux atomes de deux électrons dans des orbitales hybrides. 

La construction des NBO permet de déterminer une énergie ELewi de la fonction d’onde. 
La différence entre E et Ec.wi permet d’évaluer la qualité du schéma de Lewis proposé. Une 


analyse perturbatrice des interactions « donneuses-acceptrices », appelée analyse perturbatrice 
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au second-ordre est effectuée entre les NBO occupées (donneuses) type de Lewis et les NBO 
inoccupées (acceptrices), non-Lewis [54]. Elle traduit l’interaction « donneuse-acceptrice ». 
Elle mesure le transfert d’énergie engendré par la combinaison d’une NBO (1) occupée 
(donneuse) et une NBO (j) vacante (acceptrice). 

L'énergie de stabilisation E liée à la délocalisation de i vers j, est explicitement estimée 
par l'équation suivante [55]: 


F(j)? 
(2) — Per EE 
E AEÿ= qji cie] (L.12) 


Où qi est l'occupation de l’orbitale donneuse, F (ï, j) est l’opérateur de Fock et ei, #j sont les 


énergies des orbitales NBO [56]. 


IIL.10.4 L'analyse de populations de Mulliken 


La charge totale d’une molécule est naturellement répartie sur tous les atomes qui la 
constituent. Chaque atome est ainsi porteur d’une charge partielle, aussi appelée charge 
atomique, qui dépend de sa nature et de son environnement. Si la notion de charges partielles 
est très utile pour interpréter la réactivité d’une espèce chimique, elle n’en reste pas moins 
délicate à formuler d’un point de vue théorique puisque celles-ci ne correspondent à aucune 
grandeur observable. II existe bien des données expérimentales qui leur sont régulièrement 
associées, comme le déplacement chimique en RMN 13C ou le moment dipolaire, mais ces 
quantités ne donnent pas accès à des valeurs « exactes ». En revanche, il existe de nombreux 
modèles théoriques permettant de calculer les charges partielles. Tous ces modèles reposent sur 
une partition de la population électronique entre les atomes du système. Le choix du schéma de 
partition est crucial puisque c’est lui qui définit, la qualité de la distribution de charges obtenue. 

Parmi les méthodes de calcul des charges partielles, nous citons : 

+ Les méthodes Basées sur le Moment Dipolaire. 

* Les méthodes Basées sur la Densité Electronique. 

+ Les méthodes Basées sur les Orbitales Moléculaires, cette dernière est basée sur deux types 
d’analyses 

* Analyse de Population Naturelle (NPA). 

* Analyse de Population de Mulliken. 

En 1955, Mulliken propose une première approche pour le calcul de charges partielles. 
L’approche qu’il propose à l’époque se distingue par sa simplicité numérique et conceptuelle 
et reste l’une des plus utilisées aujourd’hui. La charge de Mulliken portée par l’atome 1 est 


définie comme la différence entre la population électronique de l’atome isolé (zi) et celle de 
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l’atome au sein de la molécule. 

Cette dernière est définie comme la somme des électrons situés dans les orbitales 
atomiques centrées sur 1 et la moitié des électrons situés dans les orbitales moléculaires 
construites à partir d’orbitales atomiques centrées sur 1. 


Sa formulation mathématique est donc : 
qi= ZilDuei(DS}4t + > Eveu(DS)UV] (I.13) 


Où: S: est la matrice des recouvrements, 
D: est la matrice densité et u: désigne une fonction de base quelconque (la somme des fonctions 
Lei correspond à la base décrivant l’atome i). 

Les résultats de l’analyse de Mulliken sont très dépendants des bases utilisées pour la 
description des orbitales atomiques. Mulliken lui-même reconnait que cette instabilité constitue 
une sérieuse limitation à sa méthode [57]. Une autre faiblesse de l’analyse de Mulliken provient 
de la répartition arbitraire des électrons situés dans les orbitales moléculaires. En effet, ces 
derniers sont attribués pour moitié à un atome et pour moitié au second, sans tenir compte de 
l’extension spatiale des orbitales moléculaires. Par conséquent, plus la base est étendue, moins 
l’analyse de Mulliken est fiable. Les composés présentant des liaisons fortement ioniques sont 
également mal décrits [58]. En revanche, l’analyse de Mulliken se montre peu dépendante de 
la corrélation électronique puisqu'elle conduit à des résultats similaires quel que soit la méthode 
de calcul utilisée (Hartree-Fock, B3LYP, MP2, QCISD) [591. 

IIL.10.5 L’analyse QTAIM (Quantum Theory of Atoms in Molecules) 


La théorie quantique des atomes dans les molécules (QTAIM) est une méthode d'analyse 
topologique de la densité de l'électron et de ses dérivées (Vp, V?p) développée par Bader [60]. 

La densité électronique peut être assimiler à un gaz à densité non uniforme, donc plus 
condensé dans une région de l'espace et dispersé dans une autre région de l'espace. La 
détermination de la forme où la topologie de la densité de charge volumique est assurée par un 
calcul d'un ensemble de points critiques dans lesquelles le gradient (Vp =0) s'annule, autrement 
(p=Constante). 


IIL.10.5.1 Gradient de la densité et Bassin atomique 


Pour déterminer la trajectoire du gradient de la densité électronique (Vp), on procédé 
delà manière suivante. On part d’un point quelconque de l’espace; en ce point le gradient est 
calcule et un “pas” selon la trajectoire vers la densité croissante est suivi. Ce processus ainsi 


répète avec différents points de l’espace, nous permet de tracer les trajectoires du gradient vers 
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tous les points critiques. Or les maxima de la densité électronique se trouvent au niveau des 
noyaux ; tout point de l’espace est donc relie a un des noyaux par le biais d’une trajectoire. Les 
noyaux sont alors appelés les attracteurs du gradient et l’ensemble des points dont la trajectoire 


finit sur un noyau forme le bassin atomique associe [611]. 


7 


Séparatrice 


Figure IIL.2 Exemple de deux bassins séparés par une séparatrice. 


Comme toutes les trajectoires du gradient convergent vers les noyaux, des bassins atomiques 
apparaissent dans le dessin du gradient (Figure IL.2). Ces bassins sont délimites par des surfaces 


interatomiques appelés séparatrices, qui satisfont l’équation : 
V p(r ). n(r )}=0 (IL. 14) 


Ou n(r) représente un vecteur normal a la surface. Cette procédure permet ainsi de réaliser un 
découpage de l’espace moléculaire en atomes. Ces derniers sont reliés par des liaisons 


caractérisées au moyen de leur point critique de liaison BCP. 


IIL.10.5.2 Points critiques 


Le découpage de l’espace moléculaire en régions disjointes repère plusieurs points 
stationnaires importants (maximum, minimum, point selle) dénomme Point Critique (CP pour 
critical point). 

Quel que soit le type de points critiques, ce sont toujours des points pour lesquels le gradient 
de la densité électronique est nul (Vp(rc) = 0). 

Pour classer les différents types de points critiques, il faut trouver un critère dédifférenciation. 
Ce critère est obtenu via le Laplacien de la densité électronique : si la dérivée seconde de la 


densité électronique en ce point est négative, alors il s’agit d’un maximum local; si au contraire 
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elle est positive, alors c'est un minimum local. Le calcul des dérivées secondes de p(r. ) serait 
donc nécessaire. 

En effet, considérant les dimensions (x, y, z), nous avons neuf dérivées secondes possibles qui 
sont représentées sous forme d’une matrice symétrique nommée, matrice Hessienne de p(rc ). 
Les points critiques sont alors distingues par un couple dealeurs (r, s) ou r représente le rang de 
la matrice Gessienne (nombre de valeurs propres non nulles) et s la signature (somme 
algébrique des signes des valeurs propres).Ces valeurs propres sont notées comme suit (A1, À, 
À3) avec comme convention (A1< À2 < À3.) 

Pour les structures moléculaires stables, tous les points critiques sont de rang 3. Dans ce cas, 
on trouve quatre différentes possibilités: 

*(3, +1) : trois valeurs propres positives => un minimum local de p(r ), il correspond a un point 
critique de cycle (RCP : Ring Critical Point) 

*(3, -3) : trois valeurs propres négatives => un maximum local de p(r ), il désigne la position 
du noyau atomique (NCP : Nuclear Critical Point) 

*(3, +3) : deux valeurs propres positives et une négative => 1l est un minimum selon deux 
coordonnées et un maximum selon une autre => un point selle, il correspond a un point critique 
de cage (CCP : Cage Critical Point) 

*(3, -1) : deux valeurs propres négatives et une positive => il est un minimum selon une 
coordonnée et un maximum selon les deux autres => un point selle, il correspond a un point 
critique de liaison (BCP : Bond Critical Point) 


La figure (IIL.3) illustre bien les différents types de points critiques CP cités précédemment. 


€ 
A 


, CE L. ca, +) 
\/ 


CP (3,+1) / 


Figure IIL3 (a) Carte du contour de la densité électronique, (b) Carte stéréoscopique 


(topologique), (c) graphe moléculaire du Benzène. 
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Pour une distribution électronique, le nombre de points critiques présent dans un système 


moléculaire doit vérifier la relation de Poincare-Hopf suivante : 
Nnce — Ngcp +Nrcp — Nccp = 1 (IL. 15) 
IIL.10.5.3 Nature de l’interaction interatomique 


Comme nous l’avons mentionne précédemment, l’analyse de la densité et de son 
gradient ne permet pas de donner des informations sur les appariements d’électrons. Untel 
manque a été résolu en considérant le Laplacien de la densité électronique V?p(r.) qui 


représente la somme des trois valeurs propres de la matrice Hessienne [62,63]. 


Vip(r.) =A + À + À (III. 16) 
La valeur du Laplacien en un point critique indique si la densité est localement concentrée 
(VZp(r.)<0), ou même diluée (VZp(r.)>0), toujours par rapport au voisinage de ce point. En 
général, les liaisons chimiques sont classées en deux catégories : 
*Interaction a partage d'électron (Shared Interaction) : le Laplacien de la densité est négatif 
au point BCP : VZp(r.) < 0. Ce type d’interaction se produit souvent entre les atomes lies par 
une “liaison covalente”. 
*Interaction a couche fermée (Closed Shell Interaction) : le Laplacien de la densité est positive 
au point BCP : VZp(r.) > 0. Ce type d’interaction correspond aux interactions “ioniques” et 
“van der Waals”.De nombreuses applications de la théorie AIM, ont contribue dans le 
développent d'autres relations importantes entre le Laplacien de la densité électronique aux PC 
et les paramètres topologiques énergétiques. L'une des relations principales est la forme locale 


du théorème du viriel [64,65] : 
=Vp(r)= 2 G(re) + V(re) (L.17) 


ou G(r.) et V(r.) représentent respectivement la densité d’énergie cinétique, potentielle. G(r.) 
est une quantité positive alors que V(r.) est négative ; le balancement entre la densité d'énergie 
cinétique G(r.) et la densité d'énergie potentielle V(r.) révèle la nature de l'interaction 
interatomique établisse entre les deux atomes impliques. 

Le théorème de Viriel, appliqué à un point critique de liaison BCP, implique des 
interactions ayant un Laplacien VZp(r.)<0 sont dominées par une réduction locale de l'énergie 


potentielle. Inversement, celles possédant une valeur positive du Laplacien (VZp(r.)> 0), sont 


dominées par un excès local de l'énergie cinétique [66]. Cet équilibre entre G(r.) et V(r.) 


89 


Chapitre III Méthodes théoriques de la modélisation moléculaire 


fournira des informations sur la déplétion ou la concentration de la densité électronique, et donc 
sur le type de liaison. 

Un autre critère a été également utilise pour caractériser l’interaction interatomique. Il 
s’agit du rapport “énergie cinétique/énergie potentielle” (G(r.)/ V(r.) [67]. Lorsque ce rapport 
est supérieur a 1, la nature de l'interaction est purement non covalente. 

Une autre alternative a été propose par Roza et al [68] pour évaluer la rigidité d’une interaction 
au point critique de liaison BCP (3, -1). Elle consiste a combiner le Laplacien de la densité 
électronique VZp(r) avec l’énergie électronique totale H(r). 

Les interactions type HBs pourraient être classées selon les critères suivants: 

+ Les liaisons HBs fortes sont caractérisées par VÆp(r) <0 et H(r) < 0 

+ Les liaisons HBs de force moyenne sont caractérisées par VZp(r) >0 et H(r) < 0 


+ Les liaisons HBs faibles sont définies par VZp(r) > 0 et H(r) > 0 


IIL.10.5.4 Notion de chemin de liaison (Bond Path) 


Bader définit des chemins de liaison, correspondant aux trajectoires reliant deux 
attracteurs et passant par un point critique appelé point critique de liaison BCP. Ce point 
correspond a un minimum local de la densité électronique suivant le chemin de liaison et a un 
maximum suivant le plan perpendiculaire. 

La présence d’un “bond path” ou un BCP entre deux atomes est la condition nécessaire 


et suffisante pour l’existence de liaison chimique entre deux atomes en question [69]. 


IIL.10.5.5 Ellipticité de la densité électronique 


L’ellipticité de liaison, représentée par #, est un autre paramètre intéressant qui fournit 


une mesure quantitative de l’anisotropie de la densité électronique au BCP. 
L’équation de l'ellipticité s’écrit comme suit : 
e=(u/h)-1 (III. 18) 


Ou À1 et A2 représente les courbatures dans le plan perpendiculaire a la direction de la liaison, 
ils déterminent la symétrie de la densité électronique au niveau du point BCP. 

*Si M = À, donc & = 0 => la liaison présente une symétrie cylindrique, c’est le cas des liaisons 
simples (C-C) dans l’éthane et triples (C—=C) dans l’acétylène. 

*Si >, donc € est élevée —> la liaison présente une symétrie elliptique, cas des liaisons 
doubles (C= C). 


Par conséquent, £ peut être utilisée comme un outil de mesure de la délocalisation. 
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IIL.10.6 Analyse des interactions non-covalentes (NCT) et méthode du gradient réduit de 


la densité de charge (RDG) 


L'analyse microscopique de toutes les propriétés physiques et électroniques d'un 
complexe de transfert de charge est basée essentiellement sur le comportement de l'ensemble 


de liaisons covalentes mais aussi sur les différentes interactions non covalentes. 


Par définition, une liaison covalente est la mise en commun des électrons de valence de deux 
atomes, cette liaison est très forte et présente une énergie de l'ordre de centaine de KkJ /mol [701]. 
Dans ce contexte, on peut parler aussi de la liaison ionique qui est généralement formée par la 
différence d'électronégativité entre atomes, c'est une liaison purement polarisée. 

Il existe dans la pratique une autre famille de liaisons qui est très différente que celles ioniques 
ou covalentes, il s'agit d'interactions fortes et faibles à effet liant et anti- liant dans le complexe. 
La nature de ces interactions est purement électrostatique. L’importance de cet effet dépend de 


l’aspect quantitatif d'atomes dans la structure. 


IIL.10.6.1 Moyen d'analyse et de visualisation des interactions non covalentes (NCI) 


Le principe d'analyse des différentes interactions non covalentes est relatif à l'ensemble 
des forces d'échanges entre ions et dipôles dans le système étudié. Généralement, la répartition 
de dipôles ainsi que de charges ponctuelles fait partie de la densité de charge volumique p qui 
présente la grandeur rationnelle de cette analyse [71,721]. 

Le moyen le plus adapté pour la quantification des forces (énergies) d’interactions non 
covalentes sera essentiellement basé sur les variations ainsi que la répartition de p dans l'espace 
du complexe. À titre d'exemple, pour décrire les forces de Van Der Waals, il suffit d'analyser 
les zones de faible densité et cela afin de visualiser et de quantifier approximativement les 


interactions faibles [73]. 


IIL.10.6.2 Principe de la NCI 


Le principe de cette méthode est basé sur l’exploitation du Gradient Réduit (RDG) de 
la densité de charge électronique noté s(p). Ce terme considère conjointement le gradient Vp 


ainsi que la valeur de la densité de charge p du système [74]. 


(p) = _ (IL. 19) 


2(37?)3p3 
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La visualisation de l'interaction non-covalente peut être établie dans les espaces de faibles ainsi 
de fortes valeurs de densité, et suivant l'aspect relatif du gradient réduit, on peut envisager deux 
situations possibles : 

e _S(p) — « pour lequel p est prédominante par rapport à Vp (p— 0), cette situation correspond 
à des régions de faible densité. 

e _s(p) — 0 dont Vp— 0 est prédominant par rapport à p , c'est à dire que même pour des 
valeurs élevées de p, le gradient réduit reste relativement très faible. Les points de l'espace 
concernés correspondant à cette situation, sont généralement les points critiques de liaisons 
(BCP). 

Un moyen classique qui permet d'identifier la présence ou l'absence d'interaction ente électrons 

partagés est le Laplacien |V?p|, qui caractérise la déplétion ou la concentration de densité de 

charge dans un espace donné du système, mais l'inconvénient de ce moyen mathématique, est 
au niveau de la caractérisation du type d'interaction en jeu. Pour cela, la NCI présente l'avantage 
de définir le type de contribution ou d'interaction non covalente à partir de l'exploitation du 

signe de la deuxième valeur propre À2 de la matrice Hessienne [75, 76]. 

La visualisation des espaces d'interactions NCI est non seulement liée aux limites du 
gradient réduit (p), mais aussi pour des valeurs intermédiaires, et suivant la valeur du RDG, on 
peut caractériser tous les contributions possibles dans le système, à savoir les interactions non- 
covalentes faibles et fortes ainsi que des liaisons. 

Selon les valeurs du (p) ainsi que le signe de À, on peut distinguer trois types 


d'interactions NCI [72] représentés généralement par des 1sosurfaces de différentes couleurs. 
i) Interaction forte stabilisante (type D) : 


Cette interaction concerne les régions d'espaces à valeurs élevées de p. La force de 
contribution dans ce type d'interactions est purement liée au gradient Vp. De plus, la nature de 
ce type de force peut être identifié par le signe négatif de À ((42) < 0). Cette interaction à effet 
attractif est généralement engendrée entre les groupements donneurs et accepteurs à l'intérieur 
du système, elle est bi-centrique, et hautement localisée. Ce type de contribution est appelée 


aussi liaison hydrogène. 
ii) Interaction forte déstabilisante (type ID) : 


Appelée aussi interaction de répulsion stérique, son effet purement déstabilisant est 
complètement opposé à celle de type I. L'identification de cette force est définie par (p signe 


(2) > 0). Il s'agit d’une contribution répulsive entre les pairs d'atomes ainsi que les contraintes 
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externes, autrement dit, la répulsion stérique est due à une somme de contributions d’atomes 


proches et multicentriques. 
iii) Interactions faibles ou de Van Der Waals (type II) : 


Une interaction de Van Der Waals est une contribution délocalisée située dans une zone 
d'espace de faible densité de charge p. De même, le gradient réduit RDG pour ce type de force 
peut avoir toutes les valeurs possibles dans l'espace réel, et ceci suivant la plage de variation de 


VPp. 


Contrairement au principe d'identification d'interactions de type I et IL, qui sont purement lié à 
psigne (2), celles de Van Der Waals ne sont identifiées que par la faible valeur de p. De plus, 
elles sont considérées comme facteur important qui contribue d'une façon significative à 


l'énergie totale du système. 
11.10.6.3 Visualisation NCI 


Le moyen de visualisation de l'ensemble d'interactions de notre composé SiC simulé 
dans ce travail est effectué par le code Critic [75], qui permet l'exploitation des valeurs propre 
de la matrice Hessienne ainsi que la densité de charge électronique p munie de son gradient Vp. 


Il permet aussi d'attribuer une couleur à chaque type d'interaction (figure IIL.4) : 


* Interactions stabilisantes (type I) en bleu 
* Interactions déstabilisantes (type Il) en rouge 


* Interactions faibles (type ID) en vert 


Plot gradient isosurfaces in real space and color by sign(2:,)2 


increasing ? =( increasing » 


d= 0 


es < > 


Strong, attractive Weak interactions: | Strong, repulsive 


interactions: interactions: 
van der Waals, 


H-bonding including dispersion steric crowding 


Figure IIL.4 Echelle de couleur RVB pour la représentation des interactions NCI suivant un 


ordre croissant de densité de charge [291. 


93 


Chapitre II Méthodes théoriques de la modélisation moléculaire 


Les iso-surfaces d’une représentation NCI en trois dimensions de l'ensemble des 
interactions NCL, présentent la signature de l'interaction ou de liaison en jeu dans le système 
[77]. Cette visualisation permet dans un cadre général d’identifier le type de liaison suivant le 
signe de À , sans avoir la valeur du RDG. Pour cela une version plus pratique, qui assure une 
représentation détaillée de toutes les contributions électrostatiques, elle montre les variations 
du RDG en fonction de p signe (2). Son avantage est d’identifier l’interaction dans chaque 


point de l’espace à travers son indice de courbure de densité de charge totale A2 


IIL.10.7 L'analyse de décomposition de l'énergie (EDA) 


L'analyse de décomposition en énergie (EDA) est une méthode quantitative puissante 
permettant de décrire les différents types d'interactions intervenant dans le système [68]. Elle a 
été formulée par Kitaura et Morokuma puis par Ziegler et Rauk [77,78]. La méthode consiste à 
décortiquer les interactions qui constituent une liaison chimique entre des fragments d'une 
molécule. En EDA, l’énergie de liaison totale (AEbonding) est décomposée en quatre termes, 
énergie électrostatique (AEelectrostatic), énergie cinétique (AEkinetic), énergie de Coulomb 


(AEcoulomb) et énergie de dispersion (AEdispersion) 


AE bonding = AE eict + AEkinetic + AEcoulomb + ÂAEdispersion (IL.20) 


AËEelct : correspond à l'énergie d'interaction entre les différentes distributions de charge non 
perturbées des fragments. 

AËEkinetic : c'est le terme d'énergie cinétique. 

ÂEcoulomb : contient les contributions énergétiques stériques et orbitalaire. 


AE dispersion : est le terme de dispersion provenant de la correction de Grimme 


IIL.11 Logiciels utilisés 
IL.11.1 Gaussian09 


Gaussian09 (version 2009) [791 est le avant dernier programme de structures 
électroniques de de la série Gaussian. Cette dernière est un logiciel de chimie numérique, créé 
en 1970 par John Pople et collaborateurs [80], Il peut être utilisé par les chimistes, les ingénieurs 
chimistes, les biochimistes, les physiciens et autres pour étudier les molécules et les réactions 
dans le cadre d’une large gamme de conditions, y compris les espèces et les composés qui sont 
difficiles ou impossibles à observer expérimentalement comme la courte durée de vie des 


structures intermédiaires et de transition. À partir des lois fondamentales de la mécanique 
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quantique, Gaussian prédit les énergies, les structures moléculaires, et les fréquences de 
vibration des systèmes moléculaires, ainsi que de nombreuses propriétés moléculaires issus de 


ces types de base de calcul. 


IIL.11.2 Gauss View 


GaussView [81] est une interface graphique complète pour rendre l'utilisation de Gaussian 
plus intuitive et visualiser les résultats. Avec GaussView on peut construire les systèmes 
moléculaires qui nous concernent rapidement et efficacement, en utilisant la fonction de 
construction des molécules. On pout également l'utilisez pour mettre en place et lancer vos 
calculs Gaussian, et pour visualiser vos divers résultats. 

GaussView comprend un excellent constructeur de molécules, permettant une 
construction rapide, même pour des grosses molécules : 
- Construction de molécules par atomes, cycle, groupe et acide aminé. 
- Importation de molécules d'autres sources, en les ouvrants tout simplement. 
- Vous pouvez également ajouter automatiquement des hydrogènes aux structures 
provenant de fichiers PDB, avec une excellente fiabilité. 


- Rotation en 3 dimensions même pour de très grosses molécules. 


IL.11.3 AIMAII 


AIMAII [82] est un progiciel permettant d'effectuer des analyses quantitatives et 
visuelles des systèmes moléculaires par la méthode QTAIM (théorie quantique des atomes dans 
les molécules), à partir de données de fonction d'onde moléculaire.Deux des composants 
d'AIMAII (AIMExt et AIMInt) sont des dérivés très fortement modifiés et étendus de deux 
programmes (Extreme et ProaimV) du progiciel AIMPAC développé et mis à jour par des 
membres du groupe de recherche de Richard F.W. Bader 
L'intérêt soutenu pour QTAIM et sa mise en œuvre pratique ont finalement abouti au 


développement de AIMAII, qui représente un effort de temps perdu depuis plusieurs années. 


IIL.11.4 Hyperchem 7.0 


Hyperchem 7.0 est un logiciel de modélisation moléculaire développé par Autodesk, 
INC, et distribuer hypercube, INC (Ontario, Canada) [83]. Hyperchem peut faire des calculs de 
la mécanique moléculaire et de la dynamique moléculaire, il offre aussi la possibilité de faire 


des calculs semi-empiriques et même quantique simple. 
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IIL.11.5 Logiciel ChemDraw Ultra 8.0 


ChemDraw Ultra version 8.0.3 est un logiciel développé en 2003 par Cambridge Soft. 
Ce logiciel permet à l’utilisateur de créer des modèles 3D, de visualiser, d’analyser et 
déterminer les propriétés des molécules. C’est un logiciel exhaustif qui offre à l’utilisateur un 


outil exceptionnel de modélisation et une interface intuitive et facile à utiliser [84]. 
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IV.1 Introduction 


Dans ce chapitre nous allons étudier l’importance des différentes interactions non- 
covalentes des complexes de transfert de charge (TC) entre le B-cyclodextrine (B-CD) comme 
donneur d'électrons et la 2,3-dichloro-5,6-dicyano-1,4-benzoquinone (DDQ), et le 
tétracyanoéthylène (TCNE), comme accepteurs d’électrons x (figure IV.1). Ces deux 
complexes ont été étudiés à l’état solide avec une stæœchiométrie 1:1 et synthétisés puis 
caractérisés par les techniques spectroscopiques usuelles telles que IR-FT, et la RMN du 'H. 
L'étude par spectroscopie IR-FT et RMN du 'H la masse aussi est intéressante a confirmé que 
la réaction a précédé la formation d'un complexe de transfert de charge (CTC). Les distinctions 
remarquées dans les bandes vibrationnelles observées du complexe et de leurs constituants pour 


les groupes fonctionnels, c'est-à-dire C = O, O-H, CN, C-H et CCI participant à la formation 


de liaisons H, confirment la formation du complexe actuel. Les complexes CT ont été formés 
par transfert de proton de f-CD vers DDQ et TCNE montrant des liaisons N°--H-O° et C=0*-- 
H. Comprendre l'interaction entre donneur et accepteur capable d'interpréter l'interaction 


donneur-accepteur et le mécanisme d'action entre ces composés. 


IV.2 Produits chimiques et Instruments 


Tous les produits chimiques ont été achetés auprès de Sigma-Aldrich; leurs puretés et 
leurs numéros CAS sont répertoriés dans le tableau (1V.1). L'eau a été distillée deux fois et 


déionisée avec une conductivité de 5 LS/cm à 25°C. 


Tableau IV.1 Puretés et numéros CAS des produits chimiques utilisés dans ce travail. 


produits chimiques CAS Pureté 
B-cyclodextrine 7585-39-9 >97% 
TCNE 670-54-2 >98% 

DDQ 84-58-2 070 


Les IR-FT des réactifs et des complexes formés sous forme de pastilles de KBr sont 
enregistrés en mode transmission dans le domaine des nombres d’onde de 400 à 4000 cm! en 
utilisant un spectromètre Perkin Elmer IR-FT. 

Le spectre RMN du l'H a été recueilli sur un spectromètre Bruker DRX- 250 à 400 MHz. 


Les mesures ont été effectuées à température ambiante en utilisant le DMSO-d6 


(diméthylsulfoxyde, de) comme solvant et le TMS (tétraméthylsilane) comme référence interne. 
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Figure IV.1 Structures chimiques de B-cyclodextrine, DDQ, et TCNE. 


IV.3 Préparation des complexes solides 


Les deux complexes solides DDQ et TCNE avec B-CD: [(B-CD)(DDQ)] et [(B- 
CD)(TCNE)] ont été préparés a un rapport molaire de 1:1. Pour chaque complexe, 1 mmole de 
B-CD dans l'eau (10 ml) a été ajoutée goutte à goutte à 1 mmole de chaque accepteur dans le 
méthanol (10 ml). 

Les mélanges ont été agités pendant une nuit et laissés s'évaporer lentement à 
température ambiante. Un changement de couleur a été observé lors du mélange de la solution 
de f-CD avec l'un quelconque des accepteurs. Ces nouvelles couleurs observées sont rouge- 
brique foncé pour [(B-CD)(DDQ)] et jaune pâle pour (B-CD)(TCNE)]. Ces changements de 
couleurs représentent une preuve solide de l'interaction supramoléculaire entre le donneur (f- 


CD) et chacun des accepteurs-r (DDQ et TCNE). 
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IV.4 Résultat et discussion 
IV.4.1 Infrarouge 


La preuve de l’interaction du donneur avec chacun des accepteurs est la présence des 
principales bandes dans le spectre du complexe résultant. La comparaison des spectres du 
donneur et de chacun des accepteurs seuls avec le spectre de leur complexe correspondant 
permet de relever des modifications en fréquence et/ou d’intensité des bandes suite à la 
formation des complexes. 

Les bandes de spectres IR-FT de la G-cyclodextrine donneuse libre, des accepteurs libres 
(DDQ, TCNE) et de leurs produits CT solides correspondants sont présentées sur la figure IV.2 
et la figure IV.3, tandis que les attributions de leurs bandes spectrales IR-FT caractéristiques 
sont rapportées dans le tableau IV.2. La formation des complexes CT lors de la réaction de la 
B-cyclodextrine avec TCNE et DDQ est fortement étayée par l'observation des principales 
bandes infrarouges du donneur et des accepteurs dans les spectres du produit. 

Le spectre IR-FT des bandes distinguées dans le donneur libre (B-CD) montre la 
vibration suivante : Une large bande correspond à l'étirement vibrationnel v (0H) observé à 
3283,58 cm!, et v (C-H) vibrations à 2925,37 cm !. Le spectre IR de la B-cyclodextrine est 


caractérisé par d'autres modes d'élongations à 1026,86 cm'!, attribués en v (C — O). 


IV.4.1.1 Complexe [(B-CD) (DDQ)] 


Le spectre IR-FT du complexe [(B-CD) (DDQ)] (figure IV.2) a indiqué que la bande 
résultant de la vibration v (C=N) de l'accepteur libre DDQ change de fréquence et diminue 
d'intensité après la complexation TC. À titre d'exemple représentatif, la vibration v (C=N) du 
DDAQ libre est observée à 2229,42 cm'!, tandis que dans son complexe, v (C=N) apparaît à la 
valeur de nombre d'onde inférieure (2217,48 cm !). Il est clair que v (C=N) de DDQ diminue 
au cours de la complexation. La bande caractéristique du groupe v (O-H) observée à 3283,58 
cm! dans la donneuse libre (B-CD) est déplacée à 3273,34 cm’! dans le complexe, et son 
intensité diminuée. Ces observations confirment clairement que le groupe OH du donneur f- 
CD et le groupe C=N de l'accepteur DDQ ont participé au processus de complexation. 

Le groupe cyano (C=N) est un groupe électro-attracteur qui existe dans le DDQ dans un 
système de liaisons conjuguées. La densité électronique rn*—-CN semble augmenter et accepter 
plus facilement un proton du donneur à cause du processus de retrait d'électrons et du système 
d'électrons conjugués. Ainsi, le mode d'interaction entre f-CD et l'accepteur DDQ se produit 


également par la migration d'un ion H* (Proton phénolique -OH ) vers l'un des groupes cyano 
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de l'accepteur DDQ pour former un ion positif (-C=N* H), qui s'associe à l'anion (—O°) pour 


former des paires d'ions [1,2]. 


1415.77 


1639.23 
1168.44 — 


1644.34 


| 1410.66 
3273.34 


B-CD/DDQ 


Nombre d'onde(cm-1) 


Figure IV.2 IR-FT du donneur libre (B-CD), DDQ libre, et le complexe de transfert de charge 
correspondant [(B-CD)(DDQ)|]. 


Le pic vibrationnel pour la formation du complexe CT observé à 3011,82 cm! devrait 
montrer une vibration d'étirement v (O°...H-—-N*). La présence de cette liaison (O°...H-N*) 
affirme que le nouveau complexe a été formé entre B-CD et DDQ. Cette liaison formée par le 
transfert d'électron de l'azote de DDQ à l'hydrogène de B-CD démontre que la liaison hydrogène 
est généralement faible. Les groupes de bandes attribuées aux vibrations v (C=O) et v (C=C) 
apparaissent respectivement à 1639,23 cm! et 1552,12 cm! dans le DDQ libre, avec un 
décalage vers les nombres d'onde supérieurs à 1644,34 cm! et 1582,94 cm'|, respectivement. 
La liaison (C=0*...H-0) a été formée en raison de l'interaction entre le donneur fB-CD et 


l'accepteur DDQ via une liaison hydrogène et un pic obtenu à 2966,54 cm‘. 
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Tableau IV.2 Fréquences infrarouges (em'!) et attributions provisoires pour la B-cyclodextrine, 


la DDQ, la TCNE et les complexes de transfert de charge formés. 


CTC 
B-CD DDQ TCNE . Assignements 
B-CD/DDQ : B-CD/TCNE 
3283.58 3273.34 : 3268.23 v (O-H) 
2925.37 2921.96 2920.25 v (C-H) 
1415.77 1410.66 1414.07 v (CH) of CH 
1026.86 1023.45 1021.74 v (C-0O) 
3011.82 3003.83 v (N*--H-0°) 
2966.54 v (C=0*--H) 
H- bonding 
2229.42 2217.48 v (CEN) 
1639.23 1644.34 v (C=0) 
1552.12 1582.94 v (C=C) 
1168.44 1151.38 v (CC) 
701.06 697.65 v (C-CD 
2255.01 : 2224.30 v (CEN) 
1577.82 _ 1589.76 v (C=C) 
1154.79 : 1163.32 v (CC) 


IV.4.1.1 Complexe [(B-CD)(TCNE)] 


Le spectre IR-FT du complexe [(B-CD)(TCNE)] (figure IV.3) montre que la vibration v 
(C=N) est observée à 2255,01 cm! dans le TCNE libre est décalée à 2224,30 cm! dans le 
complexe , et son intensité en diminué. La bande caractéristique du groupe v (0-H) observée à 
3283,58 cm! dans le donneur libre (f-CD) est décalée à 3268,23 cm! dans le complexe, et son 
intensité diminue. 

La liaison N°... H-0") a été formée en raison de l'interaction entre le donneur et 
l'accepteur via une liaison hydrogène et un pic obtenu à 3003,83 cm'!. Les bandes 
caractéristiques de C-C et C=C observées à 1154,79 cm! et 1577,82 cm! de TCNE libre, 
respectivement, sont décalées à 1163,32 cm ! et 1589,76 cm'!. 

Généralement, la présence de divers groupes donneurs de liaisons hydrogène de la f- 
cyclodextrine, comme l'alcool (-OH) et (-CH) facilite la liaison hydrogène et les interactions 
CT avec les molécules acceptrices (TCNE et DDQ). D'autre part, la présence d'électrons de 
paires isolées dans les atomes d'oxygène de DDQ et du groupe de retrait C=N dans TCNE fait 
de ces deux régions des régions acceptant les électrons et contribue par conséquent à la 


formation de CTC. 
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Figure IV.3 IR-FT du donneur libre (B-CD), TCNE libre, et le complexe de transfert de 
charge correspondant [(B-CD) (TCNE)|]. 


IV.4.2 RMN du ‘'H 


L'étude spectroscopique RMN'H est l'outil le plus important pour établir les complexes 
d'inclusion [3] et complexes de transfert de charge. Les complexes d'inclusion avec f-CD ont 
été largement étudiés et parmi de nombreuses méthodes spectroscopiques, celle de la RMN a 
été souvent utilisée [4-8]. 

Dans la structure de CD, on sait que les hydrogènes H-3 et H-S sont situés à l'intérieur 
de la cavité conique, en particulier, les protons H-3 sont placés près du bord de l'enrouleur 
tandis que les protons H-S sont placés près du bord le plus étroit et les autres protons H-1, H-2 
et H-4 sont situés à l'extérieur des cyclodextrines [9]. Si une molécule invitée est incorporée 
dans la cavité B-CD (H-3 et H-S), les constantes d'écran des protons B-CD à l'intérieur de la 


cavité devraient être sensibles au changement d'environnement, mais celle des protons 
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extérieurs (H -1, H-2, H-4 et H-6) ne devraient pas. Cela devrait entraîner des changements de 
déplacement chimique des protons à l'intérieur. 

Le tableau VI.3 montre les déplacements chimiques observés pour H-1, H-2, H-3, H-4, 
H-S et H-6 pour la B-cyclodextrine pure [10] et ses complexes de transfert de charge avec DDQ 
et TCNE. 


Tableau 1V.3 Les déplacements chimiques (8 ppm) des protons CH de la B-cyclodextrine avant 


et après la complexation avec DDQ et TCNE en solution D20. 


Compound H-1 H-2 H-3 H-4 H-5 H-6 H', 4" 
B-CD 5.047 3.606 3.801 3.459 3.781 3.624 / 
B-CD /DDQ 5.042 3.556 3.932 3.325 3.817 3.839 3.323 
5.364 
AG -0.005  -0.050 +0.131 -0.134  +0.036  +0.215 / 
B-CD/TCNE 5.043 3.591 3.944 3.327 3.823 3.851 3.323 
AG -0.004 -0.015 +0.143  -0.132  +0.042  +0.227 / 


Les changements de déplacement chimique les plus importants ont été observés pour H- 
3 (0,131 ppm) suivi de H-5 (0,036 ppm) dans le DDQ et pour H-3 (0,143 ppm) suivi de H-5 
(0,042 ppm) dans le TCNE. On observe que les valeurs de déplacement chimique pour les 
protons de la B-cyclodextrine (H-3 et H-5) après formation de CT accru, ce phénomène prouve 


que la formation d'inclusion et le complexe de transfert de charge. 


IV.4.2.1 Complexe [(B-CD) (DDQ)] 


Le spectre et les déplacements chimiques en RMN du ‘H pour le complexe de transfert 
de charge formé entre le B-CD et le DDQ sont représentés sur la figure IV.4. Le spectre 'H- 
RMN du complexe [(B-CD)(DDQ)] montre un nouveau pic observé à 3.323 ppm attribué au 
proton (*NH). La migration de l'ion H* du groupe (-OH) dans le B-CD à l'un des deux groupes 
cyano dans l'accepteur DDQ a entraîné la formation d'un ion positif (-C=N*H), qui est associé 
à l'anion O. 

Ces résultats confirment que le (-OH) de f-CD et (-C=N) de DDQ sont principalement 
impliqués dans la formation du complexe de transfert de charge entre B-CD et DDQ. Dans le 
même temps, on observe un nouveau signal à 5,364 ppm qui disparaît par deutération et son 
intégration avec le proton hydroxyle OH équivaut à un proton. Ceci est attribué à (C = O*...H) 
à la suite d'un transfert partiel de protons du groupe (-OH) du $-CD vers l'un des deux groupes 
carbonyle de l'accepteur DDQ. 

À partir des résultats obtenus, un mécanisme a été proposé pour la réaction du B-CD 


avec la DDQ, et celui-ci est montré dans la figure IV.5. 
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H1 
[(B-CD)(DDQ)] 


H" (5.364) ! 


pbigsatt is sn | RD T Ter TT Pr TT T TT 1 
75 7.0 6.5 eo s5 5.0 45 40 35 3.0 25 2.0 15 ppm 


Figure IV.4 Spectre et déplacements chimiques en RMN du 'H pour le complexe [(f- 
CD)(DDQ)|. 


Complexe de transfert de charge de f-CD:DDQ 


Figure IV.5 Mécanisme proposé pour la réaction du fB-CD avec la DDQ. 


IV.4.2.2 Complexe [(B-CD)(TCNE)] 


Le spectre et les déplacements chimiques en RMN du ‘H pour le complexe de transfert 


de charge formé entre le f-CD et le TCNE sont représentés sur la figure IV.6. Il est bien évident 
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dans le complexe [(B-CD)(TCNE)] que le nouveau pic observé à 3,323 ppm dans le spectre de 
ce complexe est attribué au proton de (*NH). Ce résultat est également confirmé par la 
formation du complexe CT entre B-CD et TCNE comme il est présenté sur la figure IV.7 


(mécanisme proposé). 


3.624 
7 ses H6 
re 
5.047 
o 
617 H5 
3.944 Ho OH\.3:616  H1 | H3 Fe 
| H4 
-CD)(TCNE | 
ICE x )] AE | H'(3.323) 
9 8 7 - - 4 3 2 1 con 


Figure IV.6 Spectre et déplacements chimiques en RMN du 'H pour le complexe [(f- 
CD) (TCNE)|. 


Complexe de transfert de charge de B-CD:TCNE 


Figure IV.7 Mécanisme proposé pour la réaction du f-CD avec la TCNE. 
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V.1 Introduction 


L’étude théorique des complexes moléculaires joue un rôle important dans la 
description de la stéréochimie de ces composés. Les résultats théoriques sont en fait 
considérés comme un complément des résultats expérimentaux. Les calculs de la 
modélisation moléculaire, permettent de déterminer à la fois la géométrie optimale et 
les paramètres essentiels comme les longueurs de liaisons, les angles de valence, les 
angles de torsion, les charges Mulliken et les coordonnées des atomes. La modélisation 
moléculaire permet aussi de calculer l’énergie stérique: paramètre caractéristique de la stabilité 
des complexes et évaluer ses propriétés physicochimiques de ces complexes. 

Dans ce chapitre nous allons discuter l’étude de la théorie fonctionnelle de la densité 
(DET) pour les deux complexes [(B-CD)(DDQ)] et [(B-CDI(TCNE)], en utilisant la 
fonctionnelle B3LYP avec l’ensemble de base 6-31G (d,p). Après cela, les énergies des 
orbitales moléculaires de frontières ont été estimés. L’analyse des bandes orbitales naturelles 
(NBO) et la théorie quantique de l’atome dans une molécule (QTAIM) ont été appliquées 
comme des approches puissantes pour l’évaluation des interactions intermoléculaires entre les 


molécules DDQ, TCNE et B-CD. 
V.2 Méthodologie de calcul 


La méthodologie que nous adapterons, consistera à appliquer les méthodes de la 
chimie quantique aux complexes moléculaires, dans le but' d’en déduire avec une très grande 
précision les différentes propriétés chimiques et physiques. 

Dans notre travail, tous les calculs sont effectués à l'aide du programme Gaussian 09 [1] et 
Gauss-View [2]. Pour aborder l’étude théorique des complexes de transfert de charge de 
stæchiométrie 1:1, nous avons commencé nos calculs par une minimisation de l'énergie, en 
utilisant la méthode semi empirique PM6, ensuite, afin d'affiner les résultats, nous avons 
réalisé d'autres optimisations énergétiques par l’utilisation de la méthode DFT/B3LYP. Les 
calculs QTAIM ont été effectués sue ces structures optimisés à l’aide du logiciel AIM2000 [3] 
et du logiciel AIMALL [4]. Cependant, l’analyse de la NBO a été réalisée au moyen du 
programme NBO 3.1 qui est implémenté dans le logiciel Gaussian 09 [1]. 

La figure V.1 illustre les structures des molécules libres optimisées par la méthode B3LYP/6- 


31G (d,p) avec numérotation des atomes. 
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Figure V.1 Structures moléculaires du donneur (c): B-CD et accepteurs (a): DDQ et (b): 
TCNE optimisés (dans le vide) par la méthode B3LYP/6-31G (d,p). 


L'énergie de stabilisation (AE) et l'énergie d'interaction (AEix) lors de la 
complexation entre le donneur (fB-CD) et les accepteurs (DDQ et TCNE) peuvent être calculées 


pour les structures d'énergie minimale selon les équations suivantes : 
ÂEstb= Ecomplex-(ED_opt+EA_opt) (V. 1) 
ÂEint = Ecomplex-( Ep sp+ EA p) (V2) 


Ebpet Ea représente l’énergie du donneur et celle de l’accepteur respectivement optimisé à l’état 
libre. 

Les résultats correspondants sont indiqués dans le tableau V.1. Les énergies de 
complexation ont été trouvées égales à -2963679,103 Kcal/mol pour le complexe [(f- 


CD)(TCNE)] et à -3614768,103 Kcal/mol pour [(B-CD)(DDQ)|]. 


D 
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Tableau V.1 Valeurs d’énergie d’interaction (AE) des complexes [(B-CD)(DDQ)] et [(B- 
CD)(TCNE)] calculées par la méthode B3LYP/6-31G (d,p). 


Energies (kcal/mol) [(B-CD)(DDQ)] [(B-CD)(TCNE)|] 
AE (kcal/mol) 3.800 4.000 
ÂEïint (kcal/mol) 3.600 3.800 


D’après les résultats du tableau V.1, nous avons constaté que les valeurs de l’énergie 
d’interaction (AE) du complexe [(fB-CD) (TCNE)] sont plus élevées que le complexe [(B-CD) 
(DDQ)]. 


V.3 Etude des complexes de Transfert de charge: CTC1: [(B-CD) (DDQ)] et CTC2: [(B- 
CD) (TCNE)] 


V.3.1 Géométrie et optimisation 


En plus de l’énergie stérique, les nombreux résultats de calculs d’optimisation nous ont 
fournis les charges Mulliken et les distances interatomiques (élongation de 
haison). Donc nous avons effectué une comparaison des paramètres structuraux (les longueurs 
de liaison et les angles) des molécules libres (donneur: B-CD et accepteurs: DDQ et TCNE) et 
leurs complexes optimisés par la méthode B3LYP/6-31G (d,p).(Voir tableaux V.2 et V.3). 

À partir de la figure V.1, l'interaction entre les deux réactifs (donneur et accepteur) s'est 
manifestée par plusieurs variations dans les longueurs de liaison, les angles de liaison et angles 
dièdres dans les deux molécules. En configuration tête-à-tête, 1l était évident que les 
modifications structurelles ont été réparties sur un grand nombre de liaisons dans les deux 
molécules. 

Dans la molécule DDQ, les longueurs des liaisons relativement sensibles C148=0155, 
C149-C1161, C151=0154 et C150-C1160 ont été embardés de -0,0004 à -0,0054 A°. Pour 
C153-C156,C148-C149,C150-C151,C148-C153 et C151-C152 les changements situés dans 
la gamme 0,0014 à 0,0063 ont été également remarqués. 

Dans la molécule TCNE, les longueurs des liaisons C152=N153, C148-C150 
CIS4=NI15S, C149-C156, CIS6=N157, C149-C154, C148-C152 et C148=C149 ont été 
embardés de 0,0011 à 0,0041 A°. Dans le plan opposé, plusieurs modifications dans presque 
toutes les liaisons constitutives ont été subies dans la molécule B-CD. Dans [(B-CD) (DDQ)], 
les principaux écarts de longueur de liaison ont été observés en 057-H79 et 071-H84 avec des 
écarts de 0,0144 et 0,0156 À, respectivement. Les déviations des liaisons adjacentes 073—H82, 
057-H79 et 076-H80 étaient de 0,0146 À, 0,0149 À et 0,0154 À, en [(B-CD) (TCNE)]. 
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Pour les angles et en CTC1, les déviations les plus importants de 3,5256 A° et -93,475 
A° ont été constatés respectivement pour C153-C156=N157 et C148-C153=C152. De même, 
dans CTC2, les déviations des angles de liaison C149-CIS6=N157, C156-C149=C148 et 
C148-CI150=N151 étaient de 1,0864, 1,3461 et 3,5506 A°, respectivement. Il a également été 


noté que les angles dièdres ont subi de grands changements. 


Par exemple, dans CTCI, les écarts pour C148-C153-CIS6=N157 et C148-— 
C153=C152-C158 étaient très remarquable et égale 128.877 et 354.886 A°. Pour le CTC2, les 
angles dièdres de C149=C148-C152=N153, C149-C148-C150=N151 et C154-C149=C148- 
C150 sont -15,756, 324,377 et 357,359 A° respectivement. 


Tableau V.2 Paramètres géométriques (longueurs de liaisons, angles valentiels et angles 
dièdres) pour les molécules libres: accepteur (DDQ), donneur (fB-CD) et leurs complexes 


(optimisé dans le vide) calculés avec la méthode B3LYP/6-31G (d,p). 


B-CD DDQ CTCI AÔ 
Longueurs de liaison (A') 
O71-H84 (0.9831) 0.9675 0.0156 
O60-H83 (0.9833) 0.9766 0.0067 
O61-H81 (0.9829) 0.9676 0.0153 
O57-H79 (0.9832) 0.9688 0.0144 
073-H82 (0.9830) 0.9674 0.0156 
O076-H80 (0.9831) 0.9838 -0.0007 
077-H78 (0.9831) 0.9814 0.0017 
C35-H123 (1.0956) 1.0971 -0.0015 
C42-H100 (1.1007) 1.1004 0.0003 
C53-H114 (1.0956) 1.0952 0.0004 
C12-H134 (1.1006) 1.0988 0.0018 
Longueurs de liaison (A) 
C148-C149 (1.4956) 1.4939 0.0017 
C148-C153 (1.5034) 1.4981 0.0053 
C153-C156 (1.4268) 1.4254 0.0014 
C148=0155 (1.2136) 1.2140 -0.0004 
CIS6=N157 (1.1625) 1.1624 0.0001 
C149-C1161 (1.7184) 1.7197 -0.0013 
C150-C151 (1.4956) 1.4913 0.0043 
C151-C152 (1.5034) 1.4971 0.0063 
C152-C158 (1.4268) 1.4267 0.0001 
C158=N159 (1.1625) 1.1622 0.0003 
C150-C1160 (1.7184) 1.7250 -0.0066 
C151=0154 (1.2136) 1.2190 -0.0054 
Angles valentiels (°) 
C153-C156=N157 (178.7427) 175.2171  3.5256 
O155=C148-C153 (120.2425) 120.7284  -0.4859 
C1161-C149=C150 (123.0121) 122.9483  0.0638 
C148-C153=C152 (27.7541) 121.2300  -93.475 
Angles dièdres (°) 


C148-C153-C156=N157 (177.9728) 49.0958 128.877 
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O155=C148-C153=C152 (178.7619) 178.9474  -0.1855 
C1161-C149-C148-C153 (-178.8004) -179.4562  0.6558 
C148-C153=C152-C158 (180.0000) -174.8860 354.886 


Tableau V.3 Paramètres géométriques (longueurs de liaisons, angles valentiels et angles 
dièdres) pour les molécules libres: accepteur (TCNE), donneur (B-CD) et leurs complexes 


(optimisé dans le vide) calculés avec la méthode B3LYP/6-31G (d,p). 


B-CD TCNE CTC2 AÔ 
Longueurs de liaison (A) 
O71-H84 (0.9831) 0.9822 0.0011 
O60-H83 (0.9833) 0.9798 0.0035 
O61-H81 (0.9829) 0.9797 0.0032 
O57-H79 (0.9832) 0.9683 0.0149 
073-H82 (0.9830) 0.9684 0.0146 
O76-H380 (0.9831) 0.9677 0.0154 
077-H78 (0.9831) 0.9792 0.0039 
C35-H123 (1.0956) 1.0929 0.0027 
C42-H100 (1.1007) 1.0968 0.0039 
C28-H131 (1.0956) 1.0948 0.0008 
Longueurs de liaison (A') 
CISO=NI151 (1.1627) 1.1624 0.0003 
C148=C149 (1.3723) 1.3682 0.0041 
C152=N153 (1.1627) 1.1616 0.0011 
C148-C150 (1.4296) 1.4282 0.0014 
C154=N155 (1.1627) 1.1609 0.0018 
C148-C152 (1.4296) 1.4265 0.0031 
C156=N157 (1.1627) 1.1607 0.0020 
C149-C154 (1.4296) 1.4270 0.0026 
C149-C156 (1.4296) 1.4277 0.0019 
Angles valentiels (°) 
C149-C156=N157 (178.9316) 177.8452  1.0864 
C156-C149-C154 (116.9834) 118.5475  -1.5641 
C156-C149=C148 (121.5198) 120.1737  1.3461 
C148-CI1S0=N151 (178.9172) 178.9172  3.5506 
Angles dièdres (°) 
C149=C148-C152=N153 (-179.8478) -164.0913 -15.756 
C149=C148-CI150=N151 (-179.8211) 144.5562 324.377 
C154-C149=C148-C150 (179.9921 ) -177.3670 357.359 


Pour compléter cette prédiction structurelle, l'analyse de population Mulliken a été 
exploitée pour calculer la charge atomique afin de savoir quels atomes dans les réactifs ont été 
affectés sur le complexe de formation et dans quelle direction le processus de transfert de charge 
se produira. Dans les phénomènes complexes de transfert de charge, l'estimation des charges 
atomiques de Mulliken joue un rôle très impératif dans l'application des calculs de mécanique 
quantique aux systèmes moléculaires [5,6]. Ces charges atomiques de Mulliken sont présentées 


dans le tableau V.4. Comme le montre la Figure V.1, C1 et C2 (nommés respectivement C148 
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et C149 dans les complexes CT) dans le TCNE libre contiennent deux valeurs positives avec 


une base de 6-31G (d,p). La valeur correspondante est +0,19728 e. 


Tableau V.4 Distribution des charges atomiques de Mulliken de certains centres dans les 


complexes CTC1: [(B-CD)(DDQ)] et CTC2: [(B-CD)(TCNE)|]. 


Atom In CTCI1 In CTC2 Atom In CTCI1 In CTC2 
B-CD 

O1: -0.57002 -0.58833  -0.58276 CS53: 0.14756 0.15264 0.14700 
O2: -0.58355 -0.60864  -0.61586 057: -0.59970 -0.57474  -0.59545 
O3: -0.58328 -0.61274  -0.59593 O60: -0.59985 -0.60693  0.62577 
O4: -0.56997 -0.59049  -0.59464 O61: -0.59925 -0.55547  0.61009 
O5: -0.56982 -0.59949  -0.57994 O71: -0.59958 -0.57765  -0.62370 
C6: 0.15590 0.16762 0.16178 076: -0.59950 -0.62135  -0.57499 
C7: 0.15566 0.16080 0.15438 077: -0.59979 -0.61086  -0.62550 
O8: -0.58320 -0.62704  -0.60628 H79: 0.35807 0.34348 0.35030 
C9: 0.15594 0.158359 0.16130 H80: 0.35795 0.36385 0.33387 
O10: 0.58320 -0.61864  -0.59439 H81: 0.35737 0.31589 0.36118 
C11: 0.15523 0.17049 0.15769 H82: 0.35772 0.34862 0.33696 
C12: 0.15519 0.16197 0.17148 H83: 0.35809 0.36235 0.36999 
018: -0.52484 -0.53521  -0.54180 H84: 0.35773 0.33660 0.37087 
019: -0.56996 -0.59817  -0.58914 H100: 0.10004 0.11410 0.14451 
C20: 0.15330 0.15148 0.15624 H114: 0.12458 0.15728 0.13612 
C28: 0.14764 0.15003 0.15999 H120: 0.10807 0.11225 0.11092 
C29: 0.15532 0.16076 0.15535 H121: 0.09995 0.11189 0.11604 
C30: 0.36030 0.37344 0.360833 H122: 0.10027 0.106683 0.10221 
C35: 0.14717 0.12530 0.12481 H123: 0.12433 0.10666 0.16831 
C40: 0.15294 0.16101 0.13489 H130: 0.10695 0.11199 0.10763 
C41: 0.02404 0.03147 0.01987 H131: 0.12420 0.14735 0.13715 
C42: 0.15495 0.16298 0.15327 H132: 0.12791 0.12241 0.14286 
C50: 0.36019 0.36759 0.35333 H133: 0.11266 0.13150 0.12333 
C51: 0.02404 0.05274 0.01920 H134: 0.10010 0.11749 0.11061 
DDQ 

C148: 0.46637 0.47337 O155: -0.38247 -0.44853 

C149: -0.15244 -0.12062 C156: 0.29467 0.32784 

C150: -0.15237 0.12507 N157: -0.42986 -0.50695 

C151: 0.46629 0.51369 C158: 0.29466 0.299907 

C152: 0.08457 0.20164 N159: -0.42990 -0.45136 

C153: 0.08450 0.11236 CI160: 0.11920 0.09431 

0154: -0.38243 -0.48374 CI161: 0.11918 0.08147 

TCNE 

C148: 0.19728 0.24401 N153: -0.41851 -0.47098 
C149: 0.19728 0.25476 C154: 0.31987 0.39097 
C150: 0.31987 0.36363 N155: -0.41851 -0.51204 
N151:-0.41851 -0.47808 C156: 0.319387 0.41144 
C152: 0.31987 0.35004 N157: -0.41851 -0.55140 


La figure V.2 illustre les liaisons hydrogène intermoléculaires des complexes [(f- 
CD)(DDQ)] et [(B-CD)(TCNE)]. 
Comme le montre la figure V.2, deux liaisons hydrogène sont observées dans lastructure 


optimisée de complexe [(B-CD)(TCNE)], ceci indique la présence de transfert de charge entre 
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le donneur f-CD et l’accepteur TCNE. 

La première liaison hydrogène est formée par l'interaction de l’atome N155 du TCNE 
et l’atome H79 du donneur B-CD avec une distance de 1.79 À. 

La deuxième liaison hydrogène est située entre N157 (C=N) et H80 avec une distance 


de 2.06 À en utilisant la base 6-31G (d,p). 


Figure V.2 Liaisons hydrogène intermoléculaires dans les complexes (a): [(B-CD)(TCNE)] et 
(b): [(B-CD)(DDQ)1. 

Il a été constaté que les charges négatives du N155 et N157 dans l'accepteur (TCNE) 
diminuent de manière significative lors de la complexation avec la molécule donneuse. Ceci est 
associé à une diminution de la charge positive de H79 à partir de +0,35841 e jusqu'à +0,35030 
e et à une diminution de la charge positive de H80 à partir de +0,35827 e jusqu'à +0,33387 e 
dans une base de 6-31G (d,p). D'autre part, la distance de liaison C148....C149 du TCNE a été 
diminué de 1,3723 A° à 1,3682 A° à une base de 6-31G (d,p) définie à l'état libre et à l'état 
complexé (tableau V.3). De même, les charges positives de C148 et C149 de TCNE dans le 
complexe [(B-CD)(TCNE)] passent de 0,19728 e à 0,24401 e et 0,25476 e, suggérant qu'une 
quantité appréciable de charge électronique a été transférée du ligand au accepteur TCNE. 

La longueur des liaisons intramoléculaires C154....N155 et C156....N157 est 
également diminuée de 1,1627 A° (dans le TCNE libre) à 1,1609 A° et 1,1607 A° 
respectivement pour une base de 6-31G (d,p), ceci est due à l'effet de la liaison hydrogène. 
D'autres interactions ont été établies entre N151....H100, N153....H123 et N157....H82 avec 


des distances interatomiques : 2,61 A°, 2,42 A° et 2,21 A° respectivement. 
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Dans le cas du complexe [(B-CD) (DDQ)] formé dans le vide (figure V.2), trois liaisons 
hydrogène ont été établies : d’abord la liaison hydrogène pouvant se former entre N157 (C=N) 
et H84 de f-CD avec la distance interatomique de 2,16 A°. 

La deuxième liaison hydrogène est 0154....H79 de distance 2,10 A° et la troisième liaison 
hydrogène est formée entre O155 (C=O) et H131 de la B-CD avec une distance de 2.09 A°. 

D'autres interactions sont également trouvées entre C1160....H134 et N159....H83 où 
les distances de liaison sont respectivement de 2,80 A° et 3,83 A°. 

Il a été constaté que les charges négatives du O154 et O155 dans l'accepteur (DDQ) 
diminuent de manière significative lors de la complexation avec la molécule donneuse. Ceci est 
associé à une diminution de la charge positive de H79 à partir de +0,35841 e jusqu'à +0,34348 
e et à une augmentation de la charge positive de H131 à partir de +0,12429 e jusqu'à +0,14735 
e dans une base de 6-31G (d,p). 

A noter également que les longueurs de liaison de O57....H79, 073....H82 et 
076....H80 sont significativement modifiées après complexation (Tableaux V.2 et V.3) où les 
distances de liaison sont changées à partir de 0,9832 A°, 0,9830 A° et 0,9831 A° vers 0,9683 
A°, 0,9684 A° et 0,9677 A° dans [(B-CD)(TCNE)] et vers 0,9688 A°, 0,9674 A° et 0,9838 A° 
dans [(B-CD)(DDQ)] respectivement. 

Un autre point important est que les charges de Mulliken correspondant aux deux atomes 
de chlore dans la molécule accepteur DDQ sont entièrement modifiées, passant d'une valeur 
identique de 0,11920 e et 0,11918 e à 0,09431e et 0,081476 e respectivement pour C1160 et 
CH61. 

Des résultats similaires ont été trouvés dans l'analyse des longueurs de liaisons avant et 
après complexation du TCNE et du DDQ avec la molécule biologique f-CD (voir tableaux V.2 
et V.3). Une diminution de la longueur de la liaison O57....H79 de 0,9683 A° et 0,9688 A°, 
respectivement dans les complexes [(B-CD)(TCNE)] et [(B-CD)(DDQ)], et les charges 
atomiques de ces atomes (057 et H79) ont été observés comme étant diminués lors de la 
complexation. On peut noter qu'il existe une forte interaction de liaison hydrogène entre 
N155....H79 et O154....H79. Dans les complexes [(B-CD)(DDQ)] et [(B-CD)(TCNE)|}, les 
deux azotes (N157 de DDQ et N157 de TCNE) ont montré une augmentation des charges 
atomiques de Mulliken (-0,50695 e et -0,55140 e ). 

D'autre part, les charges des atomes H80 et H84 de B-CD étaient de 0,36385 e et 0,33660 
e dans CTCI, mais dans CTC2, elles étaient respectivement de 0,33387 e et 0,37087 e. Ces 


découvertes confirment clairement le CT entre les deux molécules accepteur-donneur. 


us 
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V.3.2 Analyse des orbitales moléculaires frontières (FMO) 


Afin de comprendre la stabilité globale et la réactivité chimique des deux complexes, 
nous avons étudié les orbitales moléculaires frontières. Les indices globaux de réactivité tels 
que le potentiel électronique (1), la dureté (n), adoucisseur chimique (S) et l'indice 


d'éléctrophilicité global (&) ont été calculés en utilisant les équations suivantes : 


u = (Exomo+ELumo) /2 (V3) 
n = (ELumo-Exomo) /2 (V4) 
S = 1/n (V5) 
@ = pu? /2n (V.6) 


L'écart d'énergie (HOMO-LUMO) est un indice de stabilité importante pouvant 
caractériser la réactivité chimique et la stabilité cinétique de la molécule [7]. 

Quand l’écart énergétique HOMO-LUMO est élevé, l’écoulement des électrons à 
l’état d’énergie plus élevée est difficile (stable), ce qui rend la molécule dure et moins 
réactif. D'autre part, un gap énergétique HOMO-LUMO plus faible, a pour 
conséquence un écoulement facile des électrons, ce qui rend la molécule molle, ceci 
est gouverné par le principe de HSAB (Hard Soft Acide and Bas) [8]. 

Les résultats des paramètres d'énergies HOMO, LUMO et de réactivité u, n et «© des 


complexes de transfert de charge sont résumés dans le tableau V.5. 


Tableau V.5 (LUMO+1, LUMO, HOMO, HOMO-1) énergies, potentiel électronique (11) (eV), 
dureté (n) (eV), adoucisseur chimique (S) (eV) et l'indice d'électrophilicité global (œ) (eV) du 
complexes de transfert de charge calculés par la méthode B3LYP/6-31G (d,p). 


Paramètres IB-CD)DDQ)] IB-CD\XTCNE)] 
ELumo+ (eV) -2.367 -1.717 
Ezumo (eV) 4,357 -4.,628 
Exomo (eV) -7.051 -6.922 
Exomo-1 (eV) -7.140 -6.998 
Exomo —ELumo (eV) -2.694 -2.294 
u (eV) 5.704 5775 
n (eV) 1.347 1.147 
S (eV) 0.742 0.871 
œ (eV) 12.077 14.538 


Les diagrammes orbitalaires moléculaires obtenus par la méthode DFT/B3LYP 


avec la base 6-31G (d,p) pour les différents composés sont comparés sur les deux 
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figures : figure V.3 pour le complexe [(B-CD)(DDQ)] et la figure V.4 pour le complexe 
[(B-CDTCNE)]. 


LUMO+1 


ELumo = -4.357 eV 


ExoMo = ELumo = -2.694 eV 


Figure V.3 Représentation des Orbitales Moléculaires (Oms) du complexe [(B-CD)(DDQ)] des 


Principales excitations avec leurs énergies 


=] 
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LUMO+1 


ExHoMo = -6.922 eV 


Figure V.4 Représentation des Orbitales Moléculaires (Oms) du complexe [(B-CD)(TCNE)] 


des Principales excitations avec leurs énergies 


Les figures V.3 et V.4 récapitulent l'écart entre l’orbitale moléculaire la plus haute 
occupée (HOMO-1) et l’orbitale moléculaire la plus basse vacante (LUMO+1). 

D’après les valeurs énergétiques des orbitales HOMO-I et LUMO+1, on remarque que 
le complexe [(B-CD)(TCNE)] donne la plus grande valeur d’énergie qui est égale -5.281 eV. 
En revanche le complexe [(B-CD)(DDQ)] donne la plus petite valeur de E qui est égale 
à: -4.773 eV. 

En plus, l'écart énergétique HOMO-LUMO des deux complexes [(B-CD)(DDQ)] et [(B- 
CD) (TCNE)|] sont respectivement de -2,694 eV et -2,294 eV. Nous pouvons conclure, à partir 


=] 
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des résultats rapportés dans le tableau V.5, que le complexe [(fB-CD) (DDQ)] avait le plus faible 
écart d'énergie (HOMO-LUMO), ce qui indique que le complexe [(B-CD)(DDQ)] est le plus 
stable et sa densité électronique a changé plus facilement que le complexe [(B-CD)(TCNE)|]. 

Toutes les valeurs de potentiel (1) des différents complexes de transfert de charge est 
négatif sous vide, indiqant que le processus d'encapsulation est spontané. Cependant, la stabilité 
d’une molécule est liée à sa dureté chimique (n), dans notre système elle est d’environ 1.347 
eV et 1.147 eV pour [(B-CD) (DDQ)] et [(B-CD) (TCNE)] respectivement, ce qui signifie que 
le complexe [(B-CD)(TCNE)] est moins stable (donc plus actif) que le complexe [(B-CD) 
(DDQ)] car il s’oppose à moins de transfert de charge et aux changements de sa densité 
électronique. 

L'indice global d'électrophilie (©) est un indice important d'électrophilie et de grandes 
valeurs sont des caractéristiques de la plupart des systèmes électrophiles [9]; ainsi le complexe 
[(B-CD) (TCNE)] est plus électrophile que [(B-CD) (DDQ)] pour la méthode B3LYP/6-31 G 
(d,p) dans le vide. 


V.3.3 Analyse NBO (Natural Bonde Orbital) 


Pour mieux comprendre l'interaction intermoléculaire des complexes [(B-CD)(DDQ)] et 
[(B-CD)(TCNE)], une analyse orbitale des liaisons naturelles [10] a été réalisée avec DFT/6- 
31G (d, p) méthode, qui fournit beaucoup d'informations sur la nature et la force des interactions 
intermoléculaires du point de vue des interactions orbitales locales et leur énergie de 
perturbation de second ordre (E°). Plus la valeur E° est grande, plus l'interaction entre les 
donneurs et les accepteurs d'électrons est intense [11]. 

La délocalisation de la densité électronique entre les orbitales occupées (liantes ou non- 
hantes) et les orbitales vacantes (anti-liantes) correspond à une interaction stabilisante donneur- 
accepteur [12]. Pour chaque donneur (i) et accepteur (j), l'énergie de stabilisation E 


[10] associée à la délocalisation 1 — j est estimée en fonction de l'équation. (V. 7): 
p0 = q =) (VD) 
Ej—Ei 
Où qi est l'occupation orbitale donneuse, si et ej sont les éléments diagonaux et Fij est l'élément 
non diagonal de la matrice Fock NBO. 
Le tableau V.6 montre les interactions donneur-accepteur et leurs énergies de 
stabilisation E° correspondant aux structures optimisées des complexes obtenues en utilisant 


la fonctionnelle B3LYP/6-31 G (d,p) dans le vide. 
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Résultats théoriques 


Tableau V.6 Les orbitales donneuses partielles d'électrons, les orbitales acceptrices d'électrons 


et les énergies correspondantes E(? (kcal/mol) pour les complexes [(B-CD)DDQ)] et [(B- 
CD)(TCNE)] par DET/6- 31G (d,p). 


Complexes Donneur NBO(i) Accepteur NBO(j) E® (kcal/mol) 

DDQ/ f-CD B-CD donneur et DDQ accepteur 
LP(1) O 60 BD*(1) C 158-N 159 0.46 
BD(1) C 12-H 134 BD*(1) C 150-CI 160 0.19 
BD(1) O 71-H 84 BD*(2) C 156-N 157 0.29 
LP(1) O 18 BD*(1) C 150-CI 160 0.10 
DD donneur et B-CD accepteur 
LP(1) N 157 BD*(1) O 71-H 84 2.70 
BD(1) C 156-N 157 BD*(1) O 71-H 84 0.07 
BD(2) C 156-N 157 BD*(1) O 71-H 84 1.63 
LP(1) O 155 BD*(1) C 53-H 114 3,17 
BD(2) C 148-0O 155 BD*(1) C 53-H 114 71 
LP(1) CI 160 BD*(1) C 12-H 134 0.62 
LP(1) O 154 BD*(1) O 57-H 79 2.57 
BD(2) C151-O 154 BD*(1) O 57-H 79 1.00 
BD(1) C 150-CI 160 BD*(1) C 12-H 134 0.09 
LP(1) CI 160 BD*(1) C 7-H 128 0.27 

TCNE/ f-CD B-CD donneur et TCNE accepteur 
LP(1)O61 BD*(2) C 156-N 157 0.12 
LP(2) O 61 BD*(3) C 156-N 157 0.58 
BD(1) O 61-H 81 BD*(2) C 156-N 157 0.10 
LP(1)057 BD*(2) C 154-N 155 0.13 
BD(1) O 57-H 79 BD*(2) C 154-N 155 0.77 
LP(1) O 76 BD*(1) C 156-N 157 0.11 
BD(1) O 76-H 80 BD*(1) C 156-N 157 0.21 
BD(1) C 28-H 131 BD*(3) C 152-N 153 0.08 
TCNE donneur et B-CD accepteur 
LP(1) N 153 BD*(1) C 35-H 123 2.93 
BD(2) C 152-N 153 BD*(1) C 35-H 123 0.08 
LP(1) N 155 BD*(1) O 57-H 79 9.44 
BD(1) C 154-N 155 BD*(1) O 57-H 79 0.24 
BD(2) C 154-N 155 BD*(1) O 57-H 79 1.40 
LP(1) N 157 BD*(1) O 73-H 82 2.53 
BD(3) C 156-N 157 BD*(1) O 73-H 82 1.20 
LP(1) N 157 BD*(1) O 76-H 80 7.29 
BD(3) C 156-N 157 BD*(1) O 76-H 80 0.59 
LP(1) N 151 BD*(1) C 42-H 100 1.42 


On peut voir à partir de ces résultats qu'un grand nombre d'interactions non-liantes-anti-liantes 
et liaison-anti-liantes se produisent entre la molécule donneuse (B-CD) et les accepteurs (DDQ 
et TCNE) (Figure V.2). 

Les énergies d'interaction de ces contacts sont comprises entre 0,07 kcal/mol-3,17 
kcal/mol pour [(B-CD) (DDQ)] et 0,08 kcal/mol-9,44 kcal/mol pour [(B-CD)(TCNE )] . Les 


interactions sont expliquées en détail comme suit. 
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Dans le complexe [(fB-CD) (DDQ)], une interaction importante est observée entre LP(1) 
0155 de DDQ et BD*(1) C53-H114 de B-CD et correspond à 3,17 kcal/mol. Une autre 
interaction importante est observée entre LP(1) N157 et BD*(1) O71-H84 correspond à 2,70 
kcal/mol. De plus, nous distinguons également une interaction non-liante-anti-liante entre les 
deux partenaires, cette interaction est établie entre BD(2) C156-N157 et l'anti-liaison vacant 
BD*(1) O71-H84 avec une valeur d'énergie de stabilisation (1,63 kcal/mol). Lorsque DDQ agit 
comme donneur, présence de liaison H intermoléculaire correspondant aux deux orbitales 
donneur-accepteur d'électrons : BD(2) C151-0154 vers BD*(1) O57-H79, LP(1) CI160 vers 
BD *(1) C12-H134, BD(2) C148-0155 vers BD*(1) C53-H114 et LP(1) O154 vers BD*(1) 
O57-H79 entraîne respectivement une stabilisation de 1,00 , 0,62, 1,71 et 2,57 kcal/mol. 

Dans le cas du complexe [(B-CD) (TCNE)], on note également qu'il existe une forte 
interaction de liaison hydrogène entre N155 et H79 qui peut être caractérisée par l'interaction 
orbitale LP (1) N155 vers BD*(1) O57-H79 et la valeur énergétique correspondante de cette 
interaction est de 9,44 kcal/mol. Une autre interaction importante est observée entre LP (1) 
N157 et BD*(1) 076-H80 correspond à 7,29 Kcal/mol. Deux autres interactions sont observées 
: la première entre LP (1) N151 et BD*(1) C42-H100 et la seconde entre LP (1) N153 et BD*(1) 
C35-H123 avec une énergie de stabilisation égale respectivement, à 1,42 et 2,93 Kcal /mol. 

D'après les résultats des calculs de NBO, nous pouvons conclure que le transfert de 
charge entre les orbitales occupées et vacantes de donneur/accepteur et les interactions 


hydrogènes intermoléculaires sont les contributions majeures de la stabilisation du complexe. 


V.3.4 Propriétés optiques non linéaires (NLO) 


Des investigations théoriques ont été menées afin de comprendre l'origine 
microscopique du comportement non linéaire de la molécule étudiée. L'étude implique le calcul 
des tenseurs de moment dipolaire (1), de polarisabilité (&) et de première hyperpolarisabilité 
(B) de la molécule. 

Le logiciel Gaussian 09 fournit six composants polarisables et chacun représente 
respectivement à Gxx, Oxy, Oyy, Oxz, Oyz, Ozz . Pour une hyperpolarisation, nous obtiendrons dix 
composants hyper polaires et que chacun représente respectivement à Bxxx, Bxyy Bxyy, Byyy Byzz 
Bxyz, Byz Bxzz, Byzz, Be. Pour l'hyperpolarisabilité, nous obtiendrons dix composants hyper 


polaires et chacun représente respectivement à Bxxx, Bxxy, Bxyy, Byyy, Bxxz, Bxyz, ByYz, Bxzr, 


Byzz, Bzzz. 
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En utilisant les composantes x, y et z, l'amplitude du moment dipolaire statique total, la 
polarisabilité isotrope et l'hyperpolarisabilité de premier ordre peuvent être calculées par les 


équations suivantes [13] : 


Ho (Ex? + pi? + p27) 7 (V.8) 

Œ = Oxx + Oyy + Ozz /3 (V.9) 

B = (BE + BB +82)! (V.10) 
Brot = LBrax + Bryy + Brez ) + (Byyy + Brxy + Byez Y° + (Bazz + Baxe + Byyz)° 172 (V.11) 


Tous les identifiants NLO obtenus à partir des équations ci-dessus ont été calculés aux niveaux 
de théorie DFT/6-31G (d,p) pour fournir ses propriétés non linéaires et sont répertoriés dans le 


tableau V.7. 


Tableau V.7 Moment dipolaire (1), polarisabilité moyenne (a: la.u = 0,1482 x 10°?* esu) et la 
première hyperpolarisabilité moyenne (B: 1a.u = 8,6393 x 107% esu) des complexes [(fB-cd)( 


DDQ)] et [(B-CD\TCNE)] par la méthode DFT / 6-31G (d,p). 


Paramètres [(B-CD)DDQ)] [(B-CD)(TCNE)] 
u, D 

Ux 1.1151 -0.9884 
Uy 2.4721 -1.8645 
Hz -4.3297 -0.2459 
u 5.1089 2.1245 
Polarisabilité 

Oxx -532.7838 -508.5661 
Oxy 1.3120 -13.1103 
yy -517.5080 -495.4679 
Oxz -10.4001 -0.6392 
Oyz 6.0024 -7.5920 
Ozz -587.2783 -526.7995 
© (a.u) -545.8567 -510.2778 
a[10"%* esu] -80.8959 -75.6231 
La première hyperpolarisabilité 

Bxxx 27.6487 -29.6212 
Bxxy 22.9674 -71.9374 
Bxyy 71.7486 3.0716 
Byyy 73.4107 30.8195 
Bxxz 25.3357 74.1887 
Bxyz 29.1582 -27.5290 
Byyz 81.4490 60.5808 
Bxzz -41.4961 -57.7181 
Byzz 78.8013 1.0967 
Be -375.7601 -187.7123 
B (au) 326.1719 107.2645 


B [10 esu] 2.8178 0.9266 
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Le moment dipolaire le plus élevé a été déterminé à 5,1089 D pour [(B-CD) (DDQ)|], 
tandis que le plus bas a été prédit à 2,1245 D pour [(B-CD) (TCNE)]. La polarisabilité de [(f- 
CD) (DDQ)] et [(B-CD) (TCNE)] est calculée respectivement à (-80,8959x10* et - 
75,6231x107*) esu. 

L'ampleur de l'hyperpolarisabilité moléculaire $ est l'un des facteurs clés du système 
NLO. La première valeur B d'hyperpolarisabilité statique calculée est égale respectivement à 
2,8178x10* esu et 0,9266x10"* esu pour [(B-CD) (DDQ)] et [(B-CD) (TCNE)] et qui est 
environ 7 et 2 fois supérieures à celles de l'urée (pour l'urée B est de 0,3728x10"* esu) [59]. 
Ainsi, la valeur élevée du moment dipolaire et l'hyperpolarisabilité du premier ordre 


représentent le potentiel NLO du complexe actuel. 


V.3.5 Analyse QTAIM 


Les applications de l'analyse QTAIM (Quantum Theory of Atoms in Molecules) sont 
accrues pour l'étude des interactions intra et inter moléculaires, en particulier les liaisons 
hydrogènes. L'analyse QT AIM a été faite par le logiciel ADF [14]. Les paramètres topologiques 
importants; La densité électronique totale p(r) et son Laplacien V?p(r) ont été obtenus à partir 
de calculs QTAIM afin de donner plus d'informations sur les propriétés des points critiques des 
haisons (BCPSs) [15-17]; notamment les points critiques des liaisons hydrogènes (3,-1), ils 
doivent être positifs et compris entre 0,002 et 0,04 et entre 0,024 et 0,139 pour p(r) et V?p(r) 
respectivement selon la théorie de Bader [15]. 

Les interactions intermoléculaires entre la molécule (B-CD) et les accepteurs (DDQ et 
TCNE) identifiées via une analyse AIM [66] sont représentées sur les figures V.6 et V.7. Les 
valeurs des densités électroniques totales p (r) et son Laplacien V?p(r), de la densité d'énergie 
cinétique et de l'énergie totale aux BCPS sont présentés dans le tableau V.8. 

Comme le montrent le tableau V.8 et la figure V5, 1l existe trois interactions de liaison 
hydrogène différentes dans le complexe [(B-CD) (DDQ)] : O....H, CI....H et N....H. Au BCP 
important entre O155 et H114, la densité de charge p(r) et son laplacien V2p(r) sont 
respectivement de 0,0164 a.u. et 0,0666 a.u. Ce BCP a également la valeur la plus élevée de 
l'énergie totale H(r) = 0,0482 a.u. Deux autres BCPS sont observés: C1160....H128 et 
0154....H79 avec une densité de charge p(r) égale respectivement à 0,0047 a.u. et 0,0142 a.u. 

La présence de cinq BCPS dans [(B-CD) (DDQ)] révèle qu'il existe des interactions 


intermoléculaires entre DDQ et B-CD. 
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Tableau V.8 Paramètres topologiques des complexes [(B-CD) (DDQ)] et [(B-CD) (TCNE)] 
calculés par la théorie quantique des atomes dans les molécules (QTAIM), DET/6- 31G (d,p). 


Complexes BCPs p(r) Vip(r)  V(r) G(r) H(r) -G/V 


DDQ/ B-CD O155...H114  0.0164 0.0666 -0.0162  0.0644  0.0482 3.9753 
CI160....H128 0.0047 00178 -0.0024  0.0034  0.0010 1.4166 
CI160....H134 0.0061 0.0253 -0.0036  0.0050  0.0014 1.3888 
O154....H79 0.0142 0.0647 -0.0153  0.0157  0.0004  1.0261 
N157....H84 0.0131 0.0592 -0.0129 0.013383  0.0009 1.0697 
TCNE/ $-CD N151....H100 0.0070 0.0256 -0.0040  0.0054  0.0014  1.3500 
N153....H123  0.0099 00352 -0.0073  0.0080  0.0007  1.0958 
N157....H82 0.0124 00549 -0.0119 0.012838  0.0009 1.0756 
N157....H80 0.015838 0.0755 -0.0170 0.0179  0.0009 1.0529 
N155....H79 0.0193 0.0978 -0.0223 0.0233  0.0010 1.0448 


p(r): Densité d'électrons à point critique, V2p(r): point critique laplacian, V(r): Point critique Densité d'énergie 
potentielle, G(r): forme lagrangienne de densité d'énergie cinétique, H(r): point critique Énergie totale (toutes les 


unités sont exprimées en a.u). 


Figure V.5 Graphes moléculaires QTAIM des complexes (a): [(B-CD)(DDQ)] et (b):[(B- 
CD\(TCNE)] montrant différents types d'interactions. 


Dans le complexe [(B-CD) (TCNE)], cinq liaisons (N....H) d’hydrogènes 
intermoléculaires ont été observés. La quantité maximale de densité électronique (0,0193 a.u) 
attribuée à l'interaction (N155....H79) est liée au laplacien (V2p(r) = 0,0978 a.u) et (H(r) = 
0,0010 a.u). De même, le BCP lié à l'interaction entre H100 et N151 a également une densité 
de charge plus faible (0,0070 a.u) et un laplacien égal à 0,0256 a.u. avec H(r) = 0,0014 au. 

Comme le montre le tableau V.8, toutes les densités de charge au BCPS sont 


relativement faibles, et leur laplacien est positif ; cela indique que les interactions sont des 
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interactions typiques à coque fermée. 

Les résultats respectifs du NBO et de l'AIM semblent concordants avec respect de la 
nature des principaux sites réactifs et de leurs énergies d'interaction, bien que les bases 
théoriques des deux approches théoriques soient différentes. Les figures V.6 et V.7 illustrent 
une corrélation graphique entre les valeurs d'énergie de perturbation des interactions entre les 
deux orbitales situées aux deux extrémités des chemins de liaison, et les densités électroniques 
des points critiques correspondants. Plus l'énergie de perturbation est grande, plus la densité 
électronique est grande. Par conséquent, les résultats AIM obtenus étaient en accord avec ceux 


de l'analyse NBO. 
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Figure V.6 Corrélation graphique de [(B-CD)(DDQ)] entre les énergies de perturbation (E(?? 


en kcal/mol) et les densités électroniques p(BCP) des points critiques correspondants. 
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Figure V.7 Corrélation graphique de [(B-CD)(TCNE)] entre les énergies de perturbation (E°) 


en kcal/mol) et les densités électroniques p(BCP) des points critiques correspondants. 


V.3.6 Détection des interactions non covalentes 


La technique du gradient de la densité réduite (RDG) peut être exploitée pour détecteur 
les interactions non covalentes développées au sein des complexes de [(B-CD)(DDQ)] et [(B- 


CD)(TCNE)]. 


Ainsi dans sa représentation 2D (Figure V.8), les liaisons hydrogène sont localisés dans 
la partie où sign (A2).p est de -0.020 unités atomiques, les liaisons de van der Waals pour ses 


valeurs nulles et les répulsions stériques pour des valeurs positives. 


24% Partie : Chapitre V Résultats théoriques 


2.000 


1.800 


—_ 


600 


400 


200 


Reduced density gradient 
e 
Oo 
Oo 


0.800 
H-Bond 
0.600 ; 
p 0 
0.400 A2<0 NN 


Fa 
‘'Hhtéragti 


0.000 Le 
-0.050 -0040 -0.030 -0.020 -0010 0.000 0010 0.020 0.030 0.040 0.050 


sign (À\,)p (au) 


Figure V.8 Représentation 2D de l'analyse RDG du complexe [(B-CD)(DDQ)] 


Dans la représentation 3D ( Figure V.9) en analyse RDG visualise par l'intermédiaire des 
couleurs la prédominance des différentes interactions la couleur verte qui entoure la molécule 
de DDQ que les forces de van der waals sont responsables de la stabilisation du complexe de 


transfert de charge 


Figure V.9 Représentation 3D de l'analyse RDG du complexe [(B-CD)(DDQ)] 


Cu 
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GENERALE 


Conclusion générale 


Conclusion Générale 


Dans ce travail, nous avons effectué une étude expérimentale et théorique sur les 
complexes de transfert de charge entre la molécule du f-cyclodextrine en tant que donneur et 
deux molécules en tant qu’accepteurs d’électrons x : DDQ et TCNE. L’étude expérimentale 
confirme la formation des complexes stables de stœchiométrie 1:1. 

Les complexes sous forme solide ont été synthétisés et caractérisés par les techniques 
spectroscopiques IR-TF et RMN du proton. Ces caractérisations confirment l'existence d'une 
interaction par transfert de charge associée à la migration du proton de f-CD vers les accepteurs 
(DDQ et TCNE) pour les complexes [(B-CD) (DDQ)] et [ (B-CD)(TCNE)] montrant des liaisons 
hydrogènes intermoléculaires (N*... H-0O") et (C=0*...H-0). 

Les études de théorie fonctionnelle de la densité (DFT) ont confirmé le mode 
d'interaction entre le fB-CD et les accepteurs (DDQ et TCNE) en comparant les résultats 
théoriques obtenus avec les données expérimentales. 

Une analyse NBO nous a permis, par la suite, de calculer les valeurs des 
énergies de stabilisation E? et d’effectuer une analyse des charges naturelles. Ces 
résultats ont bien confirmé la formation des complexes de transfert de charge. 

D'un autre côté, l'écart énergétique Exomo-Lumo dans chaque complexe a confirmé qu’un 
transfert de charge intermoléculaire a été effectué au sein du système moléculaire 
étudié. 

Enfin, la prédiction des propriétés de liaison s’est appuyée sur deux approches : 
l’analyse topologique QTAIM et l’analyse orbitalaire NBO, celles-ci nous ont permis 
d’expliquer la nature et la rigidité des liaisons hydrogènes notamment altérées dans les 
structures des deux complexes étudiés. 

L’ensemble des calculs théoriques réalisés démontre une parfaite adéquation avec les 
résultats expérimentaux. 

Comme perspectives a ce travail, nous envisageons de compléter ce travail en utilisant: 

> Complexe de transfert de charge entre le a-cyclodextrine et des accepteurs z 
en présence et en absence du solvant. 


> Applications des nouvelles techniques d’analyses computationnelles. 


